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I. INTRODUCTION 

There are two formal ways to deal with the dynamics of General Relativity which we call the Einstein frame and 
the Jordan-Ehlcrs-Kundt frame (JEK frame, for short): 

• The Einstein frame (1915) corresponds to a second order differential equations relating the curvature tensor to 
the energy-momentum tensor; 

• The JEK frame (1960) relates the derivatives of the conformal Weyl tensor to the derivatives of the energy- 
momentum tensor using Bianchi's identities. 

Although the demonstration of the equivalence of JEK formulation to general relativity (GR) have been realized 
by Lichncrowicz in the early 60, its role in the development of applications have been less active than one could 
expect. A possible cause of this is the fact that almost all introductory books on GR do not present the JEK frame as 
an alternative formulation of gravitation. Indeed, only few advanced books — cf. Zakharov (1973), Choquet-Bruhat 
(2009), Hawking-Ellis (1977) — show an overview on this. In particular, as a direct consequence of this is the fact that 
the great majority of analysis of perturbation theory is made ignoring completely the possibility of use JEK frame to 
develop a consistent and worldwide method for this. 

The main goal of the present review is to make a little more popular the use of the JEK frame in the realm of 
perturbation theory of Friedman universes. Indeed, the Lifshitz-Bardeen method and the JEK frame, under the same 
initial conditions, give the same results for the perturbations in the linear regime, as we shall see in Sec. [HID 3]. 
The main interest on JEK rests on its unambiguous way to deal with perturbation within the standard cosmological 
FLRW scenario. 

In this paper we concentrate our attention to the (cosmological) perturbation scheme presented in Sec. [Ill], al- 
though we will describe some examples of well-known solutions — Schwarzschild, Kasner, Friedman (singular and 
nonsingular) — according to the JEK frame (cf. Section [II]) to show how this method could be used to obtain new 
solutions of general relativity (details of this discussion in Sec. [IV]). 

Finally, the JEK-framc is called alternatively the Quasi-Maxwcllian version of general relativity. The reason for this 
name is manifest from its strike similitude to Maxwell's equations of electrodynamics. In the appendix, this similarity 
is used to exhibit an example of modification of general relativity by extending further such analogy to the case of 
the electrodynamics inside a dielectric medium. 

A. Definitions, Notations and a Brief Mathematical Compendium 

We shall list below all the definitions that should be used in this review: 

• The structure of the space-time is represented by a Riemannian geometry g^ v {x a ), with Lorentzian signature 

(+,-->-); 

• The Levi-Civita tensor rj^ a p = y/—g where g is the determinant of g^ v and e^ap is the completely 
antisymmetric pseudo-tensor, 60123 = 1; 

• The Christoffel symbols arc defined by = \g aX {gp\^ + g^\,p - gpn,\); 

• The geodesic equation is 

d 2 x» M dx a dx? 
dX 2 a ? dX dX ' 

• The Riemann tensor is defined by 

pa pa pa i pa ft pa pr . 

The traces of Riemann tensor define the Ricci tensor = R a and the curvature scalar R = R a a ; 

• The decomposition of the energy-momentum tensor into irreducible parts, with respect to a normalized observer 
field V a , is given by 

T^ u = pV^Vy - ph^ + V^q v ) + -k^v, 

where p is the energy density, p is the isotropic pressure, is the heat flux and 7r M „ is the anisotropic pressure. 
We use parentheses '()' for symmctrization and brackets '[]' for skew-symmctrization; 



The Wcyl tensor is defined by 



where 
and 

g a 0tiv = g a fi.g0u — gavgp^- 

The dual is denoted by 

Remark that W^ v = W a ^; 
• The electric and magnetic parts of the Weyl tensor are, respectively, 

E aP = -W aflf}v V»V v 

and 



• The tensor defined by h^ v = g^ v — V^V;, projects tensorial quantities in the rest space E of the observers. Note 
that h flu V v = and h flu h u \ = h^\. 



Einstein's equations (EE) are given by 



Rjxv ~ ~^Rg^u + ^g^u — ~kT lxv , 



where A is the cosmological constant and k = 9>ttGn/c a which we shall set equal to 1, unless stated otherwise; 
The covariant derivative of can be decomposed into its irreducible parts, that is, 



where = V a - a is the expansion coefficient, 



f 6 

0>„ = -jh^h^V^.p) - g^Atf 



is the shear tensor, 

f 



t*V = 2 h %hu V [a;i3] 

is the vorticity and a M is the acceleration; 
One can use the quantities defined above to obtain the evolution equations of the kinematical quantities: 
• Raychaudhuri equation 

e + L. + 2 a 2 + 2uj 2 -a»., fl = ~(p + 3p), (f) 
where 2cr 2 = a^a^ and 2uj 2 = lj^lj^ and X = X_ a V a (this last definition will be used throughout the text); 



• The evolution of the shear tensor is 

h^h^&au + \h af 3{a x .,\ - 2a 2 - 2oj 2 ) + a a a p - Xh^h^ '(a„.„ + a„ ;M ) + 

(2) 

2 1 

+ Tj6>(7 a /3 + CTa^a^p + UJa^UJ^ {} = —E a/3 — ^7T Q/J ; 

• The evolution equation for the vorticity tensor is given by 

1 2 

h^h^Co^ - -h a ^hf} y {a Kl , - a u -^) + -9u> aj3 - ap^a + a a ^p = 0; (3) 

• These kinematical quantities must satisfy three constraint equations: 

^e^x - K 7 + w Q 7 );a/i 7 A - a" (ax, + luxv) = -qx; (4) 

uj a , a + 2cu a a a = 0, (5) 

where io a = — \rf l ^ s uip 1 Vs and 

- ^h ir e h X ) a r]/' r ' / V 1 ,(a a p + w Q /j) ;7 + a (r cj A) = H rX ] (6) 

• The conservation law of the energy-momentum tensor expressed in terms of its components is the conservation 
equation 

p + ( P + p)9 + q^V/j, + g% - n^a^ = 0, (7) 
and the generalized Euler equation 

(p + p)a a - p^h^a + q^a + 6q a + q v Q av + q v u av + -K a v - v + ^ v o ^ v V a = 0. (8) 

These formulas summarize the kinematical part of the QM-approach. Now we shall focus on the dynamical equations 
in terms of the Weyl tensor. 

B. Quasi-Maxwellian Equations 

The quasi-Maxwellian equations [62, 76, 118, 136] are obtained from Bianchi's identities written in terms of the 
Weyl tensor, i.e. 

W afStlv . v = -R^ a ^ - — g^ a R>®. (9) 
2 12 

Substituting the Einstein equations, we get 

W a ^^ v - v — — -T^ a ' }13 ^ + -g^ a T'P\ (10) 

From practical analysis, it is worth to project these equations with respect to a vector field V a and its orthogonal 
hyper-surface of spatial metric h^ u . There are four possibilities to do such decomposition and, therefore, four linearly 
independent equations: 

• Projection VpVph a a gives 

h ea h^E aX ;-y + v e 0^VPH vX a»x + 'iH™u y = ltf a p, a + (Lq* + 
1, 1 1 11 
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(12) 



Projection r)' TX a pV fl V\ yields 

h £a h x ^H aX;l - rf^ v V^E vX G^ x - 3E"u v = (j> + p)ui* - \il^ x V x q a , p + 

Projection h^ a 7] T ^ x a pVx gives 

h^h x H^ + 6H tX - \H v ^h x \V^ v + v Xl/ ^if pra V t ,V T H ai e^+ 

-a„£/5 ( V )7a!/ % + ^^^ ;Q 4V )7Q/3 ^7 = -|<? (e w A ) + ^ £ V^+ ( 13 ) 

Projection Vph^ T h^ a yields 

h ll e h u x W v + 6E eX - ^E^h^^V^ + v Xl/tll V tpTa V^V T E a7 9 uP + 

-\{p + p)(J cX + \q {t a^ - \hf*'h x > a q ll . a + \h^h^ X TT a ^ + \n^a x ^ + 

-l^uP* + le^ x . 



(14) 



Eqs. (11)-(14) are the QM equations. Now let us show the consistence of the QM-formalism and its equivalence to 
the dynamics of general relativity. 

C. Equivalence between QM Equations and GR 

The QM formalism is supported by the theorems that we shall mention in this section. Although our quoted 
references are mostly interested into its formal aspects, we shall focus here on the physical results. 

Following the same steps given by Lichnerowicz (1960) in order to prove the equivalence between QM equations 
and GR, we start considering a manifold M with n + 1 dimensions endowed with a hyperbolic metric satisfying 
Einstein's equations. Suppose the existence of a hyper-surface E with local equation f(x a ) = 0. One assumes that 
the discontinuity of the derivatives of when it crosses E is provided by Hadamard's conditions, i.e. 

\giw,a,p]s = a^f, a f.j3, (15) 
where the amplitude of the discontinuities a^„, under local coordinates, transforms as 

where = f.^. The discontinuity of the Jacobian matrix J£ is defined by 

=*"W 7 , (17) 

where t a is an arbitrary vector. 

As a consequence, it follows that the discontinuity relations for the Ricmann tensor is 

[RaPnv\z = -^{a afl lpl\ + ap\lalfi — a a \lpl» ~ a /3^al\), (18) 

and for the Ricci tensor is 

[Ra[i]Y, = -^9 pa {a ap ljil a + ap a l a l p — a pa lpl a — a a fjl p l a ). (19) 
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The validity of Einstein's equations for an empty space-time on E implies that [i? M „]s = 0, if and only if, the null 
vector l a is an eigenvector of the matrix a M „, that is 

a a pl p = \l w (20) 

where a = g^ u a^ v . This result shows that the coefficients of discontinuity are not arbitrary. Now let us analyze the 
discontinuity relation imposed by the Bianchi identity for the Riemann tensor. First of all, consider f(x a ) = as 
being the local equation of the hyper-surface S. From the definition of l a , it has null vorticity, i.e., l[ a -p\ = 0. The 
cyclic identity applied to the discontinuity of the Riemann tensor implies that 

]=0, (21) 

and 



l P [R p ^}=0. (22) 



The covariant derivative of Eq. (21) yields 



{lp[Ra/3\p} + l\ [RaP/ip] + lp,[RaPp\]);v = 0. (23) 

Making the contraction between the indices p and v, and then, considering that the Einstein equations for an empty 
space-time (R^v = 0) hold on S, we obtain 

2l P [R a pAp.];p + l P ;p[R al 3Xp]=0. (24) 

It means that if [Rapxp] vanishes at some point x of S, then it vanishes along the whole isotropic geodesic passing by 
x. 

In the general case, where l a is considered either space-type or null-like, we have the following result: Let Q, be 
an oriented hyper-surface in the space-time intersecting T,(x° = 0) and defining a 2-surfacc U . If one gives Cauchy's 
data g^ v and Qpv,\ on O, such that crossing on E the second derivatives admit the discontinuities [g a /3,oo] = (d a /3)s, 
then a a p on U must satisfy the condition 



(<V,-|<^)r| y = 0. (25) 



Equivalently, if one gives for all points x of U the tensor [R a ppv\u admitting as fundamental vector (l^u-, which 
contracted to the Riemann tensor is zero, then the Cauchy data considered correspond to the solution of Einstein's 
equations such that the curvature tensor, when crossing S, admits a discontinuity [R a ppv\- The tensor [Rappv] is 
necessarily the solution of Eq. (24) corresponding to the initial data [R a ppv\u- 

The results above for non-empty Einstein's equations = —kT^ are easily proven due to the continuity of T Mi , 
through £ and can be seen in Lichncrowicz, 1960 or Novello & Salim, 1985. From these considerations we can state 
the following lemma: 

Lemma (Lichnerowicz). Bianchi 's identity together with the convenient Cauchy data represented by Eq. (25) are 
equivalent to Einstein's equations. 

This is the main result that will be used in this review. Next we apply this method of dealing with Quasi-Maxwellian 
formalism by analyzing some special solutions of the equations of general relativity. 

II. PARTICULAR SOLUTIONS OF GR FROM QM EQUATIONS 

In order to present specific examples of how QM equations work to, we reproduce some known solutions of the 
general relativity theory using the quasi-Maxwcllian framework. Our task is simplified if we use Gaussian coordinates, 
where g^^ = 5°^ and, moreover, a foliation described by the observer field = 5$. This coordinate system will be 
used to deduce all solutions presented here. 
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A. Schwarzschild Solution 



The Schwarzschild metric in a Gaussian coordinate system takes the form: 



ds 2 = dT 2 - B(T, R)dR 2 - r 2 (T, R)dfl 2 . 
A geodesic observer in this metric is = Sq. The expansion coefficient 9 is given by 



(26) 



1 B 4r 

2 l B + V 



(27) 



where Y(T, R) = dY/dT. Using the metric (26), we write the corresponding shear tensor a^ v and the electric part of 
Weyl tensor E^ v in the matrix form: 



[a i j ]=f(T,R) 



( 1 

-I 







1 

"2 



(28) 



1 
1 n 



[E* j ]=g(T,R) | - 2 




1 

7 



(29) 



where 



and 



12 ^ i? ) = - 2 f + 5- 4 S + 2 ;f + 2 7^ + 4 ;-^- 4 ^ + 4 W 



where Y'(T,R) = dY/dR. The magnetic part of Weyl tensor i? a ^, the vorticity w a/ 3 and the acceleration are 
identically zero. The set of quasi-Maxwellian equations (11)-(14) reduces to the form 



3— g = 0, (32a) 



r 



r 



g + 3-g = (32b) 



r 



The evolution equations of the remaining kinematical quantities are provided by the Raychaudhuri equation and 
the shear evolution 



^+y + ^/ 2 =0, (33a) 

f + Y + l ef = ~ g - (33b) 



Finally, the only non-trivial remaining constraint equation is 



/' + 3-/- V = 0. 
r 3 



(34) 
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This is nothing but the Schwarzschild solution in Gaussian coordinates. Indeed, functions B(T,R) and r(T,R) are 
obtained by imposing on the Cauchy surface T = Tq the Lichnerowicz condition 



R^v — 



(35) 

To 



Then it follows 



^(T,R) = TTm , (36a) 

r E (T, R) = - y/F + r H /r, (36b) 
r' E (T, R) = w'(R)^F + r H /r, (36c) 

where F(R) and w(R) are arbitrary functions and r# is an arbitrary constant 1 . These expressions play the role of 
initial conditions on the Cauchy surface such that the functions Bqm and Tqm must be equal to Be and te on To, 
and then, they are evolved to the whole space-time. 
From Eqs. (32), it follows 

.9 = -4' ( 37 ) 
r 

where k is another arbitrary constant. From Eq. (34), we have 

Bqm = TTW (38) 

where h(R) is an arbitrary function. One can write Eqs. (33) in terms of the functions B and r, as follows 



3g, (39a) 

(39b) 

Substituting Eq. (38) in (39b) yields 

r' r' 

- + 2- = (40) 
r r 

which can be integrated with respect to R resulting in 

(41) 

r 

where a(T) is an arbitrary function of T. Using the subtraction between (39a) and (39b), and then, manipulating 
Eqs. (38) and (37), it follows that a(T) is a constant (a(T) = -k/2). From Eq. (41), we get 

k 



B 


1 B 2 


-2t = - 
r 


B 




i) 


IB 2 


r 




2B^ 


+ 4- = 0. 


B 


r 



r 2 =2y(R) + -, (42) 



where y(R) is an arbitrary function. Substituting the constraint (38) in the definition of g — Eq. (31) — and again 
taking the subtraction between Eqs. (39), we obtain 



We introduce the subscript indexes E and QM to distinguish the solution of the Einstein equations valid only on the Cauchy surface 
and the solution propagated by the Quasi-Maxwellian equations, respectively. 
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Substituting (42) into (43) yields 



r r r' 1 h 1 r 2 h 
3 7+o— + — ~ — = 0. 



A _ M - . it\ 

2 r "2 X ^ /> 



where x(T) is an arbitrary function of T. Eq. (42) gives 



dr 



y/h + k/r 



= / dT, 



which we integrate to get 



(43) 



(44) 



(45) 



\J (hr + k)r 



h 2h 3 ^ 2 
Partially deriving this equation with respect to R, yields 



In (k + 2hr + 2^/{hr + k)hr) = T + b(R). 



(46) 



' QM 



= b\ h + -. 



We still must determine the arbitrary functions in (36). In order to do this, we set 



(47) 



F = a 2 - 1, w(R) = aR, 
where a is a constant parameter. Then, considering the Jacobian matrix J a p 



(48) 



[J" 



dx a 



( a/ A -(a 2 - A) /A \ 



V 








.4 



aVa 2 - A 



1 
I J 



(49) 



where x a = (t, r, 9, (f>), x 13 = (T, i?, 8, <f>) and A = 1 — Th/i", wc map the metric given in Gaussian coordinates by (26) 
into the well-known Schwarzschild solution in the usual Schwarzschild coordinates 



ds 2 



1 - r JL ) dt 2 



(50) 



if wc make the identification r# = 2MG/c 2 . The parameter a is interpreted as the mechanical energy of a test particle 
obtained by integration of the geodesic equations. 

Due to the similar functional form of Eqs. (38), (42) and (47) in comparison to Eqs. (36), the determination of the 
arbitrary functions on the hyper-surface Tq is trivial: 



Be(T ,R) = Bqm{Tq,R) 
r E (T , R) = r QM (T 0l R) = 
r' E (T ,R) = r' QM (T ,R) = 



» h(R) = F(R), 
k = k\, 
b{R) = w{R). 



(51a) 
(51b) 
(51c) 



Therefore, the metric (26) becomes 



ds 2 = dT 2 



a 2 -I 



2M 
r(T,R) 



dR 2 -r 2 (T, R)dn 2 



(52) 
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The Schwarzschild internal solution is more involved to be obtained from the above steps, because the Cauchy 
surface Tq for this case is different from the spherically symmetric surface r — r$ — used in Schwarzschild coordinates 
to apply the match conditions between the internal and the external parts. Besides, Gaussian observers do not 
decompose the energy- momentum tensor as a perfect fluid contrary to the Schwarzschild observers = y/goo^Q, 
which turn calculations more cumbersome. Indeed, the energy-momentum tensor associated to the Gaussian observers 



6% is 



Tl? = (PG + Pg)V»V» - p G g^ + V {ll q v) + n^. (53) 
The energy-momentum tensor decomposed in terms of the observer field u„ = e^ v ^ T,R '{a, 1,0,0) is given by 

V = (P + P)«V-pv. (54) 

where v = v{T, R) must satisfy the Tolman-Oppcnheimer-Volkov equation and a is the external parameter associated 
to the co-moving test particle. The physical properties of the fluid in both representations are linked by the following 
expressions 

PG = {p + p)a 2 e~ v -p, (55a) 
PG = -\[{p + v){l-o?e-»)-2 l pl (55b) 
q i = (p + p)ae-"S i 1 , (55c) 

tt%- = H(i - a 2 e-n diag(l, -1/2, -1/2). (55d) 

Substituting these equations in the quasi-Maxwellian equations, one exactly find the Schwarzschild internal solution 
containing some arbitrary functions. Match conditions must be used in order that this solution is joined to Einstein's 
equations on the hyper-surface. Choosing the Cauchy surface T = T$, one fixes all arbitrary functions obtaining the 
Schwarzschild stellar solution with such procedure. 

B. Kasner Solution 

Different from the case of Schwarzschild metric, we shall use the Hadamard method in order to get the Kasner 
solution. We set for the Bianchi-I anisotropic metric the form 

ds 2 = dt 2 - a 2 (t)dx 2 - b 2 (t)dy 2 - c 2 {t)dz 2 . (56) 

In this case there is solely one non-trivial QM-cquation: the "time" evolution of the electric part of the Weyl tensor 
that reads 

p\ + q E eX 3 a ^e E X)^ + h eX^ v flu ^ = ^ ( 57 ) 

where we consider an observer field = Sq and = 0. The non-trivial equations of the kincmatical quantities are 

9+ — + 2(T 2 = 0, (58) 



and 



2 2 2 

cr a p + E a p - -h a pa + -9<r al3 + a a ^p = 0. (59) 
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To obtain the Kasner solution, we set 

a(t) = tP\ b(t) = c(t) = t? 3 , (60) 

where p\ , p2 and p% must satisfy 

Pi +P2 +P3 = 1, 

(61) 

(Pi) 2 + {P2? + (P3? = 1. 

It is necessary to consider the Kasner solution on a Cauchy hyper-surface t = to as initial condition for the Eqs. 
(57)-(59). For this case, one can reduce time evolution of E^. apv and in some algebraic constraints, that is: 

W v = -26E» V , 

cr% = -e<j* v , (62) 
6 = -2d 2 . 

Note that these expressions arc valid only on the Cauchy surface t = to. Using relations (62) in Eqs. (57), (58) and 
(59), these equations become just algebraic expressions for the Kasner background. As a consequence, one can define 
three special variables for the Kasner background: 

X a /3 = E a p — ^h a pa 2 — ^0a a/3 + cr QAl cr M ( 3, 

Y aj 3 = 9E a p + |cr''( Q ,£ , / 3) AI — h af) a iL v E v M , (63) 
W = 2cr 2 - |6> 2 . 

These quantities are identically zero for the Kasner solution (on to) and they represent the minimal set of variables 
which contains all the information about such metric, because they come from the nontrivial equations of the QM- 
formalism. Once X a p, Y a p and W are zero on to, the QM-equation will propagate these quantities to the hyper-surface 
in the vicinity of to retaining their null values, due to their tensorial features. In other words, we obtain the validity 
of Kasner solution for the whole space-time, according to the theorems of Sec. [IC]. 



C. Friedman Solution 



Consider the isotropic metric given in the Gaussian coordinate system: 



ds 2 = dt 2 + gij dx l dx 3 = dt 2 - a 2 (t)[d X 2 - cr 2 {x)dn 2 } 1 (64) 

in which gij = —a 2 (t)'fij(x k ). The material content of this universe is represented by a perfect fluid, with energy 
density p, pressure p and equation of state p = Xp, where A is a constant. 

A direct inspection of Eqs. (11) - (14) shows that the quasi-Maxwellian equations for this metric are identically 
zero, because the metric is conformally flat (we use an observer field — S^). The only kinematical surviving 
equation is the Raychaudhuri equation 



\d 2 = ~(l + 3A)p. (65) 



The energy-momentum conservation reduces to 



p+(p + p)6 = 0, 
h a "p, M = 0,=^p = p(t). 



(66a) 
(66b) 
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Thus, 



and 



3- = ~(l+3A)p, (67) 
a 2 



p + 3(l + A)p- =0. (68) 
a 



Eq. (68) can be integrated to yield 

p = Po a- 3 ( 1+A \ (69) 
where po is a constant. Substituing (69) in (65), we obtain the Friedman equation 

a 2 e I ,. 
a a S 

where e is a constant of integration. 

Finally, we can exhibit the scale factor a(t) in terms of a quadrature equation 



/ 



da 1 . . .„ , 

-.(t-t ), (71) 



Vpoa- (1+3A) +3e 73 



where io is a constant of integration. 

The initial conditions necessary to solve this problem are a(t), a(t), a(t) and a(r) on the Cauchy surface. On the 
other hand, we have A, e, po and to- Instead of specifying each initial condition of the Cauchy problem, one can 
cquivalcntly fix each free parameter. It can be done if we write the Ricmann, Ricci and curvature tensors in terms of 
the 3-geometry of the background h^, as: 

Raf3^u = % {hun hp v — h av hfj/j), 

R - 2e a . 

n-pv — 1 "<8v> 

a a 

where we use the following relation that holds for the 3-geometry: 

V ' a V /3 Xry — Vp V Q Xj = —R X 1 f}a X\, 

The symbol (") means a projection on the hyper-surface defined by . 

The explicit expression of R is obtained from the Friedman metric as follows 

- 6e = Ra 2 = -4— + - 2—^. (72) 
a cr z a 

The only three possible solutions for this equation are listed in Table I, which joins the solutions of Friedmann 
equation for different values of A and e. The constant ao, which is written in terms of po and to, has different values 
for each solution in the table and it is commonly interpreted as the "current size of the Universe" . 



D. Nonsingular Solutions 

There are many proposals of cosmological solutions without a primordial singularity. Such models are based 
on a variety of distinct mechanisms, such as cosmological constant, non-minimal couplings, nonlinear Lagrangians 
involving quadratic terms in the curvature, modifications of the geometric structure of space-time, non-equilibrium 
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p 


A 


e 


6 


a(n) 


tin) 


4,-2 

3* 












V 


3 + -2 


1/3 





3,-i 
T 


<W /3 


V 


4(1 -COST,)" 3 





1 


3 sin 77 
a o (1 - cos?;) 2 


ao(l — cos 77) 


a (ij — sin 77) 


3 1 

2 • 4 

an sm ?7 


1/3 


1 


3 cos 77 
a'o sin 2 7? 


ao sin rj 


ao(l — cos 77) 


-4 (cosh n- 1)~ 3 





-1 


3 sinh 
a ° (cosh 77 - l) 2 


ao(cosh 77 — 1) 


ao (sinh 77 — rj) 


3 1 
a 2 , sinh 4 77 


1/3 


-1 


3 cosh 77 
a 2 , sinh 2 77 


ao sinh rj 


ao(cosh 77 — 1) 



TABLE I: Fundamental quantities of Friedman Universe. (Units system k=c=l.) 



thermodynamics, among others cf. dc Sitter (1917), Murphy (1973), Novcllo and Salim (1979), Salim and de Olivcra 
(1988), Mukhanov and Brandcnberger (1992), Mukhanov and Sornborger (1993), Novello et al. (1993), Moessner and 
Trodden (1995), Saa et al. (2000). Recently, an inhomogeneous and anisotropic nonsingular model for the early 
universe filled with a Born-Infcld-type nonlinear electromagnetic field was presented by Garcia-Salccdo and Breton 
(2000). Further investigations on regular cosmological solutions can be found in Klippcrt et al. (2000), Vcncziano 
(2000) or Acacio dc Barros et al. (1998). A complete listing of nonsingular solutions can be seen at Novello and 
Bergliaffa (2008). Here wc shall analyze some of these examples presented in the literature — cf. Rcfs. [29] and [38], 
for instance — using the Quasi-Maxwellian formalism. 

1. A WIST Model 

In the Weyl integrablc spacetime model (WIST) — cf. Novello et al. (1993), Salim and Sautu (1996) and Fabris et 
al. (1998) — as well as in string theory (Gasperini, 2003), there are models that describe the geometry coupled 
to a scalar field. In those models, there are nonsingular solutions for an FLRW geometry. In order to search for a 
simple bounce scenario in cosmology, described by an analytical exact solution, we fix our attention on the background 
discussed by Novello et al. (1993), Salim and Sautu (1996) and Oliveira et al. (1997). Basically, the model concerns 
a modified Riemannian geometry with metric g^u and an extra Weyl affinity given by 

Kp = {Ze} + l(% u ,e + 6%u>, a -g aflU >>>). 

In the Weyl manifold, the vacuum field equations can be rewritten in terms of a Riemannian geometry plus a term 
dependent of the Weyl field u> included in the non-metric part of the affinity. At this level, the field equations can 
be represented by a perfect fluid with the 4- velocity given by V M = d^uj/ui 2 , where u> 2 = g a ^uj^ a uj^ and equation of 
state p = p. Originally, the scalar field is part of the affinity. However, it is transposed to the right side of the field 
equations and it can be interpreted as a perfect fluid. In this case, its effective energy density appears as a negative 
quantity. The quasi-Maxwellian equations of motion of the background written in the conformal time are 



(a') 2 + ea 2 + — (w'a) 2 = 0, (73) 
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and 

w' = ja-' 2 , (74) 

where 7 is a constant and A 2 is the coupling constant between the scalar field and the metric tensor. It follows from 
these equations that 

(a') 2 = -«> 2 -4> (75) 
a 

where we defined a 2 , = Xj/\/6. Only solutions with three curvature e = — 1 are possible. The scale factor, solution to 
(73), is given by 

a(rj) = a oV /cosh(2? ? + 5), (76) 

where 5 is a constant of integration. The scalar factor displays a bounce produced by the scalar field that was 
introduced as the Weyl part of the affinity, due to the nonmctricity condition. 

2. Nonsingular Solution From Nonlinear Electrodynamics 

The standard cosmological model, based on Friedman-Lemaitrc- Robertson- Walker (FLRW) geometry with Maxwell 
electrodynamics as its source, leads to a cosmological singularity at a finite time in the past as seen in Sec. [II C], 
Such a mathematical singularity itself shows that, around the point of maximum condensation, the curvature and 
the energy density are arbitrarily large, thus being beyond the domain of applicability of the model. This difficulty 
raises also secondary problems, such as the horizon problem: the Universe seems to be very homogeneous over scales 
which approach its causally correlated region, as pointed out by Brandenberger (1996). These secondary problems 
are usually solved by introducing geometric scalar fields (for a review on this approach see Kofman et al. (1997) and 
references therein). 

This section shall present that homogeneous and isotropic nonsingular FLRW solutions that are obtained by con- 
sidering a toy model generalization of Maxwell electrodynamics. Here it is presented as a local covariant and gauge- 
invariant Lagrangian which depends on the field invariants up to the second order, as a source of classical Einstein 
equations. This modification is expected to be relevant when the fields reach high values, as occurs in the primeval era 
of the Universe. Consequences of the inevitability of the singularity through the singularity theorems (see. Hawking 
and Ellis, 1973) are circumvented by the appearance of a high (nevertheless finite) negative pressure in the early 
phase of FLRW geometry. The influence of other kinds of matter on the evolution of the universe were also taken 
into account. It is shown that standard matter, even in its ultra-relativistic state, is unable to modify the regularity 
of the obtained solution. 

Heaviside non-rationalized units are used. The volumetric spatial average of an arbitrary quantity X for a given 
instant of time t is defined as 

( X ) = J™ T7 [V^gd^X, (77) 



v^v a V _ 

where V = J y/—gd 3 x l , and V stands for the time dependent volume of the whole space. An extended discussion 
about averages in cosmological models can be seen in Rcfs. [14, 53, 155, 156]. 

Average Process: since the spatial sections of FLRW geometry are isotropic, electromagnetic fields can generate 
such a universe only if an averaging procedure is performed — cf. Tolman and Ehrenfcst (1930), Hindmarsh and 
Everett (1998). The standard way to do this is just to set the following mean values for the electric Ei and magnetic 
H fields: 

(Ei) = 0, (Hi) = 0, (EiHj) = 0, (78) 



(EiEj) = ~E 2 9ij , (79) 
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(HiHj) = --H 2 9ij . (80) 
The energy-momentum tensor associated with Maxwell Lagrangian is given by 



Tpv — F^Fau + -Fg^, (81) 

in which F = F^F^ U = 2(H 2 — E 2 ). Using the above average values, it follows that Eq. (81) reduces to a perfect 
fluid configuration with energy density p 1 and pressure p 7 as 

(T M „) = (p 7 +p 1 )V fl V v -p 1 g^ v , (82) 

where 



Pl =3 Pl = -(E 2 +H 2 ). (83) 

From the Raychaudhuri equation, we can see that the singular nature of FLRW universes comes from the fact that 
both the energy density and the pressure are positive definite for all time. Thus Einstein equations for the above 
energy-momentum configuration yield 



o(t) = y/a^t-et 2 , (84) 

where a is an arbitrary constant. 

Nonsingular FLRW Universes: Nonlinear generalization of Maxwell electromagnetic Lagrangian will be con- 
sidered up to second order terms in the field invariants F and G = \rj a fi llv F a ^ F^ = —A(E ■ H) as 



L= --F + aF 2 + f3G 2 , (85) 

where a and /3 are arbitrary constants 2 . Maxwell electrodynamics can be formally obtained from Eq. (85) by setting 
a = /3 = 0. Alternatively, it can also be dynamically obtained from the nonlinear theory in the limit of small fields. 
The energy-momentum tensor for arbitrary nonlinear electromagnetic theories reads 



= —ALpF^Foaj + (GL G - L)g^, (86) 

in which L p represents the partial derivative of the Lagrangian with respect to the invariant F and similarly for the 
invariant G. In the linear case, expression (86) reduces to the usual form (81). 

Since we are mainly interested in the analysis of the behavior of this system in the early universe, where matter 
should be identified with a primordial plasma - for instance, Tajima et al. (1992), Giovannini and Shaposhnikov 
(1998), and Campos and Hu (1998) — we are led to limit our considerations to the case in which only the average 
of the squared magnetic field H 2 survives as was also done by Tajima et al. (1992), Dunne (1997), Joyce and 
Shaposhnikov (1997), Giovannini and Shaposhnikov (1998), Dunne and Hall (1998). This is formally equivalent to 
put E 2 = in Eq. (79), and physically means to neglect bulk viscosity terms in the electric conductivity of the 
primordial plasma. 

The homogeneous Lagrangian (85) requires some spatial averages over large scales, as given by Eqs. (78) (80). If 
one intends to make similar calculations on smaller scales, then either more involved non homogeneous Lagrangians 
should be used or some additional magneto-hydrodynamical effect introduced, as were done by Thompson and Blacs 



2 If we consider that the origin of these corrections comes from quantum fluctuations, then the value of the constants a and /9 are fixed — see 
Heisenberg and Euler (1936). 
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(1998), and Subramanian and Barrow (1998), should be devised in order to achieve correlation at the desired scale 
(see Jedamzik ct al., 1998). Since the average procedure is independent from the equations of the electromagnetic 
field, we can use the above formulas (78) (80) to arrive at a counterpart of expression (82) for the non-Maxwellian case. 
The average energy-momentum tensor is identified as a perfect fluid (82) with modified expressions for the energy 
density p 7 and pressure p 1 as 



and 



Pl = l -H 2 (l-8aH 2 ), 



p 7 = -H 2 (l- AQaH 2 ). 



(87) 



Inserting expressions (87) and (88) in the continuity equation for a Friedman model (66a) yields 



H 

H = — , 



(89) 



where H a is a constant. With this result, a similar procedure applied to Eq. (70) leads to 



k_Hl 

6a 2 



1 - 



8am 



(90) 



where k is the Einstein constant. As far as the right-hand side of Eq. (90) must not be negative, it follows that, 
regardless of the value of e, for a > the scale factor a(t) cannot be arbitrarily small. 
The solution of Eq. (90) is implicitly given as 



ct = ± 



a(t) 



dz 



urn 

6z 2 



6z 6 



where a(0) = a . The linear case (84) can be achieved from Eq. (91) by setting a = 0. 
A closed form of Eq. (91) for e = ±1 can be derived as 



ct 



(xi - x 3 )£ (arcsin y'^zfj, J§H§) + x ^ ( arcsin 



Xl 



X-2 — XI ■ 



Xl - X 2 
Xl - X3 



y/Xl - X 3 

where xi, X2, £3 are the three roots of the equation 8akH^ — kH 2 x + 3ex 3 = 0, and 



z = a?(t) 



(91) 



(92) 



J-(x, k) 



dz 



^J{1- z 2 ){\- n 2 z 2 )' 



£(x, k) 



1 - K 2 Z 2 

1-z 2 



(93) 



are the elliptic functions of the first and second kinds, respectively (see expressions 8.111.2 and 8.111.3 in Gradshtcyn 
and Ryzhik, 1965). The behavior of a(t) for e = ±1 is displayed in the Fig. (1). 

For the Euclidean section, by suitably choosing the origin of time, expression (91) can be solved as 



a 2 = H D ^ -(fcc 2 t 2 + 12a). 
From Eq. (89), the average strength of the magnetic field H evolves in time as 



(94) 



H 



1 



2 kc 2 t 2 + 12a 



(95) 
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FIG. 1: Plots from Eq. (91). We set a(l) = 1, kH„ = 12, and aH 2 = (0; ±1, 25.10" 4 ) as illustrative values. 



Expression (94) is singular for a < 0, as there is a time t = \J — 12a/kc 2 for which a(t) is arbitrarily small. 



Otherwise, for a > we recognize that at t = the radius of the universe attains a minimum value a„ 
given by 



which is 



a min ~ H \/8a. 



(96) 



Therefore, the actual value of a m i„ depends on H Q , which turns out to be the only free parameter of the present 
model. The energy density p 7 given by Eq. (87) reaches its maximum value p ma x = 1/64q at the instant t = t c , where 



12a 



(97) 



For smaller values of t the energy density decreases, vanishing at t = 0, while the pressure becomes negative. Only 
for times t < 10y/a/kc 2 the nonlinear effects are relevant for cosmological solution of the normalized scale factor, as 
shown in Fig. 2. Indeed, solution (94) fits the standard expression (84) of the Maxwell case at the limit of large times. 

The energy- momentum tensor (86) is not trace- free for a ^ 0. Thus, the equation of state p 7 = p 7 (p 7 ) is no 
longer given by the Maxwellian value; it has instead a quintessential-like term — see Caldwell et al. (1998) — which is 
proportional to the constant a. That is 
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aH 4 . 



Equation (98) can also be written in the form 



(98) 



(99) 



where Q(z) is the Heaviside step function. The right-hand side of Eq. (99) behaves as (1 — 64ap 7 )p 7 /3 for t > t c in 
the Maxwell limit ap 1 <C 1. 

The maximum temperature corresponding to t = t c is given by 



I c y/ 4 

\2Aaa) 



(100) 
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FIG. 2: On the left panel: time dependence of the electromagnetic energy density p 7 and pressure p.,. p max = l/64a and t c 
is given by Eq. (97). On the right panel: nonsingular behavior of the scale factor a(t). a m i n and t c are given from Eqs. (96), 
(97). The corresponding classical expression (84) is shown (dashed line) for comparison, with a = a min . 

where a is the Stefan-Boltzmann constant. 

Therefore, the consequences of the minimal coupling of gravity with second order nonlinear electrodynamics, from 
the cosmological point of view, propose relevant modification only in the primeval era of the universe. Indeed, the 
class of theories a > leads to nonsingular solutions for which the scale factor a(t) attains a minimum value. The 
regularity of this cosmological solution is to be attributed to the fact that, for some interval of time, the quantity 
p + 3p becomes negative. 

III. PERTURBATION THEORY IN QM FORMALISM 

Since the original paper of Lifshitz and Khalatnikov (1963), it has been a common practice to start the examination 
of the perturbation theory of Einstein's equations of general relativity by considering variations of non observable 
quantities such as 5g^ v . However, the main drawback of this procedure is that it mixes true perturbations and arbitrary 
coordinate transformations. A solution for this difficulty was found by looking for gauge-independent combinations 
which are written in terms of the metric tensor and its derivatives by many authors (cf. Bardeen 1980, Hawking 1966, 
Jones 1976, Olson 1976, Brandenberger 1983, Vishniac 1990, and Mukhanov 1992). 

The fundamental element of the gauge problem in the perturbation theory of RG was clear and geometrically 
detailed by Stewart's Lemma (1974, 1990): gauge invariant variables (scalars or not) are those which are identically 
null on the background. After that, in Stewart sense, Hawking (1966) used the QM equations to argue that the 
applicability of such alternative formalism of RG is restricted to the standard cosmology problem — the problem of a 
homogenous, isotropic and conformally flat case. 

Although it is correct in the strict sense, this argument led to a disadvantage of the QM formalism compared to the 
other methods based on the Lifshitz program. It justifies the wide use of the complex Newmann-Penrose formalism 
(1962). Despite of it, we can prove that some objects of this formalism are physically unobservable. 

Here we will follow a simpler and more direct path, which corresponds to choosing, from the beginning, as the basis 
of our analysis, the gauge-invariant physically observable quantities. The dynamics for these fundamental quantities 
will then be analyzed and any remaining gauge-dependent objects which one can deal with will be obtained from this 
fundamental set. 
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There are basically two fundamental approaches by which the perturbation theory can be elaborated: one of them 
makes use of Einstein's standard equations (Lifshitz 1963) and the other is based on the equivalent quasi-Maxwellian 
description (cf. Jordan 1961, Hawking 1966 and Novello 1983). In this paper, we will focus on the second approach. 

A. Perturbed Quasi-Maxwellian Equations 

We state here the perturbed linearized quasi-Maxwellian equations for gravity. They shall be used in the following 
sections to treat the dynamics of the perturbed quantities, writing all these dynamical variables in the form 

-A (perturbed) -^(background) ~t~ (^^)- 

After straightforward manipulations, we obtain the perturbed QM equations as follows 

h u a h^(SE^)' + 6(5E a P) - \{5E y [a )hP\V^ + ^if v ^r ] a ^ x V tl V T {5E eX )h iy + 
-^(SH^).„hfy^ x V T = -\{p + p){6a a f 3 ) + ^(Sq^ - \h^> [Sq^), v + (101) 

h u a h/(SH^)' + 6{5H a ?) - \{5H^ a )hP\V^ + ^^r ] aX ^V n V T {6H ei )h x ,+ 

- J(<f£ A ") ; ,W )T7A v 7 = |/^V )£r '%(<5n„ e ) ;T , 



and 



The perturbed equations for the kinematical quantities are 



and 



(102) 



(6H a ^), v h a "h^ = (p + p){5uj") - X - rf^V^ {8q a ), p , (103) 



(SE^h^hT = -(Sp), a h ae - -p(8V e ) - - Pt0 (SV°)V e + -h e a (SIl a ») ;fl + ^ (<5<f ) . (104) 



{SO? + 6(SV°) + h{59) - (Sa a ), a = - (1 + 3A) (6p) , (105) 

+ \h^{5a a ), a - h$a { _ a ).ff ) hp a h v P + \o{Sa^) = -(5E UU ) - ^(SIL^) , (106) 

+ \e{8w») = I^(<5^) ;7 V4 , (107) 

^(59),xh\ - + ^(^°)^° - + 8u a p) ia tf* = -(%) , (108) 

(Suj a ), a ^0 (109) 

(SH UU ) = -\h a { y y) {{5a ai ), x + {5u) aj )., x )Vf} e7X V e . (110) 

The perturbed equations for the conservation of the energy-momentum tensor yields 

(Sp) m + p(SV°) + 6(Sp + 8p) + (p + p){59) + (Sq a ), a = , (111) 

and 

i>(5V u ) + p,o(<5V%° - (Splfi^,, + (p + P )(Sa u ) + h ua (Sq a )' + p(Sq u ) + h na (67T af3 ), p = 0. (112) 
Let us give some examples of how this method of perturbation works by considering the solutions derived last section. 
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B. Schwarzschild Solution 

The analysis of the stability of the Schwarzschild geometry in the work of Regge and Wheeler (1957) uses the 
standard Lifshitz method. These authors analyze the spectral decomposition of linear perturbations in terms of two 
fundamental modes, called "evens" and "odds" , and conclude affirmatively about the stability of this geometry. The 
main difficulties presented 3 were solved by Vishveshwara (1970) by means of a convenient coordinate transformation. 

An important characteristic of these works, which are not conveniently emphasized in the literature, consists on 
the impossibility of application of the harmonic decomposition. In effect, "even" and "odd" modes are both obtained 
from the scalar derivation, where only the degrees of freedom of spin 1 are present. Klippert (1998) has shown that it 
is possible to develop a systematic way to appropriately apply the QM formalism to the Schwarzschild solution, even 
in Stewart sense — cf. Stewart (1974). Let us review briefly this approach constructing a manifold's basis. 

1. Construction of the Schwarzschild Basis 

The complete set of eigenfunctions of the Laplace-Beltrami operator (V 2 ). The manifold is a sub-manifold orthog- 
onal to = Sq and we deal with the Schwarzschild geometry in Gaussian coordinates. The Laplacian is formally 
written as 

V 2 = V Q V Q 

where V Q = h a @V « and h a ^ = 8 a ^ — V a V^ . The scalar component of the Schwarzschild base is a function Q(x a ) 
such that 



V 2 Q = -A S Q 



(113) 



explicitly in terms of the metric (52) yields 



1 d ( r 2 dQ 

r 2 ^a 2 -A^R \V^^A m 



d 2 Q 



AQ, 



(114) 



where A(T,R) = 1 — 2M/r(T, R). The angular variables will give origin to the spherical harmonics. Afterwards, we 
consider the particular case in which a 2 = 1 to get the easier differential equation below for the basis Q in terms of 
T and R coordinates 



3fT-i- j?\d 2 F ,5dF 



j T + /e +AV2M(-j\/2M(r + i?) 



1/.", 



F = 0, 



(115) 



where it was assumed that Q(T,R,8,(f>) = F(T, R)Yj n (6, (jf) and Y^O,^) are the spherical harmonics. We can 
completely integrate the equation above to obtain the general solution in terms of Bessel functions, as follows 



F(T, R) 



{T + R) 



1/3 



[a J a (x) + bY a (x)} 



(116) 



where we defined 



a 



V1-81-8P 



_V6c(T + R) 



2/3 



being c = — Av / 2~M[(3/2)v / 2l\7] 1 / 3 and, a and b constants of integration. J a {x) and Y a (x) are Bessel functions of first 
and second kind, respectively. 



3 At those times, everybody knows whether the divergence presented by the perturbation at the event horizon was real or only an effect 
caused by choice of coordinates. 
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2. Gauge-Invariant Variables 

From the subsection [II A] , the decomposition induced by V^ 1 = Sq leads to a degenerated shear tensor (two identical 
eigenvalues) proportional to the electric part of the Weyl tensor. We thus introduce the following geometrical objects 



Y = rr - 2CT „ a „ ,2(T 2 2a 2 h 



E a ^ a 
-to*-' 



(117) 



7 _ tt _ H a pCjP a 



We remark that these tensors present some particular algebraic features: they are symmetric, traceless, orthogonal 
to shear (X a ^a l3 a = Y a pa 13 a = Z a pa^ a = 0) and the most important characteristic is that they are null on the 
background. Therefore, they constitute a set of "good" variables (in Stewart's sense, 1974) to do perturbation theory 
for the Schwarzschild case. 

3. Dynamics 

Using the QM-equations, we can calculate the propagation equations of X^ u , and Z^ along to the geodesies 
represented by the vector field V*. It is useful to rewrite the outcome in terms of these objects, to get a closed 
dynamical system. We restrict ourselves to the exhibition of the propagation equations for the perturbations associated 
to the gauge invariant variables as follows 



5X, 



( 



a 

2a 2 



2a (^5X v )x 



SY L 



h%hZ - 2^-h a {fi a u f + J^la"? - X^)(a^ - X^ + 2^^V 



5 a 



9 

2 rr(Tn 



(118) 



SY„ 



A a 3 a : E y 



a : E 
~2a 2 ~ 



0) 6Y„ 



3a\JY 



2a 



h\ am ^V T SZ Xe . n 



5H a /3 t \ g 



a : E 
~2a 2 ~ 



<m 



a/3 



o- A ( Q <5rbj) A 



a : E 



2a 2 



+ 



5q. 



a a Hn aB + 2^-^-8a\ ) h„ 



2a 2 



(119) 



and 



a : E x, , 
"2a 5- ;7 



a : E 



(120) 



6 



2a 2 



5Z„ 



&t x ,JZ, 



These equations (118) (119) and (120) are not completely independent. They must satisfy constraint equations, 
which rewritten in terms of the geometrical objects described above are 



2a 
-(h? 



h a t a^ u a a 
2Ka< 



x a f 



2a 2 2a J 



h a f 



(a 3 /2a 2 y, a 
(^/2a 2 ) 



h a f 



2a 2 2a 2 



(121) 
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[/C^^W + §fc?0. Q + (a eA + w eA )a A - g e ] + 




(122) 



and 




+ ^#ff eA w A + (p + p)w e + 



(123) 



Applying the Lifshitz method, we can decompose the gauge invariant variables set in terms of the basis Q(x a ) and 
analyze the stability of the dynamical system for the scalar perturbation, but such analysis will not be pursued in 
this work. More details can be seen in Klippert (1998). 



Kasner universes constitute a paradigm of the Bianchi-type I anisotropic space-times. We shall follow the procedure 
used last section and present a minimal closed set of gauge-independent observablcs for an adequate basis built for 
this specific background, and subsequently employed in a dynamical system that is written in the framework of the 
quasi-Maxwellian equations. It is then found that the method can be carried out to its end and a closed dynamical 
system obtained. All three types of perturbation — scalar, vectorial and tensorial — are presented and discussed in the 
same way, but we limit here our analysis to the scalar perturbation case. 

From the seminal work of Bclinsky, Khalatnikov and Lifshitz (1970) it has been shown that — for any kind of regular 
matter satisfying the usual energy conditions (cf. Hawking and Ellis, 1977) in the neighborhood of singularity — the 
Bianchi-type I Kasner solution works as an attractor for all the other solutions. In this sense these geometries are 
good paradigms for anisotropic models, which have been extensively analyzed in scientific literature (sec Novello 1990, 



The problem of stability of anisotropic cosmological models and the analysis of perturbations has also been exten- 
sively studied in the literature using the method based on the perturbations of the metric tensor and the dynamics 
determined by Einstein's equations by Mataresse (1993-5), Miedema (1994), Noh (1995), Ibahez (1995) and Mutoh 
(1997). 

In order to apply the linear perturbation theory, we will obtain basis which arc analogous to the spherical harmonics 
bases. The dynamical system is then studied in the framework of the quasi-Maxwellian (QM) equations. However, 
a slight change in the method must be effected at this point, due to the existence of non-null tensorial quantities 
such as the shear and the Weyl tensor in the background, in analogy to the Bardeen (1980) variables. The dynamical 
equations for these extra quantities are obtained from QM equations. 



In order to make the temporal dependence of the perturbations explicit in the quasi-Maxwellian (QM) equations, 
it is necessary to obtain a basis in terms of which all perturbed quantities can be written. Since we are dealing here 
with an anisotropic background, we shall not writte in terms of the spherical harmonics. A new basis must then be 
constructed. 

In this section, the three types of basis - scalar, vectorial, and tensorial - are exhibited. Since we are considering 
a background free of any matter, an apparent difficulty appears regarding the adequacy of analyzing matter-related 
perturbations. However, it is possible to write a generalized solution for the specific Kasner model, which presents 
matter-related terms that are of a lower order than the geometrical ones. Therefore, although it has been argued that 
the matter-related terms do not contribute to the unperturbed background, they might have an important contribution 
after perturbation. 

The scalar basis: In order to be able to obtain a scalar basis {Q(x,y, z)} for the Kasner background, we will 
impose that it satisfies the equation 



C. 



Kasner Solution 



1998). 



1. The Anisotropic Basis 
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V 2 Q = n 2 Q, (124) 
where n 2 is a function of time. In Cartesian coordinates, this equation is integrated to give 

Q(x, y,z) = N exp[-i(mx + n 2 y + n 3 z)} = Ne~ miX ° , (125) 
where TV and rij are arbitrary constants (a choice made to simplify our calculations), and the following relation: 

n 2 = -/i Q(9 n Q n' 9 , (126) 

is inferred. 

Using results (124) and (125), we proceed to write vector and symmetric traceless second-order tensor bases, which 
will be used to define the corresponding perturbed quantities 

Qa = V Q Q = in a Q 7 (127) 

Qa/3 = VpQa - ^n 2 h al3 Q, 

where the properties of symmetry and tracelessness of the tensor can be directly shown. From the vector definition, 
Q a has only spatial components. The tensor Q a p written in terms of the scalar Q is 



Qa/3 = - ( n a np + —h a p ) Q. (128) 



The following properties are then obtained: 

V a Q a =n 2 Q, 

V 2 Q Q = n 2 Q a , 

\7 2 Q a p = n 2 Q aP (129) 
Q = 0, 

{Q a )'= (<J a p + Qh a p)Q p , 

2 ~ A 2 2 2 

(Qafs)' = — V(a^Ql3)m ~ Tj" & afiQ + Tthapa^ 1 'Q^v 

Let us choose a specific direction of propagation for the scalar basis. 4 Therefore, we set 

ni = n 2 = => n 2 = t- 2p3 (n 3 ) 2 , (130) 
Thus, the scalar basis and its correlated quantities take the very simple form: 

Q = Ne- m3Z , 

Qa = -m 3 (0, 0,Q), 



. t 2pi 

Q af 3 = \Q o o 

-2t 2 ^ 



(131) 



The vectorial basis: Similar to the scalar case, let us impose for the vector basis {P a } the equation 



4 This procedure was adopted in several instances in the literature. See, for example, Sagnotti (1981). 
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V 2 P Q = m 2 P Q , (132) 
where m 2 is a function of t. Integrating this equation, one finds 

P a =V a e- im ^ i '. (133) 
We choose P° Q and rrij as constants. 5 From Eqs. (132) and (133), it follows: 

to 2 = -h af) m a mp. (134) 

In order that {P a } constitute a basis, two properties must be valid. The first one concerns the fact that P Q must 
be spatial quantities. This is immediately satisfied, upon the choice V°o = 0, or 

V a P a = 0. (135) 

This property must also be preserved in time, i.e., 

(V a P a )=0. (136) 

which is identically valid for P^ =const. 

The second property, namely that no scalar quantities can be obtained from the vector P Q , implies that 

V Q P Q = 0, (137) 

and can be written as 

h ap m a P fj = 0. (138) 
This property must also be conserved in time; hence it follows 

V"P„ = =^> (T al3 m a Pp = 0. (139) 
The conditions (138) and (139) can also be written in terms of the V° a - 

t- 2pi mi V^ + t- 2p2 m 2 V^+t- 2p3 m 3 V° = (140) 

and 

Pl t- 2pl miV° + P2 t~ 2p2 m 2 Vl +p 3 t- 2 P3m 3 V° z = 0. (141) 

From V a , it is possible to construct three quantities: a symmetric, traceless second-order tensor (which we will 
denote by P a p), a pseudo- vector denoted by P*, and finally the corresponding pseudotensor P*^. The respective 
definitions for these quantities are the following: 



This choice is necessary to avoid that complex terms or explicit dependences on spatial coordinates occur when the simple derivatives 
of the vector basis are calculated. 
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PaP = V( Q Pg), 

P* = rjJ^V^V^), (142) 
The first of the definitions (142) is immediately rewritten as 

Pq,3 = -im(aP^) (143) 
and its corresponding Laplacian and time-projected derivative are proven to be given by 

V 2 P Q(3 = m 2 P afi , (144) 

(AvO'=-^^-^%^W ( 145 ) 
The pseudovector P*, from the second definition in Eq. (142), is 

P* a = -iruf^VfimvPp (146) 

It follows — from the fact that all these quantities describe the perturbations — that the same properties that define 
P a should also be valid for P*. Therefore, the pseudovector should be both a spatial and a divergenccless quantity, 
which is satisfied. These conditions must be preserved in time. This is identically valid for the first property. The 
preservation condition for the null divergence property is given by 

(V a P*)' = a ap m a P; = 0, (147) 
which is then rewritten in terms of the P„ as 

(P2 - P3)m 2 m 3 P° + (P3 - Pi)mim 3 Vy + (pi - P2)m 1 m 2 'P® = (148) 
In addition to the above condition, the following useful results are also obtained: 

V 2 P* = m 2 P*, (149) 

and 



(P*)-=--9P*-aP a Pl (150) 



The last quantity to be considered is the symmetric, traceless pseudotensor P*a- From the third definition (142) 
we get 

P*p = -im^P^, (151) 

and the relations below follow immediately: 

V a P*e = m 2 P* , (152) 



V 2 P^ = m 2 P* p , 



(153) 
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and 

(p^y = -\ Qp *fs - °\<> p hr ( i54 ) 

The most general form for the vectorial basis implies obtaining suitable V® and rrij , and replacing them in the basis 
expression (133). These quantities are defined from the conditions (140), (141), and (148) as 6 



mi = 7712 = 0, 

r° = o, 



o _ n (155) 



and, using Eq. (155) above in Eq. (133), we get 



P a = er"^ z (V o x ,V°,0), (156) 

where 7773, V® and V® are arbitrary constants. In addition, the quantities P a /3, -P^, and P*a are written, in this case, 
as 



v° 





P a p = -™ 3 e" m32 I 0^ Vl , (157) 



P* = nf 123 m 3 e- mi3Z (-t 2pi V y ,t 2 P 2 V°,0), (158) 

-t 2pi P° 

P* p = v 0123 (m 3 ) 2 e-" n3Z ( t 2 ^V Q x ). (159) 

„ t 2 P i V o t 2 P2V a q 

The tensorial basis: We will begin by defining the tensor U^ u (t, x, y, z), which is written in matrix form as 

/ \ 

\0 x C 7/ 

where (a, ip, 0, 77, /?, e, \, £, 7) are functions of all four coordinates. The choice we made for our tensorial basis enables 
us to simplify somewhat the future calculations. It is easy to write the totally contravariant tensor as: 

U» v = W a U av (161) 

with U au = U ua .^ 

In order that C/ M „ be, in fact, a basis, it has to be a solution of the equation 

\7 2 U' 1 V = fc 2 l>„, (162) 
where k 2 is a function of time. Solving this equation, we obtain an explicit form for the tensorial basis 

U^ v =W v e- lk ^\ (163) 



The condition (155) above implies making a choice of a specific direction for the basis and has been made a number of times in literature. 
See, for example, Ref. [140]. 
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where the U 11 u will be taken as covariantly constant tensors and kj are arbitrary constants which are related to the 
wave number k 2 and the components of the metric tensor through the following relation: 7 



k 2 = -h jl k 3 k t = 



fci 
a(t) 



b(t) 



k^ 
c(t) 



As in the previous cases, the tensor basis must obey the following properties: 
(J) The tensor basis should be orthogonal to V a : 



VM^ = 0. 



(7T) Scalars cannot be obtained from the tensor basis: 



= U% = 0, 



or, using Eq. (160), 



(164) 



(165) 



(166) 



a + /3 + 7 = 0. 
(///) Vectors cannot be obtained from the tensor basis: 



(167) 



Vl7u„ = 0, 



which gives 



(168) 



= o, 

Equation (169) above can be rewritten, using Eq. (163), as 



(169) 



k%a + k 2 i] + k 3 x = 0, 

hi> + k 2 p + he = o, 

k\4> + k 2 e + fc 3 7 = 0, 



(170) 



It is easily seen that all the properties above are preserved in time. 

At this point it becomes necessary to define a quantity that enables us to write pseudo-tensorial perturbations. 
Therefore, we define the dual U*^ u as 



(171) 



which can be then rewritten as 



V*\ = --k.i^VxU^ + ^ x \VxU^). 



(172) 



It follows that all the properties of the tensorial basis U^ v are equally valid for the dual tensor [/*% and that they 
are preserved in time. 

We can obtain an explicit form for L/ M „. From Eqs. (60) and (160) we have that 



7 This choice was made in order to avoid terms with spatial coordinate dependence when calculating the derivative (U^ii) 
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u 2 i =v = gng 22 u} = t 2 ^-p^?p, 

U\ = X = 9ii9 33 Ui = t 2 ^-^U, (173) 
U 3 2 = C = 92 2 g 33 Ul = t 2 ^-^h. 
Using the above results and the null trace condition (167) in condition (170), we find that 

feia + k 2 t 2 (p i -p^4> + jt 3 f 2 (pi-P3) = o, 

k 1 ^ + k 2 (3 + k 3 t 2( -P 2 -P 3 h = 0, (174) 
k 1( j) + k 2 e - k 3 {a + 0) = 0. 

However, since both the kj and the W '„ arc constant, we can see that each term in the three above relations must 
also be constants. We then choose a specific direction for the basis (Sagnotti, 1981), taking 



k x = k 2 = 0, 

h ¥= o. 



(175) 



Then, equations (174) simplify to 



i2(pi-p 3 )0 = o, 

t 2{P2- P3 ) e = 0j (176) 
a + /3 = 0, 



and give, as a consequence, 



<P = e = 0, 

j8 = -a 7 = 0. 

The tensorial basis and its dual arc, then, written in matrix form as 



(177) 




U» v = W v e- l ^ z = \ r] -a , (178) 



(t^i/j + t^r)) -2i 2 P 2 a 
= -^rf 123 k 3 I -2&* a -(t 2 P^ + t 2 P^) ) . (179) 
£2p lx -i 2 ^ 



2. The Gauge-Invariant Variables and Their Dynamics 

The dynamical system for an anisotropic background will be obtained in the framework of the quasi-Maxwellian 
(QM) equations. In the Kasner background, the QM equations are reduced to the set 



(aijJ+Eij + ^6<j r] - -(2o*)hij + cx ik a k j = 0, (180) 



{EijJ + WE l3 + ^ (i E j)fi - htjO^E^ = 0, (181) 



-Q 2 + 2(t 2 = 0, (182) 
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The natural step is, then, to write the perturbed QM equations, by making the usual choice: A^ pert ^ = An,^^ + (5 A). 
However, a modification in this method becomes necessary at this point: the three non-null quantities in the Kasner 
background {a^v, and 9) should be replaced by "artificial" quantities in order to eliminate all gauge-dependent 
terms from the dynamical system equations. These new variables are written in terms of the original, gauge-dependent 
variables. Nevertheless they constitute "good" quantities, in the sense that they are zero in the Kasner background, 
as per the Stewart lemma (1974). This procedure is analogous to the one implemented by Bardeen (1980), but the 
variables obtained in the present case are much simpler, as we shall see in the next section. 

Minimal closed set of variables for the Kasner background. The starting point to obtain the new variables is the 
set of QM equations for the Kasner background, Eqs. (180) - (182). If we employ the following relations (easily 
demonstrated and specifically valid for the Kasner background): 

{a a p)=-Ba a p, 

{E a p)=-26E a p , (183) 
6 = -6 2 , 

we are able to define the new variables that are to replace the original ones as 

X a p =E a p-^e<T a p-^-h a p+a a ^ai i p, (184) 

Y a p = 0E a p + ^a a ^E% + \^pE a ^ - h a a^E^, (185) 

W = 2a 2 - \o 2 . (186) 

These three variables are easily proven to be zero for the Kasner background and, therefore, "good" ones to be 
perturbed. They may, then, replace the shear, electric part of the Weyl tensor and expansion as the new variables in 
the dynamic system. An additional simplifying choice will be made, namely the relation between the energy density 
p and the pressure p is 

p = Xp, A = const., (187) 

even after the background is perturbed. This choice — also made for the FLRW case (see Novello, 1995a, b and 1996) — 
will be considered as valid throughout this analysis. The complete minimal closed set of variables to appear in the 
dynamical system is therefore 8 



M = {X a p,Y a p,H a p,Tr a p,q a ,a a ,u} a ,W,p}. (188) 

The next step in this analysis is to obtain the complete dynamics for the new variables, X a p, Y a p, W. The resulting 
set of equations and the remaining QM equations must, then, be rewritten in terms of the new gauge-independent 
variables in M. This constitutes the dynamical system of equations which is used to study the perturbations of Kasner 
model. Such complete dynamical system as well as the steps for its derivation can be seen in Novello et al. (2000). 
The next sections will deal with the three perturbation cases and the results obtained for each one. 

3. Scalar Perturbations 

In this case, the minimal closed set of observables M. involves practically all the original variables of the system, 
which we proceed to present here in terms of the scalar basis Q: 



8 Since the acceleration - the variable a a - does not have a dynamical equation of its own, it must be eliminated from the dynamical 
system in order to be really closed. This will be achieved by fixing a value for the function (Sa a ). 
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{SX af} ) = X(t)Q a p, [5Y afi ) = Y(t)Q a p, 

(5n a p) = n(t)Q afj , (Sq a ) = q(t)Q a , 

(Sa a ) = 1>{t)Q a , (5 P ) = R(t)Q, 

(SW) = W(t)Q, (8 P ) = p(t)Q, 



(189) 



where the spatial part of the velocity, {5Vk), is also zero in the background (but not an adequate variable, since its 
value in the background depends on the choice of an observer), and it is written as 



(SV k ) = V(t)Q k . 



(190) 



Since we are dealing with scalar perturbations, the vorticity and related perturbations, for instance the magnetic 
part of Weyl tensor, are not defined — see details in Novello (1995a) and Goode (1989). The relation between the 
shear and the anisotropic pressure is still valid, but in this case the viscosity £ is also a "good" variable (i.e., it is a 
gauge-independent variable for it is zero in the Kasner background), written as 

{so = m q 

and must be considered. From Eqs. (189) and (128) it is possible to obtain (<5£) in terms of (5ir a /3) as 



i 



(2a 2 ) 



(191) 



Once again - as in the previous cases - we take the perturbed anisotropic pressure (Sir a p) as zero, in order to simplify 
the dynamical system to be solved. Thus, 

7r(t) = £(t) = 0. (192) 
Using the results obtained for the scalar basis, we then obtain the dynamical system for the scalar perturbations case 



X + BX-Y + \q + \e^) (-^n 2 X + ^n 2 i' + ^{W + R)) [a^Q]+ 



0, 



(193) 



Y+^6Y+\l 



2„2l 



n 2 Y + 6{1 + X)R [a^Q] + (n 2 q - 2n 2 ^ + W - R)[E^Q] + 



0. 



(194) 



W + -eW - 2{a afs n a n )X + 2 



{a^n a n ) + -6n 2 



ip - -0(1 + 3X)R = 0, 



(195) 



[pi(l -pi) -p 2 (l -P2)} + - 3p 3 )(pi -P2) \ ip = 0, 



(196) 



R + 9(1 + X)R + n 2 q = 0, 



(197) 



q + 6q - Xtr 8/3 R = 0, 



(198) 
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ln 2 X-9q+ l -{W-R) 



Qa + 



0. 



(199) 



Equation (196) is satisfied in two cases: (1) p\ ^ P2 => ip = and (2) pi = p2 => isotropy plane, and the simplest 
choice here is (1), with -0 = 0. 
The system becomes thus 



X + 9X-Y + -q)Q^ 



-ln 2 X + ^{W + R) 



(200) 



Y + pY + \i 



-\n 2 Y + 6(1 + X)R [a^Q] + (n 2 q + W - R)[E^Q} + 



0. 



(201) 



W + \BW - 2{(j a ^n a np)X - |<9(1 + 3X)R = 0, 



(202) 



R + d(l + \)R + n 2 q = 0, 



(203) 



At~ 8 / 3 i? = 0, 



(204) 



ln 2 X-6q+ l -{W -R) 



Qa = 0. 



(205) 



Although Eqs. (200) and (201) cannot be factored out in the scalar basis, the remaining equations in the system, 
Eqs. (202), (205), can be separately integrated. The constraint (205) eliminates the variable X from the reduced 
system 



3 9 



1 1 



2n 2 2n 2 2n 2 



(206) 



so that the dynamics of W, Eq. (202), is written as 



W 



w 



R - 3—t 8/3 (a aP n a n p )q = 0, (207) 



and, with the relation 



;(<T a Pn a np) = -6{l-p 3 ), 



(208) 



the final dynamics for W thus becomes: 



W + (1 - P3 )6 W - (1 + 2A - P3 )6R + {3 P3 - l)6t 8/3 q = 0, 



(209) 



and the reduced dynamical system — which is closed in the variables W, R, and q — is given by Eqs. (203), (204), and 
(209). As in the previous cases, this reduced system can be solved for the following Ansatz: 



q(t) = qvt x , 
R{t) = R t v , 
W(t) = W t l 
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(210) 



where qo, Rq, Wq and x, y, w are constants to be determined. An immediate inspection of the powers of t in the three 
equations gives the relation between x, y, and w 



y = w = x + -, 



(211) 



while the rest of the equations gives the following results: 



Ro = -r(x + l)<Zo, 



and 



W = [(x + 1)(1 + 2A - pa) + A(l - 3p 3 )](w + 1 - p 3 rV 
where w =/= (p 3 — 1) must be valid, and 



(212) 



(213) 



(n 3 ) 2 = + [ 



(214) 



The results above are not applicable to the case A = 0; the details of this specific solution can be seen in [121]. 

In order to analyze the stability of the above solution, we will impose that (n 3 ) 2 be positive. From Eq. (214), we 
obtain - considering A as positive 9 - two possibilities: (1) x > —1; (2) x < —(8/3 — A), and we have the following 
results for each case, exhibited in Table II. 



Value for x 


?(*) 


R(t) and W(t) 


x < -(8/3 + A) 


stable 


stable 


-1< x < 


stable 


unstable (faster than t 2 ^ 3 ) 


x = 


constant 


unstable (faster than t 5 ^ 3 ) 


x > 


unstable 


unstable (faster than t 5 ' 3 ) 



TABLE II: Stability analysis results for scalar perturbations. 



We see that, although unstable solutions exist, the instability is — as we will see in the case of tensorial 
perturbations — not catastrophic. Another point which should be mentioned is that the matter density is unstable for 
any choice of the constant exponent x that implies (tt-3) 2 positive. 

This concludes the analysis for the scalar perturbations case (A ^ 0), and the case with zero A must now be studied. 
For A = 0, the reduced dynamical system is given by 



The case of A < gives (ng) 2 < for all values of x, and therefore will not be considered here. 
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W + (1 - P3 )9W - (1 - P3 )0R + (3p 3 - l)9 2 t 8 / 3 q = 0, (215) 
R + 0R + n 2 q = O, (216) 
g + 6>r/ = 0. (217) 



Equation (217) can be directly integrated as 



where go is a constant, while the ansatz 



q(t)=q t- 1 , (218) 



R(t) = R t y , 
W(t) = W t l 



(219) 



can be used in Eqs. (215) and (216) to give 



y = w = |, 



(w + 1 -p 3 )' 

The above results give a partially unstable solution for A = 0, with R and divergent in t, while q tends to zero for 
t — > CO. 

4- Vectorial Perturbations 

The perturbations associated to the state of motion of a fluid (without taking energy density perturbations into 
account) are described, in principle, by the following minimal closed set of variables: 10 

M = {X al3l Y a/3 ,H al3 ,q ai a ai uj a }. (220) 
The perturbed quantities can then be written in terms of the vectorial basis P a as 

(SX a p) = X(t)P aP , {SY aP ) = Y{t)P aP , 
{8H a0 ) = H(t)P*p, (5q a ) = q(t)P a , 

(5a a ) = tp(t)P a , (Suj a ) = n(t)P*, ( ] 

(5V k ) = V(t)P k , 

where the variable V(t) is, again, not an adequate one, since it depends on the initial choice of an observer. In 
order that the basis can be successfully factored out from the dynamical system equations, we are forced to choose 



The same observation regarding the variable a a that was stated for the tensorial case holds here as well. 



35 



between two possibilities: (1) Eliminate one of the components of the basis, say V®, and (2) analyze solely the case 
of a background with an isotropy plane (namely, the Kasner solution) . 

We will consider here the choice (1) above, since it gives a more general result (and also contains the specific case 
(2), as it shall be shortly seen). The dynamical system for the vectorial case can then be obtained by factoring out 
the basis, as 

ui 2 

X + 6X-Y+ —H + -q- -9 P2 yj = 0, (222) 
Y - \e(hp 2 + l)Y + ^d 2 p 2 (p 2 - 1)X + \em 2 {p 2 - 1)H - ^0(p 2 - l)q - \e 2 p 2 (p 2 - lty = 0, (223) 
H - h(5 Pl + 1/3)H -\x + ifl(3(p3 - P2) + 2/3)0 + \^{P2 - P 3 )(3 Pl + 1/3)^ = 0, (224) 



g + 6> 9 = 0, (225) 



H + \o(6p2 + 1)0 - \t- sfs ijj = 0, (226) 
m 2 t -s/3 X + Q m 2( p2 _ !/ 3 )jj + I(9m 2 (pi + 1)0 -6q+ \d 2 { Pl - V 3 )(Pi + W = ( 227 ) 

m 2 ff - 0(p! - p 3 )X + 39 2 t^ 3 p 2 {p 2 - 1)0 + ^9 = 0. (228) 

The system (222) - (228) is not closed because there is no dynamics for the acceleration tp(i). However, the 
constraint (227) can be employed to eliminate this variable in terms of the other variables in M, closing therefore the 
remaining dynamical system for a given background, specified by the triad (pi,P2,Pa). In order to obtain a specific 
solution for the dynamical system, we must make a choice for (pi,p 2 ,P3)- Here we solve the specific case of the Kasner 
solution with an isotropy plane, which is algebraically simpler to solve. This choice, applied to Eqs. (222) - (228), 
gives the following dynamical system: 

m 2 1 1 

X + 6X-Y + — H+-q+-9i> = 0, (229) 

Y + ^-9Y - \e 2 X + -Om 2 H + \oq + \pij> = 0, (230) 
o o 4 8 



H+^-BH-\x- ±-0n + 1-^ = 0, (231) 
o 2 12 12 m 



q + eq = 1 (232) 



+ |#Q - ^t-^ 3 tp = 0, (233) 



3G 



15 5 
m 2 t~ s/3 X + -9m 2 H + -9m 2 VL - 9q + —9 2 ijj = 0, (234) 
3 6 18 

m 2 H - ex - ^e 2 t s/3 n + -q = o. (235) 

The acceleration tp is given in terms of the other variables by (234), which then closes the system 

18 18 m 2 o,o 6 m 2 m 2 

t-^-TlF* ' X -IT H -'T^ (236) 

and the closed dynamical system is then written as 

X + ^(59 2 - 6tr 8/3 m 2 )X -Y+ \m 2 H - m 2 fl + ^-q = 0, (237) 

i o i ii OQ 

Y + —BY (56> 2 + 9t- s/3 m 2 )X + —9m 2 H - -9m 2 n + —6q = 0, (238) 

6 15 20 2 40 

H + He H --t- 8 ^x--en + -^ q ^o 1 (239) 

15 10 6 10 m 2 ^ ' V ' 

O + -I(i # 2 - 9t- s / 3 m 2 )n - -^t- 16 / 3 X - -—t- s / 3 H - —t-*' 3 q = 0, (240) 
69 5 5 9 59 

q + 9q = 0, (241) 



9X - m 2 H + -9 2 t s/3 fl --q = 0, 
3 2 y 



Equation (241) can be directly integrated in q(t), giving 



X(t) = X t x , 

Y(t)=Y tv, 

H(t)=H t z , 

n(t) = n t w , 



(242) 



q{t)=q r 1 , (243) 
where go is an integration constant. The rest of the system can be solved for the following ansatz: 



(244) 



where Xq, Yo, Ho, f^o, x, y, z, and w are constants to be determined from Eqs. (237) - (242), as well as the constant 
7713. The exponents are easily obtained as 

x = 0, 

V = -1, (245) 
z = w = —5/3, 

and the remaining constants must satisfy the following conditions: 
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(246) 



8[5 - 3(m 3 ) 2 ]X - 20Y + 2(m 3 ) 2 ff - 20(m 3 ) 2 ft + 29g = 0, 

8[5 + 9(m 3 ) 2 ]X - 380Y - 6(m 3 ) 2 i?o + 60(m 3 ) 2 O - 87g = 0, 

78(m 3 ) 2 X - 68(m 3 ) 2 ff + 55(m 3 ) 2 ft - 63g = 0, 

54(m 3 ) 2 X + 18(?n 3 ) 2 H - 5[11 - 9(m 3 ) 2 ]fi + 54g n = 0, 

6X - 6{m 3 ) 2 H + 4fi - 3g = 0. 

A simple but rather tedious manipulation gives four of the constants (say Xq, Yq, Hq, and f2o) as proportional to the 
fifth (qo), as well as polynomials of (m 3 ) 2 : 



where 



X 
Y -- 



Pi 
Pi 



3 



1 

TO 



§ [5 + 9(m 3 ) 2 ]X - ^(m 3 ) 2 ff + 3(m 3 ) 2 O - |g 9o 



Mi[l + 3(m 3 ) 2 ]- 1 
[450(to 3 ) 4 + 1434(m 3 ) 2 - 935] ' 

M 2 [l + 3(m 3 ) 2 ]- 1 
[450(m 3 ) 4 + 1434(m 3 ) 2 - 935] ' 

[39(m 3 ) 6 + 5(m 3 ) 4 - 5(m 3 ) 2 + 8] 



Pa 
Mi 

M 2 = 900(m 3 ) 8 + 2523(m 3 ) 6 + 3591(m 3 ) 



[450(m 3 ) 4 
4455(m 3 ) 8 - 



- 1434(m 3 ) 2 ~ 
5328(to 3 ) 6 - 



935] 

1524(to 3 ) 4 



- 3583(m 3 ) 2 + 598, 
3465(m 3 ) 2 - 240. 



The condition over (m 3 ) 2 is then obtained as 



(247) 



24[15 - ll(m 3 ) 2 ]X + 22(m 3 ) 2 ff - 220(m 3 ) 2 O + 319g = 0, 

which — upon substitution of X , H , fl and go — gives a fifth-order equation on (to 3 ) 2 that is satisfied for at least one 
positive value of (m 3 ) 2 namely, (m 3 ) 2 « 2.9252, and the consistency of the ansatz is then proven. The solution for 
the vectorial perturbation case, for a Kasner background with an isotropy plane is then written as 



q{t) = q Q t \ 
X(t) = X , 
Y(t) = Y t-\ 
H(t) = H t- 5 / 3 , 
fi(i) = rW~ 5/3 , 



(248) 



with constants Xq,Yq, Ho, and Qo given in terms of go by means of Eq. (247). Differently from the tensorial case, as 
we shall see, it is seen from Eq. (248) above that the only existing solution is stable, although the decrease is not fast. 



Tensorial Perturbations 



We will proceed in this section to analyze the case of gravitational waves in the Kasner background. Here, the 
minimal closed set A4 will reduce to the four tensorial variables: 



M = {X a p,Y a p, H a p,Tr a p}. 



(249) 
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However, from causal thermodynamics, a relation between the shear and the anisotropic pressure can be obtained as 



T(lTij )' + 7Tij = ^(Xij , (250) 

where r is the relaxation parameter and £ the viscosity parameter. We will make here the same choice made for the 
FLRW case and consider r as negligible; the viscosity will also be taken roughly as a constant, 11 which makes Eq. 
(250) 



T« = frij- (251) 

This result poses a problem after perturbation, for it would then imply that the anisotropic pressure — a "good" variable 
in the sense of Stewart (cf. Novello, 1983), one that is zero in the background and therefore gauge- independent, can 
be written in terms of the shear — which, in our case, is nonzero and (in the sense of Stewart) coordinate dependent; 
therefore, a "bad" variable. This apparent dilemma can be solved by taking the viscosity as a "good" variable itself 
(being zero in the background and therefore gauge independent as well). However, this is a scalar quantity and, as 
such, not defined for tensorial perturbations. The solution to this problem is, then, to consider the viscosity itself as 
zero, thus writing 

(tfTTy) = £(*<7«) = 0, (252) 

after perturbation, so that the consistency of the dynamical system is maintained. This further reduces the set M.: 

M = {X afi ,Y aP ,H als }. (253) 
We can, at this point, expand the perturbed quantities in M in terms of the tensorial basis U^p as 

0^) = £(„)*(*) (w) cWj, 

(«l ri i)=E(n) y (*) (n) %) i i> ( 254 ) 
(<^i) = E ( „)#(*) (n) #%)V 

We will deal here with linear equations so that we are able to suppress the summation and the extra indices from 
now on, since we can deal with each component (n) separately. The same reasoning will be applied to the vectorial 
and the scalar cases as well. 

With Eq. (254) above, a perturbed dynamical system can then be written. Starting from the original QM equations, 
and rewriting them in terms of the new variables above, we can exhibit the perturbed dynamical system for the 
tensorial case as follows: 



' a pV 1 hb{8H%).p + ^JV^ X (6H^) ;0 = 0, 

h£hZ(5Y a i}J+ 26{5Y» V ) - \<r£{5Y a v ) - ^(SY^ a ) + h^{5Y^ a ) + §££(<5X%)- 
+$E°}(&X* a ) - KE«{5XP a ) + ^^V^K^H^.p + ^6ri„'y a %h1(6H*).p+ 



(255) 



(256) 



11 In the noncquilibrium thermodynamics, both r and £ are functions of the variables of equilibrium of the system, such as the density p 
and the temperature T. 
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2o$(8Xi) = Q, (257) 

(258) 

h<$h^(5X^) ;v + rf^V p a^(SH2) = 0, (259) 

h%h^(5H^). v + n a ^ v V p cx^(8X2) = 0. (260) 

This dynamic system has to be decomposed in terms of the tensor basis and then solved for the specific Kasner 
solution. However, it is immediately seen (upon writing the perturbed terms in the tensorial basis) that the following 
restriction on the background must be accepted in order that the basis can be factored out from the equations 

Px=V2, (261) 

which implies the existence of an isotropy plane in the Kasner original, non-perturbed background. There are two 
such solutions, named Kasner and Milne, and defined by 

(1) Kasner solution 

pi =p 2 = 2/3, 
P3 = -1/3, 

(2) Milne solution 

Pl=P2= 0, 
P3 = I- 

Here we will only consider the Kasner solution, since the Milne case has already been analyzed by Novello (1995b). 
Some additional choices will be made over the tensorial basis which, while not constituting a material change on 
the /7 M „, will simplify the algebraic steps to be taken towards a final closed dynamical system 

3 



whereupon we write 



u\ = u A 2 = o 



a 

U^ u = | -a | , (262) 





V*^ v = --n 0i2i k 3 -2t 2pi a . (263) 



Proceeding to analyze the dynamical system, we find that Eq. (257) is identically satisfied, since 

(o- a pU a ) = 0, (264) 
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while Eqs. (259) and (260) arc also identically valid through similar arguments. The remaining equations can then 
be rewritten, giving 

X + %QX -Y + k 2 H = 0, 

Y + 6Y + %B 2 X + 26k 2 H = 0, (265) 
H + ^9H - t-^I^X = 0, 

which constitutes a closed dynamical system in the variables (X, Y, H) . This result is analogous to the ones obtained 
in the FLRW case, but in this case the system can be completely solved by using the relations 

9 = t~ x , 

k 2 = t- 2 / 3 (k 3 ) 2 , k 3 = const, 
and considering a simple form for the desired solution, in terms of powers of t, 

X(t) = X t x , 

Y(t) = Y ty, (266) 
H(t) = H t w , 

with Xq, Yq, Hq, x, y and w as constants to be determined. 



Value for w 




X(t) 


Y(t) 


H(t) 


w < -2/3 




stable 


stable 


stable 


io = -2/3 




null 


null 


stable 


-2/3 < w < -1/3 




stable 


stable 


stable 


to = -1/3 




constant 


stable 


stable 


-1/3 < w < 




unstable 


stable 


stable 


10 = 




unstable 


stable 


constant 


< w < 2/3 




unstable 


stable 


unstable 


to = 2/3 




unstable 


constant 


unstable 


w > 2/3 




unstable 


unstable 


unstable 



TABLE III: Stability analysis results for tensorial perturbations. 



Replacing this ansatz in the dynamical system above, Eq. (267), we obtain the following equations: 

(3a; + 7)X t( x ~ 1 ' ) - 3Y t v + 3(/c 3 ) 2 #ot ( "'~ 2/3) = 0, 

2X^ X -^ - 3(y + l)Y ty + 6(fc 3 ) 2 i7ot ( "'~ 2/3) = 0, (267) 
3X t(- x -^ - (3w + 2)H tt w - 2 / 3) = 0. 
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It is immediately seen that the only non-trivial solutions satisfy the conditions 

2 1 

y = x — 1 = w — — w = x — — , (268) 
and the dynamical system reduces to the following conditions on the triad (Xq,Yq,Hq) and the constant (fc 3 ) 2 : 

(3a: + 7)X - 3Y + 3(fc 3 ) 2 # = 0, 

2X - 3(y + 1)Y + 6(k 3 ) 2 H = 0, (269) 
3X - (3w + 2)H = 0. 

These conditions then give Xq and Yq in terms of Hq and the exponents x, y, w: 

_ (3w + 2) 
— 3 °' 

(270) 

_ 2 (^ + 2)(3 W + 2) 
y o - ^ (^+3) ^o, 

which, upon employing Eqs. (268) and (270) in the ansatz (266), gives 

X (t) = &£+£tV*H(t), 

Y(t) 2 (3^ + 2) 2 2/3 (271) 
H{t) = H Q t w . 

We also determine the constant (fc 3 ) 2 in terms of x, y, w as 

(h) 2 = ~(3w + 2) 2 . (272) 

If we take the arbitrary constant Ho as positive, an analysis of the stability of the solutions above, Eq. (271), gives 
the results exhibited in Table III. 

From this table, it is immediately seen that tensorial perturbations of a Kasner background may present — upon a 
choice of the exponent w — the same kind of instability which will be present in Fricdman-Lcmatre- Robertson- Walker 
spacetime. This instability is not of a catastrophic kind (as in the Einstein model), but is rather gradual. This is a 
reasonable development, since we are not interested in eliminating the background anisotropy upon perturbation. 

D. Friedman Universe: Scalar Perturbations 

In the case of the spatially homogeneous and isotropic FLRW cosmological model, the vanishing of Weyl conformal 
tensor suggests that the QM approach is more useful. Therefore, the variation of Weyl conformal tensor 8W a pu V is 
the basic quantity to be considered, once there is certainly no doubt that SWap^ is a true perturbation, which can 
not be achieved by a coordinate transformation. This solves ab initio the gauge problem that was pointed out before. 

From a technical standpoint, instead of considering tensorial quantities, one should restrain oneself to scalar ones. 
There are two ways to implement this: 

• Expand the relevant quantities in terms of a complete basis of functions (e.g. the spherical harmonics basis); 

• Analyze the invariant geometric quantities one can construct from and its derivatives in the Riemannian 
background structure, that is, examine the 14 Debever invariants. 

Either way, we shall see that the net result is that there is a set of perturbed quantities which can be divided into 
"good" quantities (i.e., the ones whose unperturbed counterparts have zero value in the background and, consequently, 
Stewart's lemma (cf. Stewart 1974) guarantees that the associated perturbed quantity is really a gauge-independent 
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one) and "bad" ones (whose corresponding values in the background are nonzero). One should limit therefore the 
analysis only to the "good" ones. 

This same kind of behavior occurs for the geometrical structure of the model for both the kinematic and dynamic 
quantities of matter. Therefore the "good" quantities which constitute the set of variables with which we work should 
then be chosen from these particular scalars that come from these three structures: geometric, kinematic and dynamic. 
Does that mean that the present approach effectively avoids the gauge problem? 

To answer this question affirmatively one should be able to exhibit a set of "good" variables in such a way that 
its corresponding dynamics is closed. That is, if we call Ai\A\ the set of these variables, Einstein's equations should 
provide the dynamics of each element oiM[A], depending only on the background evolution quantities (and, eventually, 
on other elements of A4[a])- This would exhaust the perturbation problem and we shall show that this is indeed the 
case. 

What should be learned from this discussion is that one should then understand the gauge problem not as a basic 
difficulty on the perturbation theory but just as a simple matter of asking a bad question 12 . One could imagine — 
which has been used a number of times in the literature (Hawking 1966, Olson 1976 and Mukhanov 1992) — that for 
FRW cosmology the perturbations of its main characteristics (the energy density, Sp the scalar of curvature SR and 
the Hubble expansion factor SO) would be natural quantities to be considered as basic for the perturbation scheme. 
However, these arc not "good" scalars, since they are not zero in the background 13 . Wc shall see in the next sections 
which scalars replace these ones. 

Consider FRW geometry written in the standard Gaussian coordinate system as in Eq. (64). The 3-dimensional 
geometry has constant curvature and thus the corresponding Ricmannian tensor Rijki can be written as 



Rijkl = —ejijkl- 

where Jijki = liklji — lilljk was defined in Eq. (72). For a while, it is necessary to distinguish covariant derivative 
in the 4-dimcnsional space-time by the symbol (;) and the 3-dimensional derivative by (||). 

Since the original Lifshitz paper, it has shown to be useful to develop all perturbed quantities in the spherical 
harmonics basis. Once we are limiting ourselves to irrotational perturbations, it suffices to our purposes to take into 
account only the scalar Q(x k ) (with Q = 0) and its derived vector and tensor quantities. We have thus 

% = (273) 
where the scalar Q obeys the eigenvalue equation defined in the 3-dimcnsional background space by: 

V 2 Q = mQ, (274) 
where m is the wave number of the scalar eigenfunction, with 

{q 2 + 1, < q < oo, e = 1 (open), 
q, < q < oo, e = (plane) , (275) 

n 2 — 1, n = 1, 2, . . . , e = — 1 (closed), 

and 

V 2 Q = f k Q Ak = y ik Q, i;k , (276) 
where the symbol V 2 denotes the 3-dimcnsional Laplacc-Bcltrami operator. The traceless operator Qij is defined as 

Qij = —Qij - l^Qlij, (277) 



Let us point out that some of the gauge dependent terms are particularly relevant, 8p including. 
However, as we shall see soon, we can construct associated "good" quantities in terms of these scalars. 
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and the divergence of Qij is given by 

Qfk = 2 Q + £) Q l - ( 278 ) 

We remark that Q is a 3-dimensional object; therefore indices are raised with 7 1 - 7 , the 3-space metric. 

In Debever (1964), the complete 14 algebraically independent invariants constructed with the curvature tensor 
were presented. Considering that we are using an adimensional metric tensor, we can classify them with respect to 
dimensionality as follows: 



Dimensionality 


Invariants 


L- 2 


h 


L- 4 


h> h, h 


L- 6 


h, h, h- hi 


L- s 


h, ho, I13 


L- 10 


In, In 



The expressions for these invariants are: 



h 


= w a ^ u w a ^ : 


h-- 




h 




h -- 


- c n^v 


h 


= W a ^»*W a p flu , 


ho 


= D^D^, 


h 


= w a ^ a w pa ^*w fwaP , 


hi 


= C aP D^D^ 


h 




I12 


= D^C^, 


h 


= c^w, 


h 3 


= D^D^, 


Ir 




hi 





in which we used the following definitions: 

= W^pCP, (279) 

1. Fundamental Perturbations of FLRW Universe 

As we observed previously, a complete examination of the perturbation theory should naturally include the analysis 
of the evolution of the Debever metric invariants associated to FLRW geometry. 
The only non identically zero invariants of FLRW geometry are given by 



44 



h = (1 - 3A)p, 

7 6 = |(1 + A)V, 

7 7 = -|(l + A)V, 
/ 8 = |i(l + A)V, 

in which we used Einstein's equations, and the stress-energy tensor is that of a perfect fluid. 

If we restrict ourselves to the linear perturbation theory, the only invariants which have non identically zero linear 
perturbation terms are I5, Iq, I7, Is, I9 and I\2- Among these the first four are nonzero in the background and the 
latter two are zero, since the geometry is conformally flat. This could lead to the conclusion that Ig and 1 12 are the 
"good" scalars to be examined. However, a direct calculation shows that the latter two invariants have zero linear 
perturbation. Indeed, it follows from FLRW geometry that the perturbation of Ig reduces to 

si 9 = c^cr^sw^p. 

Then (due to the fact that Weyl tensor is trace-free) the above quantity identically vanishes. This result depends of 
course on the fact that the source of the background geometry is given by a perfect fluid. In effect we have in this 
case 



8i a = {p+pf (v»v v - -<r) \ v a v p - -^J SW^p, 

which is zero. Likewise 5I±2, given by 
also vanishes. 

The corresponding perturbations for the remaining invariants are given by 

8I 5 = (1 - 3X)S P , 
SI e = l(l + \) 2 P 6p, 
8I 7 = -|(l + A) V Sp, 
« 8 = §(I + A)V Sp. 

It follows from these results that the perturbations of these quantities are algebraically related 14 . Besides, once all 
these scalars have a non-zero background value, they do not belong to the minimum set of good quantities that we 
are searching for. 

Corresponding difficulties occur for the standard kincmatical and dynamical variables, that is, the expansion pa- 
rameter 9 and the density of energy p suffer from the same disease. 

This is thus the bad choice for the basic variables which we should avoid. Let us now turn our attention to which 
good variables should be considered as the fundamental ones. 

Geometric Perturbations: From the previous section it follows that 



One can write these invariants in a pure geometrical way without using Einstein's equations. This does not modify our argument. 
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is the only quantity that characterizes without ambiguity a true perturbation of the Debever invariants 15 . We need 
thus to only consider the perturbed £!y since, as we shall see, any other metric quantity does not belong to the "good" 
basic nucleus needed for a complete knowledge of the true perturbations. We then set the expansion of this tensor in 
terms of the spherical harmonic basis 

SEij =E(t) Qij(x k ). (280) 

Thus E(t) is the geometric quantity whose dynamics we are looking for. 

Kinematical Perturbations: We restrict our considerations only to linear perturbation terms. The normalization 
of the 4-velocity yields that the variation of the time component of the perturbed velocity is related to the variation 
of the (0-0) component of the metric tensor, that is: 

SV = ^Sgoo- (281) 
The corresponding contravariant quantities are related as follows: 

SV° = hg 00 = -SV . (282) 
The expansion of the perturbations of the 4-velocity in terms of the spherical harmonic basis is 16 

(283) 

SV k = V(t) Q k {x l ). 

For the acceleration we set 



For the shear 



and for the expansion we set 



Sa k = (t) Q k (x l ). (284) 



San =E(t) Q^(x k ), (285) 



60 = H(t) Q(x l ) + Z(t), (286) 

where Y(t) and Z(t) are homogeneous terms that are not true perturbations. 

Let us point out that, once we are limiting ourselves to the analysis of true perturbed quantities, the important 
kinematical variable whose dynamics we need to examine is only E(£), since the other gauge-invariant quantity \I/ is 
a function of S (and E), as we shall see {j3 is just a matter of choosing of the coordinate system). 

Matter Perturbation: Since we are considering a background geometry in which there is a state equation relating 
the pressure and the energy density, i.e. p = A/5, we will consider the standard procedure that accepts the preservation 
of this state equation under arbitrary perturbations. Besides, our frame is such that there is no heat flux. Thus the 
general form of the perturbed energy-momentum tensor is given by 

ST^ = (1 + A) SipV^Vu) - \8{pg^) + SIL^. (287) 
We write Sp in terms of the scalar basis as: 



This is a consequence of the vanishing of the perturbation of the magnetic part of Weyl tensor (cf. above). 
The vorticity is of course zero, since we are limiting ourselves to the irrotational case. 
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Sp = N(t)Q(x t )+ f i(t), (288) 

in which the homogeneous term /i(t) is not a true perturbation 17 . 

According to causal thermodynamics, the evolution equation of the anisotropic pressure is related to the shear 
through 



rtlij+Uij =C(Tij, (289) 

in which r is the relaxation parameter and £ is the viscosity parameter as we saw previous section. For simplicity of 
this present treatment, we will limit ourselves to the case in which t can be neglected and £ is a constant 18 ; Eq. (289) 
then gives 



Uij=t<ru, (290) 

and the associated perturbed equation is: 

'VII,., OVt,,. (291) 

Following the same reasoning as before, SUij is the matter quantity that should enter in the complete system of 
differential equations which describes the perturbation evolution. One should also be interested in the dynamics of 
6p although it is not a fundamental part of the basic system of equations. We will examine its evolution later on. 

The "good" set Alu] has therefore three elements: 5Eij, Saij and STlij. But, since STlij is written in terms of <5<xy, 
the set M[A] which will be considered reduces to: 

M[ A ] = {6Eij,6aij}. 
So much for definitions. Let us then turn to the analysis of the dynamics. 



2. Dynamics 



In this section we will show that E(t) and E(t) constitute the fundamental pair of variables in terms of which 
all the dynamics for the perturbed FRW geometry is given, that is, A4\a] = {E(t),Y,(t)} is the minimal closed set 
of observables in the perturbation theory of FRW which characterizes and completely determines the spectrum of 
perturbations. Indeed, the evolution equations for these two quantities (which come from Einstein's equations) 
generate a dynamical system only involving E and E (and background quantities) which, when solved, contains all 
the necessary information for a complete description of all remaining perturbed quantities of FRW geometry. Such 
conclusion does not seem to have been noticed in the past. 

We remark that we will limit ourselves only to the examination of the perturbed quantities that are relevant for 
the complete knowledge of the system. These equations are the quasi-Maxwcllian equations of gravitation and the 
evolution equations for the kinematical quantities. In Vishniac (1990) and Novello (1988), this system of equations 
was presented and analyzed. 

The Perturbed Equation for the Shear: The perturbed equation for the shear Eq. (106) is written as: 



ha^hp" (Sa^y+ -6 8<j af} + —h a p Sd -x - -h^hp^Sa^ + 5a v; „] = 8M a0 , (292) 

where 



17 We will set Y = Z = fi = 0, since these homogeneous terms are just a matter of choosing of the coordinate system. Nevertheless we are 
not interested in examining pure gauge quantities such as Y, Z and fi. 

18 In the general case, § and r are functions of the equilibrium variables, for instance p and the temperature T and, since both variations 
STlij and dau are expanded in terms of the traceless tensor Qij, it follows that the above relation does not restrain the kind of fluid we 
are examining. However, if we consider § as time-dependent, the quantity STlij must be included in the fundamental set A^r^j. 



47 



M af) = R a ^ v VV V - ^R„ v V»V V hap. (293) 
Using the above spherical harmonics expansion and Eq. (291), Eq. (292) reduces to: 

E = -E - -£ E + m (294) 

The Perturbed Equation for Eiji The perturbed equation for the electric part of the Weyl tensor is given by 
(101). Using the above spherical harmonics expansion and Eq.(291), one obtains: 

^.(l^) pS .(| + |) £ _|(| + |) s+f{4 . 

This suggests that E and E may be considered as canonically conjugated variables. We shall see later on that this is 
indeed the case. 

Equations (294) and (295) contain three variables: E, E and VP. We will now show that using the conservation law 
for the matter we can eliminate \I> in all cases, except when (1 + A) = 0. We will return to this particular (vacuum) 
case in a later section. 

The proof is the following. Projecting the conservation equation of the energy-momentum tensor in the 3-space, 
that is 

T^-vh^ = 0, (296) 
and using the perturbed quantities this equation gives: 

(1 + X)p Sa k - X(5 P ). k + Xp 6V k + OT fe ^ = 0. (297) 
Using the decomposition in the spherical harmonics basis we obtain 

(1 + A)p* = X[N - pV] + 2£ f I + — ) A- 2 E. (298) 

\3 my 

Now comes a remarkable result: the right hand side of Eq. (298) can be expressed in terms of the variables E and 
E only (since we are analyzing here the case where (1 + A) does not vanish). Indeed, from the equation of divergence 
of the electric tensor — see Eq. (532 — we find 

N - pV = + £ ZA- 2 -2(1 + —) A' 2 E. (299) 
\ m J \ m J 

Combining these two equations we find that ^ is given in terms of the background quantities and the basic perturbed 
terms E and E: 

(1 + A) p * = 2 ( 1 + — ) A- 2 [-XE + iA£E + ^£E]. (300) 
\ my 2 3 

Thus the whole set of perturbed equations reduces, for the variables E and E, to a time-dependent dynamical 
system: 



E = iq(E,£0, 
E = F 2 (E,E), 



(301) 
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with 



i*i = -E — -£ £ + m#, 

and 

in which ^ is given in terms of E and S by Eq. (300). 



3. Comparison with Previous Gauge- Invariant Variables 

FLRW cosmology is characterized by the homogeneity of the fundamental variables that specify its kinematics 
(the expansion factor 0), its dynamics (the energy density p) and its associated geometry (the scalar of curvature 
R). This means that these three quantities depend only on the global time t, characterized by the hyper-surfaces 
of homogeneity. We can thus use this fact to define in a trivial way 3-tensor associated quantities, which vanish in 
this geometry and look for its corresponding non-idcntically vanishing perturbation. The simplest way to do this is 
just to let U be a homogeneous variable (in the present case, it can be any one of the quantities p, 8 or R), that is 
U = U(t). Then, we use the 3-gradient operator V M defined by 



V M = h„ x V A , (302) 



to produce the desired associated variable 



Up = VpU. (303) 

According to Ellis (1989), these quantities were discussed and its associated evolution analyzed. In the present 
section we will exhibit the relation of these variables to our fundamental ones. We shall see that under the conditions 
of our analysis 19 these quantities arc functionals of our basic variables (E and S) and the background ones. 

The Matter Variable \i- It seems useful to define the fractional gradient of the energy density Xa as 

Xa = - V Q P . (304) 
P 

Such quantity \a is nothing but a combination of the acceleration and the divergence of the anisotropic stress. 
Indeed, from Eq. (297) it follows (in the frame in which there is no heat flux) 

SXi = Sdi + j- p SI!/.?, (305) 

From what we have learned above, it follows that this quantity can be reduced to a functional of the basic quantities 
of perturbation, that is S and E, yielding 



+ ^ (M s ) 



The Kinematical Variable r^: The only non- vanishing quantity of the kinematics of the cosmic background fluid 
is the (Hubble) expansion factor 9. This allows us to define the quantity r\ a as: 



Remind that here we restrain our examination to irrotational perturbation. The formulas which we obtain are thus simpler. However, 
the method of our analysis is not restrictive and the study of generic cases can be obtained through the same lines. 
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V a = hj 0,p. (307) 
Using the constraint relation Eq. (4) we can relate this quantity to the basic ones: 




(308) 



The Geometrical Variable t: We can choose the scalar of curvature R which depends only on the cosmical time 
t like p and 6 to be the [/-geometrical variable. However, it seems more appealing to use a combined expression r 
involving R, p and 9 given by 



r = i?+(l + 3A) p- 




(309) 



In the unperturbed FLRW background this quantity is defined in terms of the curvature scalar of the 3-dimensional 
space and the scale factor a(t): 




We define then the new associated variable r Q as 

T a = Tp. (310) 

This quantity r Q vanishes in the background. Its perturbation can be written in terms of the previous variations, 
since Einstein's equations give 



We can thus, without any information loss, limit all our analysis to the fundamental variables. Nevertheless, just 
for completeness, let us exhibit the evolution equations for some gauge-dependent variables. 

Perturbed Equations for p and 8: From Eq.(lll) and using the decomposition of the perturbed energy density 
in the scalar basis (defined by Eq.(288)) we obtain the equation of evolution for dp as: 



N-^Pp+(l + \)9N+(l + \)pH = 0. (311) 

Applying the same procedure for the perturbed Raychaudhuri equation (108) and using the decomposition in a 
scalar basis (286), we obtain 



g_l J 3d+gfl g+ gtt + ( 1 + 3A ) iV = 0. (312) 

To solve these two equations we need to fix the gauge (3(t) and to use the values for E and £ which were obtained 
from the fundamental closed system found in the previous section (Eqs. (301)). All the remaining geometrical and 
kinematical quantities can be likewise obtained. This completely exhausts our analysis of the irrotational perturbations 
of FLRW universe. 
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4- The Singular Case (1 + A) = 0: The Perturbations of de Sitter Universe 



We have seen that all the system of reduction to the variables £ and E was based on the possibility of writing 
the acceleration in terms of E and S. This was possible in all cases, except in the special one in which (1 + A) = 0. 
Although no known fluid exists with such negative pressure, the fact that the vacuum admits such an interpretation 
has led to the identification of the cosmological constant with this fluid. It is therefore worthwhile to examine this 
case in the same way as it was done for the previous sections. 

At this point it must be remarked that, contrary to all previously studied cases, perturbations of this fluid must 
necessarily contain contributions which come from the heat flux or the anisotropic pressure. Indeed, if we take both 
of these quantities as vanishing, then the set of perturbed equations implies that all equations are trivially satisfied, 
since all perturbative quantities vanish, except for the cases where Sp = A Sp, with A — 0, and A + 1 = 0. We will 
analyze these cases below. 

When Sp = A Sp, for A = 0, the system is stable. Indeed, we obtain for the electric part of Weyl tensor, in the case 
that is constant in the background, the following expression: 

E(t) = E e-^. 

The other case of interest is the one in which the condition (1 + A) = is preserved throughout the perturbation. 
Looking at Eq. (311) it follows that, from the fact that p = and reminding the reader that (1 + A) = 0, temporal 
variation of the energy density exists only if we take into account the perturbed fluid with heat flux. We then write 

<z» = q(t) Qi{x k ). 

Equation (111) gives 

TTl 

N = ^q. (313) 



The projected conserved equation gives (see Eq. (112)): 



The evolution equation for the electric part of Weyl tensor gives: 



in which we used the definition 



S+-S = ~q, (315) 



Finally, from the equation that gives the divergence of £!y , we have the constraint 

^(l + ^.)s = -(N + eq). (316) 

The evolution equation for the shear provides the value of the acceleration Equations (314)- (316) constitute 
thus a complete system for the variables E, £ and q. This completes the general explicitly gauge-invariant scheme 
that we presented here even in the singular case (1 + A) = 0. Notwithstanding, just as an additional comment, it 
would be interesting to consider the perturbation scheme in the framework of Lanczos potential. This will be done in 
a later section. 
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5. Hamiltonian Treatment of the Scalar Solution 

The examination of the perturbations in FLRW cosmology, which we analyzed above, admits a Hamiltonian for- 
mulation that is worth considering here (cf. Grishchuk, 1990). In this vein, the variables E and E, analyzed in the 
previous section, are the ones that must be employed to obtain such a formulation. From the evolution equations for 
E and E (Eq. (301)) it follows that they are not canonically conjugated for arbitrary geometries of the background. 

The natural step would be to define canonically conjugated variables Q and P as a linear functional of E and E 
as 20 : 



(317) 



It should be expected that functionals of the background geometry would appear in the construction of the canonical 
variables in the functions a, (3, r\ and 5. It seems valid to remark that this matrix is univocally defined up to canonical 
transformations. We can thus use this fact to choose r\ and 5 as zero; we shall use this choice in order to simplify our 
analysis. 

The Hamiltonian % which provides the dynamics of the pair (Q, P) is obtained from the evolution equations of E 
and E (301). The condition for the existence of such a Hamiltonian is given by the equation 



~Q~ 




a 7] 




' E " 


p 




8 /?_ 




E 



a p 3 



2m£ 



1 



3e 
m 



0. 



3(1 + A)pa 2 

It then follows that the Hamiltonian which provides the dynamics of our problem takes the form 



(318) 



where h±, hi and /13 are defined as 

■(l + > 



hi 



-P 



n = ^Q 2 + ^-P 2 + 2h 3 PQ, 



2 V 2 3 



m£ 2 



(1 + X)pa 2 



3e 
m 



(319) 



(320) 



2mA 



(1 + X)pa 2 



3e 
m 



(321) 



2mA 



(1 + \)pa- 













m J 



(322) 



Let us consider the case in which ^ = 0, that is, there is no anisotropic pressure. The case where £ does not vanish 
presents some interesting peculiarities which will be left to a forthcoming section. 

We will choose /3 = a and take a as given by Eq.(318). We then define the canonical variables Q and P by setting 

Q = E, 
P = aE. 

It then follows that H is given by 



H 



-A 2 (t)P 2 + 7 2 (i)Q 2 , 



(323) 



The attentive reader should notice that in this subsection the quantity Q shall not be confused with the previous scalar basis. 
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where j(t) and A(i) are given in terms of the energy density of the background p, the scale factor a(t) and the wave 
number m as: 

A 2 (i) = 

Let us make two comments here: first of all, the fact that the system is not conservative (which means % is not 
zero) is a consequence of the fact that the ground state of this theory (Q = P = 0) corresponds not to Minkowski 
flat space-time but to FLRW expanding universe. The second remark is that the same applies to the non-positivity 
of the Hamiltonian; this is also a consequence of the non- vanishing of the curvature of the fundamental state. The 
system which we are analyzing is not closed and, so, momentum and energy can be pumped from the background. 

We notice that the Hamiltonian structure obtained in terms of the variables E and E is completely gauge-invariant 
and, as such, deserves an ulterior analysis, which we will make elsewhere. We would like only to exhibit an example 
where this pumping effect can be easily recognized; this will be achieved by applying the Hamiltonian treatment to a 
static model of the universe. 



6. Fierz-Lanczos Potential 

As it was remarked in a previous section, perturbations of conformally flat spacetimes do not need 21 the complete 
knowledge of all components of the perturbed metric tensor 5g^, although they certainly need to take into account 
the Weyl conformal tensor, since all the observable information we need is contained in it (namely, SEij and SHij). 

Let us note at this point that the tensor W a p^ v can be expressed in terms of the 3-index Fierz-Lanczos potential 
tensor — see Ficrz (1939) and Lanczos (1962) — that we will denote by Lap^, and which deserves a careful analysis. 
Indeed, one could consider 8L a p^ as the good object for studying linear perturbation theory, since, as we shall see, it 
combines both <5£.y and 5ak (which are alternative variables to describe SEij). 

Before going into the perturbation- related details, let us summarize here some definitions and properties of L a p^, 
since the literature has very few papers on this matter 22 . 

Basic Properties: In any 4-dimcnsional Ricmannian geometry there is a 3-index tensor L a p^ which has the 
following symmetries: 

Lafifi + Lp ali = (325) 

and 

Lapp + Lp^a + L^ a fj = 0. (326) 

With such L Q/ 3 M we may write the Weyl tensor in the form of a homogeneous expression in the potential expression, 
that is 

2 

+ -^L a a .\ g a j3^v, (327) 



21 The above quoted gauge problem has been widely discussed in the literature (see [120] and references therein). 

22 This tensor was introduced in the 30's to provide, in a similar way as the symmetric tensor i/? M „ does — in a more used approach — an 
alternative description of spin-2 field in Minkowski background. In the 60's Lanczos rediscovered it — without recognizing he was dealing 
with the same object — as a Lagrange multiplier in order to obtain the Bianchi identities in the context of Einstein's General Relativity. 
However, only recently (cf. Novello, 1992a, b) a complete analysis of Fierz-Lanczos object was undertaken and it was discovered that 
its generic (Fierz) version describes not only one but two spin-2 fields. The restriction to just a single spin-2 field is usually called the 
Lanczos tensor. We will limit all our considerations here to this restricted quantity. 



(1 + A) 
— A — 

1 [ n 

1m 



pa, 

2mA 
(1 + A)pa 2 



1 + 3e 
1 T m 



(324) 
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where 



T = T c 
' J a: li — ±-*a 



and 



Let us point out that, due to the above symmetry properties, Eqs. (325) and (326), Lanczos tensor has 20 degrees of 
freedom. Since Weyl tensor has only 10 independent components, it follows that there is a gauge symmetry involved. 
This gauge symmetry can be separated into two classes: 

A (1) L Q(3Ai = M a g Ptl - Mp g afi , (328) 

and 

A {2) L a ^ = w a p-, M -hv ltai p + hv l 0. a + 

+ \g m tV ; a - \am W a x -x, (329) 

in which the vector M a and the antisymmetric tensor W a p are arbitrary quantities. 

Lanczos Tensor for FLRW Geometry: The fact that Friedman-Lematre-Robertson- Walker geometry is con- 
formally flat implies that the associated Lanczos potential is nothing but a gauge. That is, we can write the Lanczos 
potential for FRW geometry as 



L a pn = N a gp^ — N/3 g at i + F a p.^ — -F^ a .p + 

+ ^F^P;a + -flVra F P X - X ~ 2^/3 • F M A ;A' ( 330 ) 

for the arbitrary vector N a and the antisymmetric tensor F a p . 

Perturbed Fierz-Lanczos Tensor: In the case we are examining in this paper (irrotational perturbations) the 
perturbed Weyl tensor reduces to the form 



SWaPnv = (VaPje 1)liv\p - 9af3-ye g t iv\p)V 1 V X 5E ep , (331) 

since the magnetic part of Weyl tensor remains zero in this case. 

It then follows that the perturbed electric tensor is given in terms of Lanczos potential as: 



- SEij = <5L i[o ; j] + SL j[ .q - -^L{QQ)lij 

- \sL {ij) + ^8L° x a . iXlij . (332) 

Although the L a p^ tensor is not a unique well defined object (since it has the gauge freedom we discussed above) 
we can use some theorems — see Novello (1987) and Lopez Bonilla (1989) — that enable one to write L a pn in terms of 
the associated kinematic quantities of a given congruence of curves present in the associated Ricmannian manifold. 
Following these theorems and choosing the case of irrotational perturbed matter, it follows that 5L a p^ (the perturbed 
tensor of FLRW background) is given by 



where 



5L a p^ = Sa^aVp] + F(t)8a [a V 0] V^ 



(333) 



54 



In other words, the only non identically zero components of SL a p^ are: 

6L okQ = -F(f)* Q k , (335) 

and 

SL 0ij = -Jl{t)Qi i , (336) 

which coincides with the previous results. 

From what we have learned in the previous section, we can conclude that this is not a univocal expression, that is, 
Eqs. (335) and (336) arc obtained by a specific gauge choice. 

Let us apply the above gauge transformation to the present case. In the first gauge, Eq. (328), we decompose vector 
M a in the spatial harmonics (scalar and vector): 

M = M^(t)Q(x), (337) 

Mi=M^(t)Qi(x), (338) 

and in the second gauge, Eq. (329), we have 

W 0i =WW(t)Q i (x), (339) 



and 



W v = -\e ijk W {2 \t) Q k (x). (340) 



y ~ a 2 

To sum up, asking what the Lanczos tensor is for the perturbed FLRW geometry is one of those questions (like 
the one about the perturbed tensor <5<7 M „) that should be avoided, since this quantity is gauge-dependent. A good 
question to be asked should be, as we remarked before: What is the perturbation of Weyl tensor? This was precisely 
the motivation of the previous section. 

E. Friedman Universe: Vectorial Perturbations 

As it has been discussed in the previous section, we will make use of the perturbation formalism of Einstein's theory 
of gravitation based in gauge-independent and evident physically meaningful quantities, such as the vorticity, shear, 
electric and magnetic parts of the conformal Wcyl tensor and so on. 

In the scalar case, the convention was simplified in order to facilitate the resulting system of dynamical equations. 
For the vectorial and tensorial cases, however, we feel that the convention set in Hawking (1966) is more adequate. 
Therefore, we will present it here. 

The metric of the background is given in the standard Gaussian form, thus defining a class of privileged observers 
V a = Sq . The projector /i M „ defines, in the 3-dimensional space orthogonal to V a , the 3-dimensional quantities with 
the symbol (~ ). Thus, X a = h a ^ Xp denotes a projection into the 3-geometry. For the same reasoning, we define the 
operator V a as the covariant derivative in the 3-geometry. The relation between the 3-dimensional Laplacian (V 2 ) 
and the 4-dimcnsional one is given as follows: 



V 2 X a = (J\ X a + hi V 2 Xt 
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We then introduce the fundamental harmonic basis of the functions projected onto the 3-surfacc 

{Q(x)} } {P Q (x)}, {U a p(x)}. (341) 
In this section, we are interested in the vector basis P a (x), which is defined by the following relations: 

V" = o, 

P U ' (342) 

Y7'2 p m p 

V r a — -^2 io , 

where the eigenvalue (again denoted by m, despite the fact that this eigenvalue and the scalar basis one have no 
relation at all) is given by 

q 2 + 2, < q < oo, e = +l(open), 

q, < q < oo, e = (plane) , (343) 

n 2 - 2, n = 2, 3, . . . , e = -1 (closed). 

From this basis, it is possible to derive a pseudo- vector and a tensor: 

P a p=V P a , (344) 

P* — Vo P* 

suited to developing pseudo- vectors and tensors. 

The following vectorial relations are useful in obtaining the dynamical equations: 

p* — _ 2 a p* 



V^P (Q/3) = -^(m + 2e)P a , 

V^(^) = ^(m + 2e)P*, ( 345 ) 

V72 p* 771 p* 

p* ;£ = ^ (m - 2e) P Q , 

i.v Xjt y V7 p _ JjA* fti/ p* 
"(a "(9) ''A 1 Vt V 7 r H "(a 'V) /*f 



where we have used the constraint relation below, 



2 



a 2 3 r V 3 

valid in the FLRW background. The following auxiliary relations are also useful 



--oP+lo) =«. ( 346 ) 
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e = -le*-i( P + 3 P ), (347) 

p = -9(p + p). 

With the above basis, we are able to expand any good perturbed quantity (again denoted by SX, where X is any 
quantity associated to the matter content, kinematics and geometry) as 

8uj a = n(t)P*, 

Sq a = q(t)P a , 
5a a = *(t)P Q , 
8V a = V(t)P a , 

(348) 

5<Jap = E(£)P( Q( 3), 

SH aP = H{t)P* {ap) , 
8E aP = E(t)P aP , 

<5n Q(3 = n(t)P M) . 

1. Dynamics 

In order to get simpler equations, we will again approximate the thermodynamic equation, 

rn Q(3 +n Q(3 =^a/3, (349) 
to the limit of small relaxation time r (adiabatic approximation) and constant viscosity coefficient £ to get 

STlap = £ 8a afj ~* n = £ E. (350) 
The vorticity can be written in terms of the 3-vclocity as 

5u a = ~SV a ~» V = -2 0. (351) 

We will denote by (x r , <&s) the fundamental dynamical and constraint equations, respectively. Introducing equa- 
tion (348) and equations (350)-(351) into the perturbed quasi-Maxwellian equations (101)-(112) and making use of 
equations (342)-(345), we get 



^ " ¥ t + \ e E + \ {p + p) E + i (m - 2e) H + \ q = °' ! ' !r,2i ' ] 



E+Q + |j £ + = 0, (352b) 

f2+-0fi+-* = O, (352c) 

H +-6H -]-E -i£E = 0, (352d) 

3 2 4 v ' 

q + -6 q + \ (m + 2e) £E + 2pO + (p + p) * = 0, (352e) 
o fx 



and 
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S + fi + 2H = 0, (353a) 
-1 ( m + 2e) £ - (m + 2e)££ + ^0(p+p)f2-^0g = O, (353b) 

^(m + 2e)H - (p + p)n + ^q = 0, (353c) 

-^(m + 2e)£ + |i(m-2e)+4 (J^J + ^ (p + - q = 0. (353d) 

By making use of equations (346)-(347), it can be easily shown that constraint (353d) is not essential, since it is 
written in terms of (353a) and (353c). We also note that we can write constraint (353b) in a simpler form as 

E-^Z + ^6H = 0. (354) 
The fundamental differential system is now written as 



E - \ i t + ^6 E + 1 (p + p) E + ^ (to - 2e) if + J 9 = 0, ( sr,:,,, i 



^ + |)s + B- i* = 0, (355b) 

fi+i(9fi + i<I< = 0, (355c) 

F+-0If--.E--£E = O, (355d) 

3 2 4 v y 



9 + r <? 9 + -4 (n» + 2e) £ S - 2p + (p + p) * = 0, (355e) 
3 <r 



and 



£ + fi + 2# = 0, (356a) 
£;-ieS + ^6»F = 0, (356b) 
1 

X 2 v 



2 (m + 2e)iJ-( /9 + p)fi + - g = 0. (356c) 



It could be argued that the acceleration *5> should be eliminated from the dynamical system by using the definition 

a a = V a = V a ., V . 

If this is done, we obtain 



= Iv + ^j P a -sr° Qa . 

However, it is easily proven (see Novello, 1995a) that 

K a = ^g 00 ), a = (5V ), a , 
which is zero in the vector basis. Then, making use of Equation (351), we have the following relation: 

* = -2fi- -on, 

3 
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which is precisely the dynamical equation (355c). The variable "J can then be eliminated only by means of losing 
some degree of freedom. This way we get physically motivated {i.e., by observation) algebraic relations between 
acceleration and other selected variables. We will restrict ourselves here to the three cases that follow. 

The first possible choice is to admit a shear-free model for the cosmological perturbation. In such a case, there is no 
shear, and, from this, the anisotropic pressure vanishes too. Thus we will refer to this case as "isotropic" throughout 
the remainder of this section. Therefore (355b) becomes 

^ = 2E. (357) 

The second possibility is to admit that no vorticity should be taken into account. As it has been known for long, 
the presence of a non vanishing vorticity usually brings together along some troubles related to causality violation. 
So, we motivate this case by eliminating the main source of causality breakdown. In this case, we have 

n = o, 

and (355c) then results in 

* = 0. (358) 

Another possibility is to impose the physical source of curvature to be a Stokesian fluid. This means that the energy 
flux (heat flux in this case) vanishes. Despite the fact that we can always set this quantity to zero by a suitable choice 
of observers, it actually represents a true restriction, for our equations are written in such a way that no observer 
changes can be performed — that is, we have already fixed the observer by imposing the particle flux to vanish. Now 
(355e) yields 

fr= -\{m + 2e) r -l- T S + 2A0n, (359) 



(P + P) 



with 



(p + p) = (1 + X) p ^0 X = const, 

a relation that eliminates VP for all but the de Sitter background. All three possibilities will have their respective 
dynamics and Hamiltonian treatment investigated in a later section. 

2. Permanence of Constraints 

Since we obtained a constrained differential system, given by Equations (355) and (356), it is useful to consider 
whether constraints are automatically preserved or not. If one differentiates the expressions (356) and inserts into the 
results the relations (346)-(347), one gets directly 



$2 = XI - § X 4 - 2 (p + p) $1 - f *2 + \ $3 - 



-e , (6\ 2 1 



(Q + 2H), (360) 



4> 



: -0 + p)X3 + ^(ro + 2e)x4 + 2 X5- ^(m + 2e) $ 2 - §0$ 3 - 



where z = 1,2,3 are the constraint Eqs. (356) and Xj> 3 = 1) 5 are the evolution Eqs. (355). Thus, it follows that 
no secondary constraint 23 (SC) appears in the case of vector perturbations. One should expect it, since this result 
reflects the fact that our basic (quasi-Maxwellian) equations are dynamically equivalent to Einstein's field equations, 
which are complete. 



23 Terminology due to Bergmann relative to Dirac's work (1950) on constrained systems. 
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3. Hamiltonian Treatment of the Vectorial Solution 



If we keep all degrees of freedom, as we mentioned before, the simplest solution for Equations (355)-(356) is then 
to consider f as a small arbitrary function of the background — i.e., ^ = ^(t) — which can also be parameterized by 
the perturbation wavelength m. 

The constraints can now be used to eliminate three of the five variables, and the most suited pair for this solution 
is (£, il). The resulting free dynamics is 



£ = 
O = 



f ) s - £n 



(361) 



directly integrable as 



E(t) = a- 2 (t)e^ t I d+ j a 2 (t')e^' 
Q(t) = a- 1 (t) \C 2 - J \a{t')m{t')dt> 



-ffi n(t') + A *(*') 



<2t' 



(362) 



where i?o is the Hubble parameter and cq a positive integration constant. 

Solution (362) can be thought of as a particular case of an arbitrary linear relation 24 between "J and the fundamental 
variables, 



9 = y(t)Q + z(t)P + g{t), 



(363) 



where 



and (Q, P) is a pair of canonical variables (as we shall see) that describe the vector perturbations, given by 



a (3 
7 5 



If 6 ~P 

A \~7 " 



(364) 



where A = a6 — ^ 0. 

The choice for the above variables is motivated by traditional results of perturbations assuming a perfect fluid law; 
within this assumption both the vorticity and the shear are essential variables: none of them may vanish, or else all 
the system turns out to be trivial (see Goode, 1989). In the more general case, such result does not apply. 

Differentiating Equation (364) and making use of the dynamics given above in (361), we get the dynamics written 
in terms of (Q,P) as 



Q = 
P = 



a 



2 

I 9 



E 
S- 



/3-f(a + /3)j + A>-/?)*, 

5-U-r + s)] n + A (7 - 6) *. 



(365) 



To ensure that we arc actually working with canonically conjugated variables, we write the Hamiltonian constraint 



Linearity is a requirement in order to preserve coherence with our basic assumption of linear perturbations approximation. For the 
understanding of the physical meaning of such a relation, see the examples given in Section [III El]. 
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a 



dQ dP_ 
dQ + dP 









P - I (a + 0) 














5- 


7 ~ 


7- 







S - g (7 + 5) 



a+^-(a-0)y+ — (j-6)z 
1 



(366) 



and set the solution of $ = as 



a 3 (t) 
a 



^,a = 6 = A 1 / 2 , 



with ao = const. It will indeed be a solution if 

(a-l3)y+( 1 -5)z = 
holds, which leads to the following three possibilities: 



(367) 



i) y 0, \/z 

ii) y = 0, z 7^ 
Hi) y = z = 



/= |, p = a(l-f), 7 = 0; 



(3 = 0, 
(3 = 0, 



7 = a; 
7 = 0. 



(368) 



For the first case, the dynamics results in 



Q 

2 



-gZ+ea-/)-/-^/ 2 



(369) 
(370) 



described by the Hamiltonian 



a ~2 , 1 



W(Q,P) = I Q 



a 3-?4 QP+ f (0 + / - 9P) - 



(371) 



For the second case the dynamics becomes 



a 6 t 



P 



v 



a 
1 



P 



(372) 



associated with the Hamiltonian 



H(Q,P) 



3 2 



z P 2 



(373) 



The third case is equivalent to the situation given in Equations (355) with new variables, and can be written as 
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= (g-§«-e+f»)o+(- 

P—%vQ + {i-$-%v)p 

The Hamiltonian associated with this case is then 

nQ,P) = -,yQ 2 + l(-l + ^y)p 2 -^-l-^y)QP + ^(Q + P)- (375) 

4- The Specific Solutions 

We proceed to study the three particular cases presented in Section [III E 1] , where the acceleration VP was eliminated 
by an explicit loss of a degree of freedom. 

In the first case (the isotropic or shear- free model), we have E = and, using Equation (357) in the system 
(355)-(356), we obtain the following results: 



p 



a n 

~2~ g ' 



(374) 



H(t) = 


c ia - 2 (t), 


E(t) = 


-^6a- 2 (t), 


n(t) = 




*(*) = 




<?(*) = 


-2C ia - 2 (t) [(m + 2e) a~ 2 (t) + 2 (p + p)}, 



where C\ is an integration constant. It is worth pointing out that a non-zero heat flux is necessary in order that 
a shear-free linear perturbation is obtained, since zero shear is a characterizing condition for no perturbation in the 
perfect fluid case (cf. Goode, 1989). 

The second case (irrotational model, = 0) gives, upon substitution of Equation (358) in (355)-(356), the results 
below: 



E(t) = C 2 a- 2 (t) exp-**, 

H(t) = -Qa- 2 (t) exp-C*, 

1 (377) 

£?(t) = C 2 (f + §) a~ 2 (t) cxp"^, 
g(t) = C 2 (m + 2e)a- 4 (t) cxp-« 4 , 

where C% is another integration constant. In this case, the fluid must be non-perfect in order that a linear perturbation 
with zero vorticity be possible. 

Finally, for the third case (Stokesian fluid), with q = 0, p = Xp and Equation (417) holding, the system (355)-(356) 
gives the reduced dynamics 



E = — 



(m + 2e) 
(P + P) 



n = 



1 (m + 2e) 
2a 1 {p + p) 



E - (1 - 3A) | 



■ 3A) |n. 



(378) 



We again seek a Hamiltonian description with variables (Q, P), using the same transformation given in Equation 
(364). Differentiating these expressions, we find that Equations (378) can be written as 
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Q = la- 



■ e + a i 



1 (to + 2e) 
2^ (p + p) 



1 (to + 2e) „ Q \ 1 



/3- [a(l-3A)+/3(l + 3A)]f} -^(-jQ + aP), 

. l (to + 2e) 



,1 (to + 2e) 
C |1+ 2^ (p + p) 



P = |7- 

+ {s- 

From Equations (379) we read the Hamiltonian constraint 



2»" WTW^ 6 }±^-^ 

[ 7 (l-3A) + <5(l + 3A)]f} ^(- 7 Q + aP). 



(379) 



$ = A 



whose solution is given by 



dQ 
dQ 



dP 
dP 



= A - 



e i 



1 (to + 2e) 
2a 2 (p + p) 



A + A (9 A = 0, 



A(i) = a (1 ~ 3A) (i)e ( + o) 



(380) 



(381) 



where Co is again a positive integration constant. We now set the Hamiltonian variables (Q, P) as given by Equation 
(364) with 

a = S = A 1/2 , 
= 7 = 0, 

where A is given by Equation (381). Therefore we finally obtain the dynamics 



1 



p 



1 (m+2e) 



.2a 2 (p+p) 

submitted to the constraint of vanishing heat flux 



1 (to + 2e) 
2^ (p + p) 



Q-{(1-3A) §} P, 
P 



(382) 



= 



2a 2 , 
(m + 2e) 



P. 



The associated Hamiltonian is then given by 



H(Q,P) 



1 (to + 2e) 2 
2^~Go + pT 3 r 



(1 + 3A)6> 



■(l-3A)^ 2 -{(l + A)^ + f 



1 (to + 2e) 
2^ (p+p) 



Q - 
QP 



(383) 



(384) 



As an example, equations of motion (382) can be explicitly integrated by taking into account Equation (383). Thus 
the system evolution follows 



P=-<|(1 + 9A)- + | 



1 (m + 2e) 



2a 2 (p+p) 



P. 



(385) 



which can be readily integrated, and we finally find 
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Q 



1 



2a 2 
(m + 2e) 



(P+P) 



(1 + 9A) 



/ 

(ff _1 + c o) 



(ro+2e) 



2a 2 (i') U+*)p(*') 



(1 + 9A) 
P = a 2 exp 



(386) 



(H " + c °) 



1 -I- 1 (m+2e) 



d*' 



Returning to the physically relevant variables we particularly find that 



1 (m + 2e) 
2^2 (1 + A)p(0 



eft' 



(ff " +°o) 



(387) 



The perturbation in vorticity appears to diverge — thus breaking down our fundamental approach of the linear 
treatment — for perturbation wavelengths such that 



m < -2e - 2 (1 + A) a 2 p. 



(388) 



However, Equation (343) shows that we always have m > — 2e, and from this fi goes to zero. Such a result could also 
be expected from the angular momentum conservation law. 

Therefore we get the minimal set of observables for the vectorial mode: 



j^vector 



However, the system is not closed, since the variable cannot be written in terms of the other ones. In order to solve 
this system we have then to eliminate one of the variables involved, thus losing a degree of freedom. 

F. Friedman Universe: Tensorial Perturbation 

Here we will proceed as in Section [III El] in order to get an ordinary differential system which describes tensorial 
perturbations in terms of good variables. 

The tensorial basis U a p{x) is defined by the relations 



U afj = 0, 

IV = 0, 
V^ 1 % = 0, 
U a p = Up a , 
V 2 U a p = -gz U a p, 
where the new eigenvalue m has the following spectrum 



(389) 



'7- 



3, < q < oo, e 
< q < oo, e 



n z - 3, n = 3, 4, 
Using the tensor basis, we can define the dual tensor 



hi (open), 
(plane), 
-1 (closed). 



(390) 



(391) 
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The tcnsorial relations below are employed in obtaining the dynamical equations system: 



U. 



a/3 



We remark that it involves the energy density p and the expansion 9, which satisfies: 



(392) 



6> = 3-. 

a 

Its inverse is usually referred to as the Hubble parameter, whose present value will be denoted by Hq. 
Now, we expand the good perturbed quantity in terms of the above basis as 



(393) 



Scrap 




5H a p 


= H{t)U* ap 


5E a p 


= E{t)U aP 







(394) 



where the time dependent functions E, E, H and II are not related to the vector components of the previous section. 

1. Dynamics 

Under the properties (389)-(392) and again making use of relation (350), quasi-Maxwellian equations are written 
as: 



E- -E 
2 



2 



S + \{m-Ze)H = 0, 



H+?-8H + E+^Z = 0, 



constrained to 



r) = H - E = 0. 



We also know that is dynamically preserved as 



(395a) 
(395b) 
(395c) 

(396) 
(397) 



where \2 and X3 are Eqs. (395b) and (395c), respectively. From this we are, therefore, authorized to insert it into 
dynamics. So proceeding we get the unconstrained coupled differential system 



E+ 6 



3 + 2 



(P+P) 



(m - 3e) \ H = 0, 



H 



(398a) 
(398b) 
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It should be remarked that the coefficient H in Eq. (398a) in the de Sitter background yields a positive 25 constant 
leading term, for times such that (1/a 2 ) ~ 0. This feature will be important in Section [IIIF2]. 

We also stress that Equations (398) have no non trivial solution unless both (E, H) are assumed to be non zero. 
That is, both variables are essential in describing tensor perturbations — it should be remembered that these variables 
constitute the electric and magnetic parts of Weyl tensor, allowing one to write "gravitational waves" for "tensor 
perturbations" . 

2. Hamiltonian Treatment of the Tensorial Solution 

The basic system given by Equations (398) can be described in the Hamiltonian language, which provides a more 
elegant interpretation of the dynamical role of our variables. The link between it and perturbation theory has worth 
on its own. We thus introduce new variables 



where we suppose 



i) - ° i) s 



A = dot ( a & I = aS - h ^ 0, 



(399) 



7 8 

which is proven a posteriori to be actually correct. Therefore we can use the set (Q, P) for (E, H) in order to 
characterize the gravitational waves. Inserting definitions (399) into Equations (398) we eventually get 

Q= Lot -a (0+1) -0} E+{$-p (§0+f) - a (I k (t + f) + + p)l -^(m-3e))} H 7 
P = {7-7 + E+{5-s(p + l)- 7 (l, [^(| + |) + ( p + p)]-^( m -3e))}if. 

We also need to show that our variables are, in fact, canonically conjugated to each other, as suggested by notation. 
That is, we again make use of the Hamiltonian constraint, 

A particular solution of Equation (401) is 



dQ dP \3 



and we then set 



A(t) =a 5 (<)e^, (402) 



a = A w , 

5 = A<- 1 - U3 \ (403) 
/3 = 7 = 0, 



where w is an arbitrary constant. 

With the choice (403), and using solution (402), system (400) becomes 



25 The Hubble constant, here translated to 6, is positive from astronomical observations, despite the fact that its magnitude is not 
universally agreed upon. Thermodynamical reasoning ensures the nonnegativeness of parameter £. 
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P= - 



1 



1)« 



p_ A (l-2w)g 



i 



2 (p+p) 

Prom this we directly read the Hamiltonian 



2a' 



(m — 3e) 



A (2a 



P 



(404) 



I A (2"-l) 

2 



p- 



A(l - 2w)Q 2 



u - 1 + w 



PQ. 



(405) 



This result shows that de Sitter (0 = const.) geometry admits a tensor perturbation Hamiltonian of a typical 
harmonic oscillator with imaginary mass, which evidences instability. This is obtained by setting the arbitrary 
constant parameter 



_ 3 {29 + e) 

U ~ 2 (50 + 3£) ' 

in the Hamiltonian (405). We thus recover the well known result of instability of dc Sitter solution. The above 
result also shows, however, that the same remark applies to arbitrary Friedman-like background with no tensorial 
perturbation in anisotropic pressure tensor, £ = 0. In such cases we set u = 3/5 to get 



H(Q, P) 
where (Q, P) are given by 



2 



(P + P) 2"( m 



3e) 



P- 



1/5 



Q 



(406) 



Q = a 3 (*)et«*P 

2 (407) 
P = a 2 (t) off P. 

Summarizing, we have found that there is a complete set of good perturbed variables for tensorial modes: 



M tensor = {^pf}. 

In this case, the system is closed and completely independent of the other modes, due to the linearity of the harmonic 
basis. 

We have also obtained the Hamiltonian formulation for all modes, according to the previous sections, and we 
address the possibility to canonically quantize the cosmological perturbations of FLRW universes. This analysis will 
be discussed with more details in the next section. 



G. Friedman Universe: Quantum Treatment of the perturbations 



In the previous Sections [HID], [HIE] and [IIIF], we have shown how to treat, in a completely gauge-invariant 
way, the evolution of the perturbations of FLRW universes. Besides, we have shown that it is possible to select a 
minimal set of observable quantities to analyze the perturbations of the FLRW universe. We have also shown that 
the complete dynamical system of the perturbed geometry is described only in terms of two quantities: E and S, 
(respectively, the electric part of conformal Weyl tensor and the shear) for scalar perturbations; E and H, (where 
the last quantity is the associated magnetic part of Weyl conformal tensor) for tensorial perturbations. For vectorial 
perturbations, in a more general case, this minimal set should be expanded to include P, S, H and the vorticity $1. 
Now we have completed the classical treatment of these gauge-independent perturbations; it is rather natural then to 
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go beyond this classical theory. Indeed, the purpose of this section is to treat these perturbations into the quantum 
framework. 

This amounts to using a semi-classical description in which the background geometry is taken in the classical 
framework and considering the perturbations as quantum variables. There are many ways to perform this task. 
Here we will follow a very natural way that consists of applying the method of the auxiliary Hamiltonian, which 
was introduced before. The problem can be stated in the following way: using the quasi-Maxwellian formulation of 
Einstein's General Relativity, we find out that the complete dynamical system reduces to the form 



[*]-"[£]■ (408 > 

in which M. is a 2 x 2 matrix containing information that characterizes the background geometry and Mj are the 
"good" observables which describe the perturbations 26 . The dot (') denotes a temporal derivative. 

For the FLRW background, this set constitutes a non-autonomous dynamical system. A direct inspection of the 
matrix M. shows that it is not trace-free; thus, this system does not have a Hamiltonian. Nevertheless, we exhibited a 
method by which we were able to obtain an auxiliary Hamiltonian, %. for this system. As we will see later, the linear 
relation between the associated canonically conjugated variables (Q, P) and the original physical quantities (E, E, 
H or 0) is not unique. This is not a drawback of this approach, but merely a consequence of the fact that the set of 
possible pairs (Q, P) is related in turn by canonical transformations. 

The existence of this Hamiltonian leads us to consider the possibility of employing the canonical method to arrive at 
the quantum version of the perturbed set. The quantum study of these perturbations in FLRW was done by Lifshitz 
(1963), Hawking (1966) and Novello (1983). However, all these previous works deal with variables which either are 
gauge-dependent or follow the general scheme introduced by Bardeen (1980) and subsequent papers (cf. Ellis, 1989), 
which has a difficult physical interpretation. This fact makes the analysis more complex. 

Let us remark that the use of our method, in which the gauge problem is incxistent, seems to be really the best 
way of making the transition to the quantum version. Indeed, alternative methods, for example, the minisuperspaces 
approach (e.g. Ryan, 1991) suffer not only from the need to fix a gauge, but also from the fact that the order in 
which this choice is made (before or after the quantization) leads to different theories. Even the schemes that consider 
gauge-independent variables (as defined by Bardeen initially) have, as a main problem, the fact that these variables 
do not have an evident physical interpretation — which makes the physical comprehension of the results obtained 
quite complex. Since the gauge-independent variables in our method are observables, and completely equivalent to 
Bardcen's variables, the advantages of quantizing the dynamical systems obtained in our procedure are evident. 

In order to perform the quantization of our system we will make use of squeezed states of the Quantum Optics, 
which was first employed in the framework of Cosmology by Grishchuk (1990), Schumaker (1986) and Bialynicha- 
Birula (1987). However, the use of these approaches suffer from the same difficulties pointed out above. So, the 
advantages of incorporating to this approach the method of gauge-independent variables are again obvious. 

All the definitions and notations that were employed to obtain the results for gauge- invariant, observable perturba- 
tions in the FLRW background (scalar, vectorial and tensorial)are equally valid here. So, we write FLRW geometry 
in the standard Gaussian coordinate system. 

For a comoving observer (one with V a = Sq ), we denote p as the energy density, p as the isotropic pressure, and 9 
as the expansion. The constraint relation 



el 



^ 3> + UJ = ° (409) 
holds, as well as the following auxiliar relation, 

p = Po a- 3 ^\ (410) 

which comes from Raychaudhuri equation for a fluid with the usual linear state equation p = A p (with A as a constant 
defined in the interval [0, 1]). The parameter po denotes the value of the energy density for a(t) = 1. 



We are considering here only the cases in which the minimal closed set of observables contains only two variables. In the more general 
vectorial case, M should be a 4 X 4 matrix, as stated above (See Section [IIIG3] for more details). 
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We will denote the viscosity by £. As before, it will be taken in the limit of small relaxation times for the adiabatic 
approximation of the thermodynamic equation (see Novello, 1995a) and, due to this choice as well as to thermodynamic 
considerations, it will be taken as a negative constant. 

As for the conventions, let us point out that we will choose the geometrical units system, h = k = c = 1. 
The constant m will denote the wave number associated to the perturbations in the FLRW background and the 
arbitrary integration constants k and b will be employed throughout the whole section, for all three perturbation 
types. Additionally, we will denote by calligraphic letters the matrixes (such as, for example, the Hamiltonian matrix 
H) and by capital letters their corresponding linear counterparts (for example, the Hamiltonian H). It should be 
remarked at this point that, when following the standard quantization procedure, we will drop any extra indication 
(such as the Hamiltonian H turning into the Hamiltonian operator H), so as not to make our notation excessively 
complicated to be read, except the case of creation a) and annihilation a operators where we will maintain the " (^) " 
symbol to avoid confusion with the scale factor a(t). 

1. The Auxiliary Hamiltonian 

Let us consider the linear two-dimensional dynamical system for the variables M\ and M 2 , which gives the dynamics 
of the evolution in the FLRW universe background for the minimal closed set of these gauge-independent linear 
perturbation quantities given by 



Mi 
M 2 



M 



Mi 
M 2 



(411) 



in which M is a 2 x 2 matrix that may depend on time through the known background quantities and Mj are the 
observables that form the minimal closed set which describe all the perturbations. As it has been pointed out, the 
variables (Mi, M 2 ) are not canonically conjugated. 
We thus define a new set of variables (Q, P) as: 



(412) 



~ Q~ 


= s 


' Mi " 




'a P~ 




' Mi " 


p 




M 2 _ 




7 8 




M 2 _ 



where a, /?, 7 and 8 are functions of time and S, the transformation matrix, has a determinant given by 



A = det(S) = «(5-/37^0, 



(413) 



and its inverse exists. 

As a consequence, the variables (Q, P) satisfy the following dynamics: 



P 



= H 



in which H depends on time through M. and the transformation matrix S. From Eqs. (411) and (412) it follows that 

n = sms- 1 + ss-\ 

If we require (Q, P) to be canonical variables, then the matrix H must be traceless: 



From the equation above, it then follows: 



TrU = 0. 



TrM +— =0, 
A 



(414) 
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which can be easily integrated. 

Thus we have a set (Q, P) of canonical variables, which has an associated Hamiltonian H that is linearly related 
to H. The above condition ensures that the set (a, (3, 7, 8) has only three independent quantities. These degrees of 
freedom are fixed by the canonical transformations, as it will be discussed later. 

The most general quadratic Hamiltonian for our system can be written as 



h P 2 + 2h 3 PQ, 



(415) 



whose matrix form, Ji, is then given as 



H 



2h 3 h 2 
-hi -2h 3 



(416) 



where the hi are given in terms of the functions a, /?, 7 and 5, as well as of quantities in the FLRW background. 
Since "H is a 2 x 2 traceless matrix, we can decompose it as H — jl . <j where fl has the following components: 

(fa - hi ) i 

2 , 2 { h i + h V , 2h 3 

The hk {k = 1, 2, 3) depend indirectly on the parameter t, through the known quantities of FLRW background. The 
vector a is built with the Pauli matrices: 



CTl 



1 

1 



02 



-i 

1 



CT3 



1 

-1 



Let us remark that the hk (written in terms of background quantities) vary according to the perturbation type. They 
will be presented in later sections. We then have: 



Mi 




E 




E 




E 


Ah 




E 









H 



for scalar, vectorial and tensorial perturbations respectively 27 . 

Let us make a canonical transformation by changing variables Q and P into Q and P: 



Q 
P 



= J 



(417) 



We then obtain the new Hamiltonian for the transformed system as a function of the previous one and the transfor- 
mation matrix J, that is, 



U = J % J- 1 + j J' 1 . 



(418) 



The guarantee that the system will keep being described by a Hamiltonian is provided by imposing that H be traceless 

Tr(j J- 1 ^) =0, 

that is, 



detJ=l. (419) 

This is nothing but the well known fact that quadratic Hamiltonians constitute the equivalence class of the harmonic 
oscillator. The group 5'L(2,IR) describes the canonical transformations on a plane. 



The observables for vectorial perturbations give a reduced dynamical system in the specific case of a Stokesian fluid (for more information 
on that issue, see Novello (1995b)). 
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2. The Scalar Case 

The Auxiliary Hamiltonian: In this section we will state the results of the auxiliary Hamiltonian method for 
the scalar perturbations of FLRW background (cf. Novcllo, 1995a, b). The resulting dynamical system for scalar 
perturbations in the general case (with non zero viscosity £) is given by: 



2m 



(1 + A)pa 2 



3e 



A 1 

2 + 3 



2toA 



(1 + X) P a 2 



3e 



E. 



(420) 



and 



E = - 
2 



2m 



A 1 

2 + 3 



e 

2 



;A 



(1 + A)pa 2 V ' mj \2 ' 3 J " 2 v ' " ,r 3 
The components of the Hamiltonian matrix (416), hk, are then written as 



(1 + X)pa 2 



i 3e \ <■ e £ 

ml 6 2 



E. 

(421) 



(l + A) 



i 

2 
1 

A 

2 
1 



p- 1 2 {l + \L)- 1 5- + 



2 1 5L + I |(1-AL) + . 

a^-d/3-a 2 (l + AL)-a/3-+/3 2 ^ + ^ 



2a/3L A 



(|(1-AL) + 



-a^ + /j7~/?(5 (1 ^ A) p + a 7 (l + A£)+Q^^ + 



P5 ( |(i-AX) + 



or 

a<5(l + AL)+/3 7 ( — -(1-AL; 



(422) 



The auxiliary quantity L, which appears in the above equation, is defined as follows: 



L = 



2m 



1 



3e 



(1 + X)pa 2 

From Eq. (423), the condition for the existence of the Hamiltonian, Eq. (414), now reads: 



A 8 . ( L 
A-S-U 1 -! 



0. 



(423) 



We will tackle the simpler case of zero viscosity for scalar perturbations in the FLRW background, which then 
yields: 



A(t) = Ka(t), (424) 

where k is an integration constant. 

Canonical Quantization: The problem to be analyzed at this point is nothing but a single harmonic oscillator 
problem with a time-dependent quadratic interaction. This problem appears in many different places, e.g. the 
equation for quantum test fields in homogeneous and isotropic expanding/contracting universes, Quantum Optics, 
etc. There are many ways to face this problem; here we will follow the standard procedure of Quantum Optics. The 
creation and annihilation operators a and are defined in the standard way as 
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6 t = -^(Q-iP), (425) 

[a, at] = i. 

Using Eqs. (415), (423) and (425), the Hamiltonian then becomes: 

H = Ho+H int , (426) 

with 

U = u>{t) (1 + 2JV), 

(427) 

w(f) = i (/n + fta), 

where JV is the number of particles (photons): TV = a' a. The self- interaction Hamiltonian is given by 

Hint =v(t) a 2 + rj*{t) (a^ 2 , (428) 

where 

V(t) = \(hi- ha) - i h 3 , (429) 

and rf is the complex conjugate of r\. 

Schrddinger equation is then easily written for the operator H as 

i — — — =Hip{x,t), 

with the wave function ip(x, t) given as 

ip(x,t) = U(t,t ) i>(x,t ), 

and U(t,to) is the evolution operator. 

Now we will proceed to solve the equation above by employing the formalism of Quantum Optics. It involves writing 
the time evolution operator as a product of the rotation and the single-mode squeeze operators, along with a phase 
factor: 

U (t, t ) = 5( r ^) R(r) , 

where is a time-dependent phase and the rotation operator -R(r) and the single-mode squeeze operator S^ip) are 
defined respectively as: 

-ir d d* 



where T, r and ip depend of time through the known quantities in the FLRW background and are defined as the 
rotation angle, the squeeze factor and the squeeze angle, respectively. It should be remarked that all these quantities 
are real. For further details and explanations the reader should see Novello et al. (1996). 
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Schrodinger equation thus reduces, after a direct but somewhat long calculation, to the following first order coupled 
differential system: 



0(t) = £0(t), 
P _ 2a; 

cosh(2r) ' / 43Q s 
f = ^ (hi — h 2 ) sin (2<^) — 2h% cos (2y>), 

r(p = \ (hi — h%) cos (2(/j) + h 3 sin (2y) — a>(t) tanh(2r), 

where uj(t) is defined by Eq. (427). 

We proceed now to solve the above system, Eq. (430). The last two equations are coupled, which makes their 
integration very much involved. However, this task can be easily accomplished if we perform a transformation in the 
Hamiltonian matrix Ji, by means of a canonical transformation matrix J , given by Eqs. (418) and (419), such that 
the transformed Hamiltonian has now: 



Q=-(h 1 + h 2 ) = ) (431) 

and we will drop the (~) symbol from now on, to keep our notation simpler. This choice can always be made, since it 
amounts only to making the right choice of the original functions a, (3, 7 and 5. The first step taken is then to choose 



where k is an integration constant that comes from the condition of existence of H, Eq. (423). 
The condition (431) may then be written as 



(1 + A) /Kfl\ 2 2 

which then gives 



^) = 5 ^{^(T^<™^ e ,^)}-', (433 ) 

where po is the matter density when a(t) = 1. The above quantity is always positive (since m > — 3e in all circum- 
stances). In this case we have the following results: 

, K 2 {l + X)p _ 2 , 

h ±= 2a a 

h 3 = \ % (434) 

r,(A - l (h h \ ih - k2 ( 1 + x )p i a 
r](t) = - h 2 ) - ih 3 - — ^—2 ^ 5 . 

Proceeding now to integrate the differential equations in r and ip, we now obtain: 



re -2i V = _ 2i / v ^ dt 



To make the integration above simpler we will choose the case of 3-curvature zero and A = 1/3 — radiation era, the 
regime at the quantum phase of the Universe (cf. Novello, 1995b) — , to find 
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v (t)dt = -- 



2 to , 
-r + -: In 



Poo a 



+ - < In a(i) + - In ( a 2 + 



4 po 
2po 



2po 



2po 



+ 



(435) 



With the above result it is easy to decouple the differential equations to find 



/„ x rx 1 / to ~ to , 

rsm(2</>) = -v2/6<- A / a 2 + H In 

a V 2p 4p 



2po . / 2p 



(436) 



r cos (2(p) 



iln 

4 



Observables: From the above construction it follows that the observables of the theory are written in terms of the 
corresponding creation and annihilation operators in the same modes. 
For the shear £ and for the electric part of Weyl tensor E we have 



£ = X(*)fi + X*(*) a\ 
E = a + ^*(t) a\ 



(437) 



where x(t) an d ^(*) are defined as 



X(<) = 



V2A 



V + iP), 



(438) 



*(*) 



v/2A 



(-7 + ia), 



with a, ft, 7, (5 the same quantities defined on Equation (412). In terms of the same solution, Eqs. (432), (433), we 
easily find that 



x(t) 



*(*) 



1 

V2a 



V2h 



(439) 



The commuter is then easily calculated to give 



[S, E] = -i 



na{t) 



(440) 



if the choice H = c = k — 1 holds 28 . 



In the general case we have: 

P = (teV )- 1/4 ^f 

in order that (Q,P) be adimensional. The commutator is then written as 
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The total noise of the observables E and E can be calculated as 29 



AE| 2 |V>) = cosh(2r) (0||AE| 2 |0), 
<^||A^| 2 |?A) = cosh(2r)(0||A£;| 2 |0), 



(441) 



where (Vll^^l 2 ! 7 /') is the total noise calculated in the time t, (0||AX| 2 |0) is the total noise in the vacuum state and 
r(t) is given by Eq. (437). 

3. The Vectorial Case 

In this case we find that the resulting dynamical system is not closed: it depends on a choice of the perturbation 
in the acceleration >J . Three different reduced dynamics were studied: 

• Stokesian Fluid: q = 0; p = Xp; 

• Shear-Free model; 

• Vorticity-Free model. 

The second and third models give very simple, directly integrable results. The Stokesian Fluid model implies, for 
the perturbed acceleration Vf, that ^ = 2X8 £1 and we obtain a closed reduced dynamics for the observables E and 
(respectively the shear and vorticity). This is the case that will be quantized here. 

Auxiliary Hamiltonian: As for the scalar case, the results of the auxiliary Hamiltonian method will be given 
here. 

The reduced closed dynamical system for the special case of Stokesian fluid with non zero viscosity and zero heat 
flux q is: 



i 



(m + 2e) 
2(1 + X)pa 2 



E-(l-3A)Sfi, 



^2(T^ S -( 1 + 3A ^ 



(442) 



where we must have q — 0, p = Xp and ^ = 2X8 Q. The condition for the existence of a Hamiltonian is then written 
as 



A(t) = K[a(t)} 3 ^ e« M «, 
where k is again an integration constant and M (t) is an auxiliary quantity, defined by 



(443) 



M{t) 



(m + 2e) 



2(1 + X)pa 2 (t) 



dt. 



We can now proceed to the same quantization formalism described by Eqs. (425)-(430). Making the choice 



/3 = 7 = 0, 
s A(t) 



4 _ _3£ (m + 2e) Ka (7+9A) ^Mjt) 
2 (1 + A)(1 -3A) p 8 



(444) 



This quantity is denned as the mean square uncertainty in the annihilation operator a. The total noise of a Gaussian Pure State (GPS) 
is conserved even if the total number of photons is not, and, it therefore is more useful to describe the quantum wave functions that are 
obtained from Schrodinger equation. See Appendix for more details on this issue. 
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we obtain the Hamiltonian coefficients hj, (j = 1, 2, 3) as: 



(445) 



It should be noted that, since the viscosity £ is negative due to thermodynamic considerations, and the fact that 
m > — 2e in all cases (see Novello, 1995b), the function a(t) is real. 

We are then able to decouple the differential system that results from Schrodinger equation, Eq.(430), to find: 



r cos {2tp) = In (b ^) , 



(446) 



with a given in Eq.(444) and b being an integration constant. 

Observables: The same method applied for the scalar observables can be employed here in the vectorial case, 
giving (for the choice h = k = c = l) 30 : 

(447) 

and, in this case, the commutator between the perturbed variables gives the following result: 

[E, fl] = i 



Finally, the relation between the total noises of E and Q in a time t and their respective values of total noises in 
the vacuum state is given by 

(*||AE| 2 |tf) = cosh(2r)(0||AE| 2 |0), 
(*||Afi| 2 |W> = cosh(2r)(0||Aft| 2 |0), 

where r is given by Eq.(446). 

The Case £ = 0: If we consider the special case of vectorial perturbations in a Stokesian fluid with zero viscosity 
(£ = 0), we can follow the same steps detailed in previous subsections to find the Hamiltonian coefficients as 

hi = -h 2 = 0, 

(448) 

h _ 1 fa 2, x 



30 For the general case, the perturbations should be written as 

S = (^) 1/4 vk^' 

1 fc 9 \ 1/4 ae« M « 
and the commutator [£, Q] should then be accordingly rewritten 
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where we have the choice 

a = arbitrary function of time, 
p = 7 = 0, 

and the following condition must hold: 

A = i (449) 
so that u>(t) is zero and the system which arises from Schrodinger equation, Eq.(430), is easily decoupled to give: 



r = In 6 



W) / (450) 

if = 0, Or (yS = — 7T, 

with 6 = const, again. Let us note that the condition on A, Eq.(449), ensures that our model of vectorial perturbations 
in a Stokcsian fluid with zero viscosity applies only to the radiation era, as expected. 
The observables £ and f2 are then written as 

(a + a f ), 



a{t) * y/2a(t) 

(451) 



if 7i = fc = c = l 31 . The commutator between S and £1 will then be 

[E,fi]=-i-^, (452) 
and the total noises of the observables in terms of their respective values for the vacuum state will be written as 



<*||A£| 2 |v]/) = cosh ( 2 In ( fc^W ) ) (0||A£| 2 |0), 
(*||Afi| 2 |*) = cosh ( 2 In ( ) ) (0||Afi| 2 |0). 



(453) 



31 In the general case we have 



with the commutator [£, Q] modified accordingly. 



^) 1/4 75fe (a + at) ' 



n = -i fe(ric 3 )^ 1 / 4 ^ffi 4 (a - at), 
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4- The Tensorial Case 

In this case we obtain a new dynamical closed system on the observables E and H (respectively the electric and 
magnetic parts of Weyl tensor), which follows as: 



E 



H = -E 



e ( 6 + £ 



(p + p)] -4,(m-3e)} ff, 



(454) 



The transformation for the variables (Q,P) is made in the same way as in the scalar and vectorial cases, Eqs.(412)- 
(416). The condition for the existence of a Hamiltonian is then written as 



A(t) = Ka 5 (t) e ft , 



(455) 



where k is again an arbitrary integration constant and H is again given in the form of Eq.(415), with the following 
coefficients: 



h 1= ^ < -76 + 5 2 + § 1 5 + 1 5 -Y 



|(f + |) + |,P + P ,-(=^M 



1 



-afi + P 2 + ?,aP + aP- a 2 



13 = JK {~afi-h-0S- 



a d 



■M ( 



(456) 



-a 7 



We are then able to perform the quantization by employing the standard method described by Eqs. (425)-(430) 
and by making the same choice oj(t) = in order to decouple the first order differential system, Eq. (430). We then 
find that 



a 4 (t) = k 2 a 10 e 2 ^ 
/3(f) = 7 (f) = 0, 



(§ + §) + l(, + ,)-^ 



(457) 



S(t) 



a(t) 5 e € *. 



Therefore, the coefficients of the Hamiltonian for the choice given by Eq. (457) are: 



hx{t) = ^a{tfe^ = -h 2 {t), 



h 3 (t) = 1 (g 



(458) 



The decoupled first order differential system then gives the following results: 



r cos (2ip) = \n(b^r e&) , 



r sin (2y>) = — J. 



$(j + $) + $(, + „-i»-*l 



1/2 



(459) 
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where b is an integration constant. We then have 32 : 



E = ^ Q= vfe (a + at) ' 



(460) 

f (a-at). 

and therefore the commutator between the above observables is given by 



H=% a(t)~ 5 e-e* P = -i -% a(t)~ 5 e"^ (a - at). 



a(tr 5 . (461) 

At 

Finally, the total noises for i? and -ff in terms of their respective total noise values for the vacuum state are written 
in the same way as before: 

(*||A£| 2 |1<) = cosh(2r) (0||AS| 2 |0), 

(462) 

(*||Ai7| 2 |*) = cosh(2r)(0||AiT| 2 |0>, 

where the function r is given by Eq.(459). From the above results it follows that the total noise in a time t is always 
greater than its corresponding vacuum value, increasing with r. 

H. Milne Background 

There is a particular class of FLRW geometries dealt with in this present section that seems worth to be explicitly 
examined. This is the case analyzed by Milne and contains a portion of Minkowski geometry. The metric is then 
FLRW-type, where the radius of the universe, the 3-curvature and the expansion are given respectively by: 

a(t) = t, 

£ = +1, (463) 



We will present only the results which follow: 



3 
7- 



Scalar Perturbations 



If wc consider the case of scalar perturbations, the vorticity should vanish, which implies that the magnetic part of 
Weyl conformal tensor will also be zero; thus we have: 



Suiij = 0, 



SHij = 0. 



(464) 



32 In the general case, Eq. (460) can be written as 

E - 



1 (c^\ 1/4 

J? KIT J ' 



H = -i/ifc 3 / 4 «£-4(a-at), 
\/2Ka 5 v ' 



and the commutator [E, H] is suitably altered. 
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Following the notations used before, the other perturbed quantities are listed below: 
Geometric Quantity: 



Kinematic Quantities: 



6Eij =E(t) Qij(x). 

5V = -SV° = \5g m = \p{t)Q{%) + \Y{t), 
SV k = V(t) Q fe (f), 
Sa k = *(t) Q fe (f), 
5(Tij = E(t) Qij(x), 
59 = B(t) Q{x) + Z(t). 

Sp = N(t) Q(x) + L{t), 

SILij = £ ScTij = £ E(i) Qjj(x), 

(5p = A 5/5, 

where we have used again the proportionality relation between the perturbed anisotropic pressure and the shear; we 
also consider the standard formulation in which the perturbed pressure is proportional to the density. The quantity 
j3(t) is gauge-dependent and Y(t), Z(t) and L(t) are homogeneous terms. 

Making use of the quasi-Maxwellian equations we obtain the system for the above quantities as: 



Matter Quantities: 



£6 £9 m 
E = — - + -E+ — E + — q = 0, 
2 3 6 2 1 



(465) 



29 2 
"3" 







E - -E 


g- 






m J 


2 



N + 6q = 0, 



(466) 



• 26> 2 „, . (9 2 (1 + 3A) 



(467) 



E + £ + |E - m * = 0, 



(468) 



(469) 



AT + (1 + A)6W - — g = 0, 
9 



(470) 



q + 8q - XN — 



2£0 2 fl e 



9 V 3 m 



- + — E = 0. 



(471) 
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The dynamical equations on the homogeneous terms Z(t) and L{t) are written as: 



L + (1 + X)9L = 0. (473) 

Let us solve this system for the special simple case where q = 0. We have then (from Eq. (471), the dynamical 
equation for q), 

- \N - ^ ( \ + ±) E = 0. (474) 

Equations (463) and (470) give 



9 V 3 to 



N(t) = N a r 3{1+X \ (475) 
where No is a constant. From (474) and (475) we obtain 



AiVn ( 1 e N 



s «>=-iH5 + y < ~" +w - (476) 



These results applied in equation (466) give 



3A\ /l e N 



= ( 1+ fJ(i + sJ r<1+3A> - < 477 ' 

Equation (465) is automatically valid if we make use of the above results for N(t), E(t) and Equation (468) 

then gives *ff(i) as: 



2m V 3 to 



A(1 + 3A) j (2 + 9A) 



Z 6 



-(1+3A) 



(478) 



It must be noted that the constant A^o cannot be zero, since this would give a trivial result. Equations (467) and 
(469) give the quantities B(t) and V(t) in terms of N(t), #(t) and E(t), B{t) respectively. Both quantities may be 
obtained if the gauge-dependent function j3(t) is chosen. They are therefore "bad" quantities to analyze. The minimal 
closed set of quantities for perturbations in Milne universe is 

M ^alar = { E ^, N, $>} . 

The homogeneous part of (Sp), L(t), is directly determined by equation (472): 

L(t) = L r 3(1+A) , (479) 

where again Lq denotes a constant. The function Z(t), whose dynamics is given by equation (472), can only be 
integrated by choosing another homogeneous term (Y(t)). That completes the solution for the case q = 0. 
We can analyze the behavior of the solution above for different values of A. The results are as follows: 

l.A>-i: 

E, E, TV and \& go to zero when t — > oo; 



81 



2. A = -i: 

E, E and & are constant; A" goes to zero when t — > oo; 

3. -1 < A < -f: 

_E, E and ^ diverge when t — > oo and N goes to zero; 

4. A = — 1 (vacuum A): 

E, E and VP diverge when i — > oo and A" is constant; 

5. A < — 1 (unphysical situation): 

E, E, A" and IP diverge when t — > oo. 

2. Vector Perturbations 
In this case the original dynamical system, Equations (355)-(356) give 

E - | E + | £ + (m - 2e) if + J g = 0, 



E + £-^* = 0, 
9+{j0g+4(m + 2e)£E = O 



Si+|o + ^* = 0, (480) 



and 



E + fi + 2i/ = 0, 

£-|e + §0F = O, (481) 
■4 (m + 2e) ii + A g = 0. 



We will present here only the three cases dealt with in Section (HIE 4): isotropic, irrotational and Stokcsian fluid. 
The results are as follows: 

Isotropic Model: For E = 0, we obtain 

E{t)=ixt- 2 
Hit) = %t~ x 

1 (482) 
#(i) = 2 fit' 2 

q(t) = ~(m + 2)fit- 3 , 

where /i is an integration constant and we used e = +1. These functions of t diverge when t goes to zero and become 
null for infinite values of t. 

Irrotational Model: In the case £1 = 0, the acceleration is also zero and 



E(i) = vt- 2 exp"£*, 

E(t) = v exp-f* t- 2 (I + 
H{t) =-%t~ 2 cxp-« 4 , 
q(t) = v(m + 2) t~ A cxp-« * . 



(483) 
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These functions also diverge when t goes to zero and become zero when t goes to infinity. 

Stokesian fluid: If we consider q = 0, the only possible solution is trivially zero. We conclude therefore that vector 
perturbations in Milne universes must have a non zero heat flux. 



3. Tensor Perturbations 



The original Equations (395) give a closed dynamical system in the variables (E, E): 



E + 
E + 

where H is given by the constraint 



E 



t (0 



+ E + £ = 0, 



-tt (m — 3e) 



S = 0, 



(484) 



£ = H. 

We have then the following set of good quantities for tensorial perturbations in Milne background: 

M A = {E,H}. 

I. WIST Model: Scalar Perturbations 

This section deals with the dynamical system of the perturbed quantities that are relevant for the complete knowl- 
edge of the system. The equations completely describe the perturbation evolution, according to the quasi-Maxwellian 
equations of gravitation. 

In this case, we will assume that the background of the model can be characterized by a source consisting of a 
scalar field minimally coupled to the gravitational field. The energy-momentum tensor for a minimal coupling scalar 
field is represented by a perfect fluid and it can be demonstrated that the general linear perturbations of this fluid 
also behaves as a perfect fluid. This property of the source simplifies the equations and the minimal set of observables 
that determine the scalar linear perturbations of the model can be obtained from the following equations: 



{E al3 j+E 



a/3 



(485) 



1 12 
(5ap V )'+ -h^(8a a ). a - -8a,( p . v ) + -68a ^ 



SEuv, 



(486) 



(p + p)Sa M = (Sp),fih*-p,p6(h%). 



(487) 



The equations are considered to be linear in the perturbations; so, to solve them we will split the perturbations in 
terms of the scalar spherical harmonics basis defined by Eq. (274) in terms of the conformal scale factor a(rf). Since 
the model we are investigating has an open three-section, the eigenvalue m can assume the following values: 



m = q + 1, 0<<7< oo. 



(488) 



With the scalar function Q, we can construct vector Q a and tensor Q^ u quantities. Using this base we can expand 
the perturbations as follows: 



(489) 
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5<T ia/ = ^2a^0^, (490) 
5a^=^V( g ) (?) ^, (491) 

1. Dynamics 

From now on we will suppress the indices q to simplify the notation. Substituting this decomposition in Eqs. (485) 
and (486), after simple algebraic calculations (see Novello, 1995a), we obtain the following dynamical system for each 
mode of the variables E and S: 



S' 



— (3 + m) — a 
pa 



E, (493) 



E' = -—E-paS. (494) 



a 

The prime means covariant derivative projected on V a . 

In terms of matrix, this dynamical system can be written in compact form as 

e)' = m (e)> ( 495 ) 

The components of the matrix M. are: Ain = 0; M.12 = —a + (3 + m)/(pa); M.21 = —p a ', M12 = —a' /a. 

2. Hamiltonian Treatment 

The examination of the perturbations of Robertson- Walker geometries, using the variables associated to the per- 
turbed metric tensor 8g^ Vl admits a Hamiltonian formulation. In this vein, it was shown in detail by Novello (1995b) 
that the present formulation using variables E and E also admits a Hamiltonian formulation. The usual way to do 
this is introducing auxiliary field variables as 



o) (7 5) {I 

The matrix S with components a, /3, 7 and 5 is univocally defined up to canonical transformations. 
In terms of the auxiliary variables, the dynamical system in matrix notation becomes 

q) =s (q 

where 



(496) 



(497) 



S = SMS' 1 - S'S' 1 . (498) 
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The requirement that (Q, P) are canonical variables implies the necessary and sufficient condition that 



A' 



trE = trM + = 0, 
where A is the determinant of S. The Hamiltonian in our case is given by 



(499) 



H = =^P 2 - ^fQ 2 ~^f{PQ + QP). 



(500) 



It seems worthwhile to remark that the matrix M. and S univocally determine H up to canonical transformations. 
This freedom can be used to simplify our analysis in each particular case. 

The background model we will investigate is asymptotic flat in the limits of the conformal time r/ — > ±oo. Convenient 
canonical transformations can be used to fix the functions of the matrix S, in order that in the limit r\ — > — oo we 
obtain, for each mode m, the Hamiltonian of an harmonic oscillator in a Minkowski spacetimc. With this choice, the 
Hamiltonian that describes the system is given by 



H = -P 2 + - 
2 2 



! tanh(2?y) + (m - 7) + 



Q 2 + [tanh(27 ? ) + 1]{PQ + QP). 



(501) 



Using Hamilton's equations, after a simple calculation, we obtain 



,(»/) + [q 2 - 3(tanh(277)2 - l)]Q m ( V ) = 0, 



(502) 



= I 2 - V eff- 



The exact solution to Eq. (502) is given by 



(503) 



Q m (v) = AF{a 1 ,b x ,c 1 ,z) + Bsinh(2r/)F(a 2 , & 2 , c 2 , z)], 
where F{ax, b\,c\, z) and F(a 2 , b 2 ,c 2 , z) are hypergeometric functions with parameters 



(504) 



The variable z is given by 



a\ = | + ^Vm - 1, a 2 



5 
1 



Cl 



1 

2' 



The asymptotic limits 7/ — > ±oo are given by 



c 2 = §. 



z = — sinh 2 (2?7). 



bi = | - ttV^X b 2 = § - iV^T, 



(505) 



Q£fa) = lim Q m (ri) = -= 



'(t?) = lim Q m ( V ) = ^He-'HiHe'". 

T/— fCXD 



(506) 
(507) 



The amplitudes d\ (m) and d 2 (m) are expressed in terms of trigonometric and gamma functions as 
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d 



1 



Vg( g 2 + l)r(f) 2 sinh(^/2) 



r U + f) [sm(7r(l + 2(? )/4) 2 -sm(^(l-i g )/4) 2 ] (508) 
, -2 + cos(tt(5 + zg)/4) 2 + cos(tt(3 + iq)/4) 2 
[cos(7r(5 + ig)/4) 2 - cos(^(3 + iq)/A) 2 \ 

Using the scalar base l m (x) the conjugate variables Q and P, describing our dynamical system, can be expanded 
in terms of traveling waves according to 

1 I' d 3 m 

Q(x,r,) = -j=J -^[Q* m (r,)l* m (x) +Q m (r ? )Ux)], (509) 



and 

d 3 m 



, . 1 f d 3 



3/2 



[P*( ?? )4(x)+P m ( ?7 )/ ro (x)]. (510) 



3. Quantum Treatment of the Perturbations 

The conjugate variables Q(x) and P(x) can be quantized following standard procedures, by transforming the mode 
functions into operators (cf. Birrcl, 1982) 

1 ( 511 ) 

P(x ' ^ = 72 f jiffi^PXA*) + d+p m i m ( x )]. 

The classical canonical variables are replaced by operators P and Q which satisfy the commutation relations [Q a , Pp] = 
ihS a p and [Q a , Qp] = [Pa, Pp] = 0, and through this we can define the operators of creation and annihilation as 

at = J^-(Qa(t)TP a ), (512) 

and from these operators the new commutation relations [a~,a^~] = 5 a p and [d~,dj] = [a+,a^] = 0. 

Through the expansion (509), we can write two different sets of functions corresponding to in-going and out-going 
waves, respectively given by 

1 f d 3 m 

Q m (x, V) = ^J ^pyj [SmCCW + a+v m l m ( X )} (513) 

and 

1 f d 3 m - 

Q OMt (x, rj) = -j=] -_[&-«^(x) + b+u m l m (x)] (514) 
The operators a m and b m are related by Bogolyubov coefficients a m and /3 m , 

tm=am«m+^*m A-m, (515) 

satisfying the normalization condition \a m \ 2 — \/3 m \ 2 = 1- The two sets of isotropic mode functions u m and v m form 
a base of in and out functions. We can write this base, according with (506) and (507), as follows: 
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1 



(516) 




uM =di(g)e^ + d 2 (g)e iOT . 



(517) 



The "in" vacuum is described by the standard Minkowski mode function and the mode corresponding to "out 
vacuum can be written as a linear combination of the basis v q : 



The existence of singularities appears to be a property inherent to most of the physically relevant solutions of 
Einstein equations, in particular to all known up-to-date black hole and conventional cosmological solutions (see 
Hawking and Ellis, 1977). In the case of black holes, to avoid the singularity some models have been proposed — cf. 
Bardccn (1968), Barrabes and Frolov (1996), Cabo and Ayon-Beato (1997) and, Mars et al. (1996). These models 
nonetheless are not exact solutions of Einstein equations since there are no physical sources associated to them. Many 
attempts try to solve this problem by modifying general relativity, for instance, Cvetic (1993), Tscytlin (1995), and 
Home and Horowitz (1992). More recently it has been shown that in the framework of standard general relativity it is 
possible to find spherically symmetric singularity-free solutions of the Einstein field equations that describe a regular 
black hole. The source of these solutions are generated by suitable nonlinear vector field Lagrangians, which in the 
weak field approximation become the usual linear Maxwell theory demonstrated by Ayon-Beato and Garcia (1998, 
1999a-b). Similarly in Cosmology many non-singular cosmological models with bounce were constructed where the 
energy conditions or the validity of Einstein gravity were violated (see Sec. [II D 2] for details). 

In 2002, de Lorenci et al. investigated a cosmological model with a source produced by a nonlinear generalization 
of electrodynamics and succeeded to obtain a regular cosmological model. The Lagrangian of such model is a function 
of the field invariants up to second order — Eq. (85). This modification is expected to be relevant when the fields 
reach large values, as is the case in the primeval era of our universe. The model is in the framework of the Einstein 
field equations and the bounce is possible because the singularity theorems are circumvented by the appearance 
of a negative pressure (although the energy density is positive definite). Recently some papers started a detailed 
investigation of the transition from contraction to expansion in the bounce of several models (see Cartier et al., 2003). 
In particular, in Einstein general relativity, models with stress-energy sources constituted by a collection of perfect 
fluids and FLRW like geometry were examined by Peter and Pinto-Neto (2001). The claim in that paper is that a 
generic result about the behavior of scalar adiabatic perturbations was obtained. The result is the following: scalar 
adiabatic perturbations can grow without limit in two situations represented by the points where the scale factor 
attains its minimum value and where p + p = 0. The first point corresponds to the moment in which the Universe 
passes through the bounce; the second corresponds to the transition from the region where the Null Energy Condition 
(NEC) is violated to the region where it is not. However, these instabilities are not an intrinsic property of generic 
models with bounce but a consequence of the existence of a divergence already appearing in the background solution, 
when the source is described as a perfect fluid. We show a specific example of a model with bounce, generated by a 
source representing two non-interacting perfect fluids, that has regular perturbations. 



We set the fundamental line element as being an FLRW metric. According to the definitions presented in Sec. 
[II D 2] , we rewrite the Einstein equations and the equation of energy conservation written for this metric below 



u q = a q v q - P* q v* q . 



(518) 



J. Nonlinear Electrodynamics: Scalar Perturbations 



1. 



The Model 




(519) 



'a ( a\ 2 e 
2-+ - + — +p 7 =0, 
a \a J 



(520) 
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p 7 + 3(p 7 +P-y)^ = 0. (521) 

Inserting Eqs. (87) and (88) in (521), it yields Eq. (89) for the magnetic field, where H is an arbitrary constant. 
With this result the equation (519) can be integrated. For the case e = the solution is given by Eq. (94). 

The interpretation of the source as a one component perfect fluid in an adiabatic regime has some difficulties (see 
Peter and Pinto-Neto (2001). The sound velocity of the fluid in this case is given by 

d Pi \ = h.= p -y (522) 

dpjj p 7 6»(p 7 +p 7 )' 

This expression, involving only the background, is not defined at those points where the energy density attains 
an extremum given by 6 = or p 7 + p 7 = 0. In terms of the cosmological time, they are determined by t = 
and He = 12a/ kc 2 . These points are well-behaved regular points of the geometry indicating that the description of 
the source is not appropriate. This difficulty can be circumvented if one adopts another description for the source 
of the model. This can be achieved by separating the part of the source related to Maxwell dynamics from the 
additional non-linear term dependent on a on the Lagrangian. By doing this, the source automatically splits in two 
noninteracting perfect fluids: 



where 



and 



Tpv — + Tp V , (523) 

= {pi +Pi)V fi V„ - p igfiU , (524) 
T% = (P2 + p 2 )V»V„ - p 2 g^, (525) 



Pi = \h\ (526) 



Pi = \h\ (527) 



p 2 = -4aH 4 , (528) 
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P2 = -y«i? 4 . (529) 

Using the above decomposition, it follows that each one of the two components of the fluid independently satisfies 
equation (521). This indicates that the source can be described by two non-interacting perfect fluids with equation 
of states pi = pi/3 and P2 — 5/32/3- The equation of state for the second fluid should be understood only formally as 
a mathematical device to allow for a fluid description. 
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2. Gauge Invariant Treatment of Perturbation 

The source of the background geometry is represented by two-fluids, each one having an independent equation of 
state relating the pressure and the energy density. Following the standard procedure we consider arbitrary pertur- 
bations that preserve each equation of state. Thus the general form of the perturbed energy-momentum tensor is 
written as 

57%, = (1 + WfftiVpVv) - XiS( Pi g^). (530) 

The background geometry is conformally flat. Thus any perturbation of the Weyl tensor is a true perturbation of 
the gravitational field. It is convenient to represent the Weyl tensor 1V Q ^„ in terms of its corresponding electric 
and magnetic parts, because these variables have the advantage that since they are null in the background, their 
perturbations are gauge invariant quantities (Hawking, 1966). 

The equations of motion for the first order perturbations arc linear so it is useful to develop all perturbed quantities 
in the spherical harmonics basis. Here, we will limit our analysis to perturbations represented in the scalar base. This 
basis, we also take from Sec. [HID]. 

In the case of scalar perturbations the fundamental set of equations determining the dynamics of the perturbations 
are: 

(SE^)'h, a hJ + {5E%yh«h/ + 6(SEf +SEf) = + Pl )Saf - \{p 2 +p 2 ) 6a f , (531) 

{SE 1 ^ + SEl^h^h^ = ±(6 Pl + 5p 2 ), a h^ - l - Pl 5v\ - ±p 2 6vl (532) 

(Sa^Y + (50%,)' + §V(<K + <K);a - UKff + S <p)^ a K) fi + ¥{5a% + K») = , s 

(533) 

= ~^ E % - 5E* V , 

-XiS(p,ph p lt )+(l + X 1 )pSal = 0, (534) 

- A 2 <5 (p,fih p „) + (1 + \ 2 )pSal = 0. (535) 

The acceleration a M , the expansion 9 and the shear cr M „ that appear in the above equations are parts of the irreducible 
components of the covariant derivative of the velocity field. 

The expansion of the perturbations in terms of the spherical harmonic basis is 33 

Sp = N(t)Q, (536) 

5V = V{t)h^ a Q. ai (537) 

5a" = Vh> MX Q, a , (538) 

5E" V = E(t)P^, (539) 

So*" = E^P"", (540) 



Note that once we are dealing with a linear process each mode can be analyzed separately. 
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3. Dynamics 

After presenting the necessary formalism we shall start the analysis of the perturbations of the bouncing cosmological 
model displayed in a previous section. Using the above expansion into the equations (531)-(535) we obtain: 



E x +E 2 



6e + fc 2 



(JVi + N 2 -p x V x -p 2 V 2 ), 



(541) 



Ei + E 2 + -e(Ei+E 2 ) = 
o 



Plhl - P2^2, 



(542) 



S x + S 2 - V x - V 2 = -Ei - E 2 



(543) 



- Ax (JVi - piVi) + (1 + Ai) Pl Vi = 0, 



(544) 



- A 2 (N 2 - P2 V 2 ) + (1 + A 2 ) p 2 V 2 = 0. 
These equations can be rewritten in a more convenient way as: 



(545) 



(546) 



(547) 



Si + -^ = -i(l + Ai)piEi, 



(548) 



S 2 + ^^2 = -^(l + A 2 )p 2 E 2 , 



(549) 



The whole set of scalar perturbations can be expressed in terms of the two basic variables: Ei and Ej. The corre- 
sponding equations can be decoupled. The result in terms of variables Ei is the following: 



Ei H — -OEi 



2 + 3A. i/l2 [2 \ 1 .., ol . (3e + fc 2 )A 4 



= 0. 



(550) 



Note that there is no sum in indices and j 7^ t in this expression. In our case the values of A^ are A^ = f ). In 



the first case the equation for the variable E\ became: 



.3' 3j 



Ei + - 8E\ 



- n2 ^ 
-0 -Pl-P2-^ 



E (l) = 0. 



(551) 



We should analyze the behavior of these perturbations in the neighborhood of the points where the energy density 
attains an extremum. This means not only the bouncing point but also the point in which p + p vanishes. Let us 
start by examining at the bouncing point t = 0. 

If we consider up to second order terms in the perturbation theory, the expansion of the equation of E\ in the 
neighborhood of the bouncing is given by: 
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E\ + a c tA + (b + M 2 )#i = 0. 
The constant a c and the parameters b and b\ arc defined as follows 



(552) 



(553) 



V6H t c 



and 



Defining a new variable f as 



1 2& • 



(554) 



(555) 



/(t) = J B 1 (t) exp 



c + 2 



(556) 



and doing a coordinate transformation for time as indicated bellow 



we obtain the following confluent hypergeometric equation (cf. Abramowitz and Stegun, 1974) 



(557) 



where 



e/+(l/2-0/ + e/ = 0, 



i(b-a c /2) 1 
4(6i - a2/4)V2 ~ 2" 



(558) 



(559) 



The solution of this equation is given by: 



/(*) = foM 



d,l/2,- l ( V / 4&i-a2)| 



(560) 



where / Q is an arbitrary constant and Ai(d, 1/2, £) is a confluent hypergeometric function. The confluent hypergeo- 
metric function is well behaved in this neighborhood and so is the perturbation E± (t) given by: 



f M(d, 1/2, -i(^^ c )L exp + iJb! - ^ < 



4 2 



(561) 



The perturbation £2 at this same neighborhood, after the same procedure we did before, results in the following 
equation: 



E 2 + a c tE 2 + (b + b 1 t 2 )E 2 =0. 



(562) 



This is the same equation we obtained for Ei, they differ only by the values of the parameters a Cl b and b\ that in 
this case are: 
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2tr 



(563) 



and 



Then the solution in this case is: 



E 2 (t) = K 



.= 4-5- m2 



2tl VGH Q t c 



5m 2 



f M [d, 1/2, -i(^^I) -) exp- ^ - - y fti - -f \t 2 



(564) 



(565) 



(566) 



Again the confluent hypergeometric function is well behaved in this neighborhood and so is the perturbation E^t). 
At the neighborhood of the point t = t c the equation for the perturbation E\ is given by 



E x + a c E x + (b + bit) E x = 0, 
where the parameters a c , b and &i in this case are given by 



(567) 



it' 



(568) 



and 



3 V3m 2 

= ~^r 6H Q t c 



, a/3 / 

b\ = —7T 



At 2 c \3H 2t c 



(569) 



(570) 



We would like to remark that this equation is different from the equations (552) and (562) obtained in the neigh- 
borhood of t = 0. We proceed doing the following variable transformation: 



Ex(t) = w(t) exp 



a r t 



The differential equation for this new variable is 



+ (6-(a c /2) 2 + M)w = 0. 



(571) 



(572) 



The solution for this equation is 



,{t) 



wq AiryAi — 



b- {a c /2) 2 + bxt 

b 2/3 



(573) 



The AiryAi are regular well behaved functions in this neighborhood and also the perturbations Ex ■ Finally, we look 
for the equation of E 2 at the neighborhood of t = to, and it becomes 
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E 2 + a c E 2 + (b + bit) E 2 = 0, 



(574) 



where the parameters a c , b and b\ in this case are given by 



9 



(575) 





) 



(576) 



and 



bi 



5V3 ( 1 



in 



2 



) 



(577) 



2t% \t c 



6H 



This equation differs from Eq.(567) only by the numerical values of the parameters a c , b and b\, so we obtain the 
same regular solution 



Summarizing, recently there has been renewed interest on nonsingular cosmology. As a direct consequence of this, 
some authors have argued against these models based on instability reasons. Peter and Pinto-Neto (2001) have argued 
that a rather general analysis shows that there are instabilities associated to some special points of the geometrical 
configuration. They correspond to the points of bouncing of the model and maxima of the energy density, where the 
description of the matter content in terms of a single perfect fluid does not apply. In the present paper we have shown, 
by a direct analysis of a specific nonsingular universe, that the result claimed in the quoted paper does not apply to 
our model. We took the example from a paper, De Lorenci et al. (2002), in which the avoidance of the singularity 
comes from a nonlinear electrodynamic theory. We used the quasi Maxwcllian equations of motion — cf. Novello et 
al., 1995a-b, 1996 — in order to undertake the analysis of the perturbed set of Einstein equations of motion. We 
showed that in the neighborhood of the special points in which a change of regime occurs, all independent perturbed 
quantities are well behaved. Consequently the model does not present any difficulty concerning its instability. This 
paves the way to investigate models with bounce in more detail and to consider them as good candidates to describe 
the evolution of the Universe. 



IV. ON THE ROLE OF INITIAL CONDITIONS TO THE EQUIVALENCE THEOREM 

In this section we briefly discuss the question that appears explicitly in perturbation theory: how do Einstein 
equations determine the Weyl tensor if they contain information only of the traces of the curvature tensor? 

Einstein equations relate the energy-momentum tensor with the traces of the curvature tensor (local quantities), 
leaving the remaining components of curvature tensor, which correspond to the Weyl tensor (nonlocal quantity), un- 
determined. However, this indetermination is apparent, because Bianchi's identities relate the traces of the curvature 
tensor with the Weyl tensor through equation (9). Nevertheless, remark that this relation involves partial derivative. 
Thus, substituting the Einstein equations into (9), we obtain the set of equations that involves the energy-momentum 
and Weyl tensors, which leaves the traces of curvature out — see equation (10). These equations arc direct consequence 
of Einstein equations, but they are not equivalent in principle. In order to get such equivalence, it is necessary to 
impose an appropriate initial conditions — see details about Lichnerowicz's theorems in [84]. 

From the mathematical point of view, the problem is automatically solved, however it does not occurs from the 
physical viewpoint, because it involves the determination of initial conditions from empirical data on a Cauchy surface. 
The empirical data that we get they determine the initial conditions of the physical system that we want to describe 
taking into account the curvature and energy-momentum tensors. Actually, these empirical data are precisely the 
curvature tensor of space-time and the energy-momentum tensor. Once the curvature tensor has 20 independent 
components, then we face the question: how can we determine the Cauchy data specifying the metric components gij 




(578) 
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and first time derivatives gty g (i, j = 1, 2, 3)? The Einstein equations are differential equations for the metric potential 
they do not give any information about the curvature tensor as initial condition, and then, this situation imposes 
a serious ambiguity on the determination of boundary conditions. This ambiguity clearly appears in gravitational 
waves process and can be consistently shown even in simple cases of perturbation theory, where the gravitational 
theory is linearized. Indeed, considering an exact solution with a perfect fluid as source of curvature, the perturbed 
equations describing gravitational waves reduce to 

6R„ V = 0. (579) 

In the QM formalism, yields 



5W a W; V = --(p + p)5V K[a Vp ] . (580) 

The equations (580) are not precisely equivalent to (579) because they consider distinct initial conditions. Besides, 
Eq. (580) has a kind of "source" for the gravitational field (in this case, due to coupling terms with shear tensor) 
absent in the Lifshitz perturbation theory. To fade away this ambiguity with clear examples, we have to resort to two 
recent papers on tensorial perturbations of isotropic metrics using the QM formalism [67, 129]. In these works, the 
authors use null coordinates to analyze tensorial perturbations in Robertson- Walker metrics and they encountered a 
gravitational wave type of solution different from those one obtained through the Bardeen method of gauge-invariant 
perturbations [9, 49]. The definition of gravitational waves in [67, 129] leads to no restriction on the equation of state 
of the fluid on the isotropic background, opposite to what is obtained by [16, 152], where the speed of sound is equal to 
the speed of light. According to what we said before, this difference can be associated to different boundary conditions 
given in each approach, which are fundamental to consistently define gravitational waves in each case. We will not 
enter the details of their works, but we would like to emphasize the question: in a nonlinear regime of gravitational 
theory, how can we translate the information contained in the nonnull Weyl tensor in terms of the perturbed potential 

In other words, GR does not contemplate in its boundary conditions any kind of information about intrinsic 
accelerations of a given configuration on the manifold, which would lead to the development of a nonlocal description 
of the space-time as suggested by [19, 32, 88, 89]. This type of information shall explicitly be given to the QM 
formalism, because these equations have partial derivatives of the energy-momentum tensor and contain higher order 
derivatives of the metric tensor in the dynamics. 

V. CONCLUDING REMARKS 

During the whole history of the quasi-Maxwellian equations of general relativity, they have been just considered as 
an alternative (and, eventually more complicated) approach equivalent to the Einstein equations. However, even Prof. 
Jordan and collaborators knew the importance of this issue, as we quoted: With this paper we start an enterprise 
which, considering the present state of the scientific development of this field, is of some urgency. We want to collect 
and describe all those exact solutions which are dispersed in the literature and thus probably, in their completeness, 
known to very few authors only. [76] After all these years of their seminal work, we summarize here all those works 
which have been done in this area showing that this formalism is actually as powerful as the standard one, once the 
well-known solutions of GR are easily regained from QM-formalism. 

In the case of conformally flat universes, singular [118] or not [116], this approach is more convenient to treat all 
types of perturbations (scalar, vectorial and tensorial). It reduces the dynamics to a pair of equations in £ and E, 
which are variables physically meaningful, providing a dynamical planar system and a reparametrization of these 
variables allows to establish a gauge-invariant Hamiltonian treatment for the perturbation. In particular, some years 
ago, there has been a renewed interest on nonsingular cosmology and we review that these models are completely 
stable based on results obtained from the QM-approach. In the case of Schwarzschild metric [79] as well as the Kasner 
case [95], it is also possible to use the Stewart Lemma [143] for applying such method and get significant results. 

There are many areas of research on this topic. Some of them concern the development of deep analogies between 
QM-equations and the Maxwell theory [57, 58]. Another topic concerns modifications on the dynamics of general 
relativity suggested by Bianchi's identities for the Weyl tensor [103]. Most of these issues are still under development. 
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VI. APPENDIX: AN EXAMPLE OF FURTHER DEVELOPMENTS IN QM-FORMALISM 

This appendix presents an example in which the QM-formalism is modified and the analogy with the electromag- 
netism is strengthened. Modifications of this kind indicate possible manners to solve problems that appear in the 
general relativity theory in a more intuitive way. We present an enlightening example next. 

Electrodynamics is the paradigm of field theory. Its theoretical and experimental properties have been simulated 
and sought for in many other theories, and in particular in the analysis of gravitational phenomena. Many works have 
been done along this line and discussed the resemblance between electrodynamics and gravidynamics 34 . However, it 
seems possible to further improve this similarity, as we shall show. 

In this vein, we review here a modification of Einstein's theory of general relativity under certain special states of 
the geometry of the space-time. Since the original proposal of Einstein's geometrization of gravitational processes, 
many physicists have discussed alternative models of gravitation. The kind of theory we shall analyze here is given 
by means of the metric properties — represented by a symmetrical metric tensor gp, u {x) — and by two other functions, 
e{x) and which arc independent from the metric 35 but have intimate connection to the space-time. 

For pedagogical reasons, we think it is convenient to limit our considerations in this presentation to the case in 
which both s and fi are constants. The meaning one should attribute to these two constants is obtained by a direct 
analogy with the dielectric and permeability constants of a given medium in electrodynamics. 

We shall simplify the model by merely stating that e and fj, can be provisionally identified with the characteristics 
of certain states of tensions, in free space-time, due to an average procedure on (quantum) properties of gravitation. 

In other words, e and /i are interpreted as the result — in a macroscopic level — of some sort of averaging microscopic 
field fluctuation 36 . This is perhaps not difficult to assume if we can say exactly how the equations of motion of gravity 
phenomena must be modified by them, as we shall do later. We remark that we are not supporting this interpretation 
but merely suggesting it as a provisional sursis of the model. 

We shall describe gravitational interaction 37 by means of a fourth-rank tensor We shall set up its algebraic 

properties and give its dynamics. It is possible to separate this tensor, for an observer moving with four- velocity V^, 
into four second-order symmetric trace- free tensors E a p , B a p , D a p and H a /3 . The principal result is then obtained by 
showing that it is possible to select a class of observers with velocity in such a way as to have equations of motion 
for Q a f3fj,L> similar to the Maxwell equations for the Electrodynamics. That is, for E a p, D a p, B a p and H a p separated 
into two groups: one containing only E a p and B a p (and their derivatives) and the other containing only H a p and 
D a p (and their derivatives). These equations have the same formal structure of Maxwell's equations in a given general 
medium. So, we arrived at the conclusion that in the present theory there is a class of privileged observers in which 
gravitational field equations admit this simple separated form. Any other observer, which is in motion with respect 
to mixes the terms E a p, D a p, H a p and B a p into the equations. This situation could be thought of as defining a 
new type of ether. However, unlike the ether of the pre-Einstein epoch, our ether is not a substance, but it is only a 
preferred frame of observation. 

In the remainder of this presentation we discuss in some detail a very particular situation of these tensors, that 
is, the case in which they can be reduced to two tensors plus two constants: they arc e and /i. Then we show that 
Einstein's theory is obtained from this for a particular set of values of e and /j,, that is, e = [i = 1. It is in this sense 
that we can call this theory a generalization of Einstein's gravidynamics. 



This resemblance is far from being accepted by all physics community. Indeed, in the final session of the 1972 Copenhagen International 
Conference on Gravitation and Relativity, A. Trautman argued that perhaps many of the difficulties of gravitational theory may be due 
to the extension of this similarity to all aspects of both fields. 

This hypothesis is made here only for simplicity. It is an over-simplification under certain drastic situations, such as, for instance, for 
very strong gravitational fields. 

We call the reader's attention to the fact that we are, perhaps, in a situation similar to that experienced by Maxwell, a century ago. 
His theory described the electromagnetic fields in the interior of substances by means of the same type of fields in vacuum and by 
characterizing the distortion produced by the matter on the fields, as given by macroscopic quantities: the dielectric constant e^ax ant ^ 
the permeability ^Max 0*Y "Max" we mean Maxwell). It took many years before Lorentz - having at the disposal the atomic theory of 
matter - made Maxwell's theory rigorously understood be averaging properties of microscopic fields in a macroscopic scale. 
We remark that in the present paper we shall limit ourselves to the sourceless case, i.e. the so-called vacuum gravitational fields. A 
further generalization to include matter is straightforward and does not present difficulties. 
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A. The Q-Field 

1. Definitions 

Let us define in a four-dimensional Ricmannian manifold a fourth-rank tensor Q a p^v described by an observer 
in terms of four second-order tensors E a p, D a p, H a p and B a p. We set, by analogy with the irreducible decomposition 
of the Weyl tensor, that 

Q a p^ = V [a D 0] \»Vl + V [a E m ^ + 5feE% - W ^B*rl - rf vp °V P B a{a Vp X . (581) 
The tensors E a p, D a p, B a p and H a p, represented below by X a p, satisfy the following properties: 

X a a = 0, (582a) 
X aP V a = 0, (582b) 
X a p = Xp a . (582c) 

We can write D a p, E a p etc., in terms of Q a p^,u and projections on like, for instance 

D a[i = -Q safl pV e V» 

and so on. 

2. Algebraic properties 
From definition (581) of Q a p^u it follows directly the properties: 

Qa^r = -QpcT, (583) 

Q a p^ = -Q a p" M , (584) 

Q a p av = Ep v -Dp v , (585) 

Q a a = (586) 

3. Dynamics 

By analogy with Einstein's equations in vacuum we impose on Q a pnv the equation of motion 38 

Q a ^ v . v = (587) 

Now, we shall use the above properties to project the system of Eq. (587) parallel and orthogonal to the rest frame 
of a selected observer with 4— velocity from the whole class of V^. We impose that the congruence generated by 
£ M satisfy the properties. 

« M = +1, (588a) 
1 
2 

Oafi = \ h {* X h) e i\;e = 0, (588c) 

t l = l». v l v = 0. (588d) 

(588e) 



w a p = ^h [a x h ?] e l x -e = 0, (588b) 



38 Indeed, as we shall see, in the case in which E a R = D a p and B a p = H a p, Qafivy, can be identified to Weyl's tensor and Eqs. (10) 
reduces to Einstein's equation in the vacuum. 
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where /i M „ is the projector in the plane orthogonal to that is 

Ku = 9n» - V" (589) 

So, the congruence generated by t p is geodesic, irrotational, non-expanding and shear-free. The reason for selecting 
such a particular class of observers will appear clear later on. Then, Eq. (587) assumes the form: 

D aii;v h" v h a £ = , (590a) 

D^h a {a h$ + hf a ne) ypr t p H Ta , v = 0, (590b) 

B ail , v h^h a e = 0, (590c) 

B^h Ka h A) „ - h£ v x) upT e p E Ta] „ = , (590d) 

in which a parenthesis means symmetrization. 

This set of equations has a striking resemblance with Maxwell's macroscopic equations of electrodynamics. Indeed, 
we can formally understand the above set as being [76] 

V • D = , (591a) 
D - V x H = , (591b) 

V • B = , (591c) 
B + V x E = , (591d) 

where the symbol — > is put over D, E etc. only to represent its tensorial character; the V operator represents the 
generalizations of the usual well-known differential operators. 

So, we can understand the reason for selecting the above privileged set of observers, given by the tangential vector 
l p . Eq. (587) takes the form (590) only for this class of frame. Any other observer which is in motion with respect to 
£ M will mix into the equations of motion of the set of tensors (E a p, B a p) with the set of tensors (D a p,H a p). So, it is 
in this sense that there is a natural selection of observers, with respect to the equation of motion satisfied by Q a pp,v 



4- e and [i states of tension 



A particular class of states of space-time occurs in the case in which there is a linear function relating the tensors 
B a p with H a p and E a p with D a p by the intermediary of two constants, e and \i. 
We set 

B a \ = fiH aX , (592a) 
D aX = sH aX . (592b) 

If we put expressions (592) into definition (581) of Q a pnvi a straightforward calculation shows that it is possible to 
write Qafinv in terms of the Weyl tensor C a p^ , and its electric and magnetic parts £ a p and H a p , if we identify the 
tensor E a g with £ a p and H a p with H Q ^. Then we can write 

Q a p = C a ^ v + (e - l)i [a 4] £ v] + (1 - riv^n^ tk p , (593) 

where 

cy = u [a e^f] + s fe £ fi - w *t x n a ^ - Tf vp °i p u a{a i 0] (594) 

and, consequently, 

S a p = -C a ^ v l mu t , (595) 



(596) 
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The resulting equations of motion (13)-(14) turn into the set: 

£ a ^\ W h^h a £ = Q, (597a) 

e£ afi hf a h$ + hfc rfj; r i P U Ta \\ v = , (597b) 

n a ^ v h^h% = , (597c) 

iMa^ ~ h a [a V ^ T £ p £ Ta ^ = . (597d) 

By the same argument that guided us to Eqs. (590) we see from the above set that we can identify e as being the 
gravitational analogue of the dielectric constant of electrodynamics, and /z as being the permeability of space-time. 

Now, we recognize in Eqs. (597) Einstein's equations for the free gravitational field for the particular case in which 
e = [l = l. 39 

So, it seems natural to interpret Eqs. (597) for the general case (e, [i different from unity) as the equations for 
the gravitational field on states of space-time that are macroscopically characterized (in the sense discussed in the 
introduction) by the two constants e and fj, 

5. Conformal behavior of Q Q ^„ 

A conformal transformation of the metric g afi is given by the map 

g^uix) — ► g^vix) = n 2 {x)g tiv {x) 

<T{x) — ► g^(x) = n~ 2 (x)g^(x) (598) 
How we can set T) a g as independent from the conformal transformation: 

VaP — VaP 

Then, it is easy to see that the electric and magnetic parts of Weyl tensor remain unchanged, 

£fj,u Cfipv(j£^£ — £fiu ? (599a) 

Uya, = -\ V pa Pa C pavX M X = %fj, v , (599b) 

where we have used conformal transformation of the velocity as usual, 

¥ = rr 1 ^ (600) 

As a consequence of the above mapping, Qap^v behaves, under the conformal transformation, as the Weyl tensor, 

Qaj3 P u{x) = ^l 2 {x)Q a p p , l/ {x) . (601) 

B. Gravitational Energy in an e — State of Tension 

There have been many discussions, since Einstein's 1916 paper, concerning the definition of the energy of a given 
gravitational field. We do not intend to discuss this subject here but we shall limit ourselves to consider one reasonably 
successful suggestion of Bel [12] for the form of the energy-momentum tensor of gravitational radiation. 

The point of departure, see [12], comes from the similitude of the equation of motion of gravity and electrodynamics. 
He defines a fourth-rank tensor T ,Q ^ A " / given in terms of quadratic components of the field (identified with the Riemann 
tensor) and written in terms of the Weyl tensor C a/3,lu . 

Bel's super-energy tensor takes the form: 



It seems worthwhile to observe that this equivalence is complete only if we impose as initial date the set R^v = on a given space-like 
hypersurface. 
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where * is the dual operator. Note that the symmetry of the Weyl tensor (C*p = C a p* = *C a p llv ) does not hold 
for QatPfiw This is related to the lack of Q a /3 P u 7^ Qfiva/3 symmetry. Indeed, we have 



2 \- -/■ivc/jftu - 2 

This T a ^^ v tensor has properties very similar indeed to the Minkowski energy- momentum tensor of electrodynamics. 
The scalar constructed with and the tangent vector for instance, takes the form 



and 



Then wc have 



and gives the "energy" of the field 



where 



U {T) = T^^lJLpl^ (603) 



U (T) = l -(Z 2 +U 2 ) , (604) 



£ 2 = £ af} £ afi ■ (605a) 
H 2 = U af 3U afs . (605b) 

In the context of the present extended theory, for a space-time in the state e — fx of tension, we are led to modify 
rpotfitiv j n £ Qap^v d enn ed j n an analogous manner by 

Q ^ = l{ Q ^C? p % + Q^*Ce p %} . (606) 
Then, the energy Ur e ,^) a s viewed by an observer V 1 will be given by 

U { ^) = ffrt^ljtplplv = \ (e£ 2 + iM 2 ) (607) 

in complete analogy with the electrodynamics case in a general medium. 

We would like to make an additional remark by presenting two special properties of 8 Q ' 3 ^ 1/ 

® a pa„ = l(l-e)£ P<T Cp p ^ , (608a) 

e = e Q % = o . (608b) 

Property (608a) states that not all traces of B"^" are null for a general state of tension of space-time and that the 
non-null parts of the contracted tensor are independent of the "permeability" [i. The second property (608b) states 
that the scalar obtained by taking the trace of Q a P^ v twice is null, independent from the state of tension of the 
space-time. 

C. The Velocity of Propagation of Gravitational Disturbances in States of Tension 

In order to know the velocity of gravitational waves in e — fi states of space-time, let us perturb the set of equations 
(590). The perturbation will be represented by the map: 

£^u — > £ pv + 5£^u , (609a) 
— > U pv + 6H P „ (609b) 
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in which b£ llv , 5T-L pv are null quantities. Then, Eqs. (590) go into the perturbed set of equations. 



dSj.p « , (610a) 

eS£ atl + X - h^^ip SH T x: P « , (610b) 

SHj-13 m , (610c) 

fiSH-^ h\ a i h f^£ TX , p « , (610d) 



where the covariant derivative is taken in the background - and we limit ourselves to the linear terms of the pertur- 
bation. 

Now, let us specialize the background to be the flat Minkowski space-time 40 . In this case the covariant derivatives 
are the usual derivation and we can use commutative property to write: 

eS£ aP + 1 h x (am paT e a 5iirx, p « , (611) 
by taking the derivative of Eq. (610b) projected in the privileged direction . Multiplying Eq. (610d) by the factor 



2H ft? ■ 

we find 

\h\ a SH lv>a - ±- h\ a V0> °n £ T £ p h s {l riu^Se^,^ « . (612) 

Substituting Eq. (611) into (611) we finally find 

6£ ap - — V 2 5£ ap = , (613) 
en 

where V 2 is the Laplacian operator defined in the three-dimensional space orthogonal to 

In the same way an analogous wave equation can be obtained for T~L ap - From Eq. (613) we obtain the expected 

result: the velocity of propagation of gravitational waves in e, \i states of tension of space-time is equal to 1/^/e/Z. 
The set of privileged observers that we dealt with here may be enlarged by somehow weakening the defining 

conditions (see Eqs. (588)). This point deserves further investigation. 
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